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EDITORIAL 


A perspective on seismology in India 

Dr Vinod Gaur kindly asked me to be a guest editor for a special issue of the 
Proceedings of the Indian Academy of Sciences to be focussed on seismology in India. 
We selected a number of prospective topics and authors and I began negotiations for 
contributions. The result is the group of excellent papers in this volume. I made no 
elfort toward a comprehensive coverage of seismology, but rather attempted to 
choose topics illustrative of a range of problems related to India. Thus many 
important topics (e.g. induced seismicity) are not represented. Some invited authors 
were unable to contribute, and probably others should have been invited to represent 
their topics, but the press of time and schedule, mostly on my part, made this 
impossible. I apologize to all those who could have made valuable contributions to 
this volume. The paper by Royer et al on the tectonic fabric of the Indian Ocean is not 
strictly seismological, but provides an important background for seismological 
interpretation, especially with regard to the structure in the Bay of Bengal. Many may 
feel that the authorship is too heavily weighted with scientists from outside India, and 
thus the overall perspective tends more to that of an outsider, and I would not 
disagree with this evaluation, but I am nevertheless hopeful that the volume will be 
stimulating to seismologists both inside and outside India. I would like to thank 
Dr Kailash Khattri for reviewing some of the contributions for me. 

I began my career in seismology during the IGY (International Geophysical Year), 
an international programme which led to the installation of modern seismographs 
worldwide, and provided the data for my thesis. Thus I am very aware of the 
importance of international co-operative projects in seismology. The continuing 
benefit of such projects is illustrated by the strong motion data now being collected in 
India as a result of the Hawaii International Workshop on Strong Motion Arrays 
(1978), which stimulated arrays in India, Taiwan, Mexico, China, the U.S., and Turkey 
(see the article on strong motion seismology in this volume). The importance of such 
arrays is graphically emphasized by the recent damaging earthquake in the 
India-Ncpal border region. Recognition of the importance of international co¬ 
operative projects is most recently illustrated by the plans for the International 
Decade of Natural Hazard Reduction, including an important meeting to be held in 
New Delhi in the summer of 1989. Of course national efforts in research will usually 
take priority over international efforts, however, most nations are realizing that inter¬ 
national co-operative projects are crucial to their own national development. 

The problems inherent in developing international co-operative projects illustrate 
some basic problems in research and I would like to share some of my views on these, 
based on my years of fruitful co-operative research in Mexico and my more recent 
visits to India. 

As a scientific discipline seismology has developed from a young science to a mature 
science. The characteristics of the types of research programmes have changed, and 
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ana purely naiionai researcn. Many oi us can rememoer wnen ulscyvcllc^ lu 
seismology were simple and inexpensive, and contrast that situation with the current 
situation where much of the funding for seismology is for large multi-institutional 
projects. 

When seismology was young, major discoveries, causing qualitatively great changes 
in our understanding, were made relatively simply by a few dedicated “basic” 
researchers and students, with minimal funding support. Research was carried out 
with individually constructed equipment and simple calculations using elementary 
first-order theory inspired by physical intuition. The individual scientist struggling to 
advance his own career, automatically helped other scientists advance their careers, 
and promoted the rapid advancement of science. Examples of discoveries during this 
stage include discoveries of the Mohorovicic discontinuity, the earth’s core and 
general structure, surface wave dispersion, the double-couple earthquake mechanism, 
free oscillations of the earth, and of course, plate tectonics. 

Seismology has now matured in many respects. Emphasis is on relatively small 
improvements in quantitative knowledge, and eliminating non-uniqueness in 
sophisticated theoretical models. These are brought about only with great difficulty 
by teams of scientists and technicians, many with an applied research and 
management orientation. There is a relatively great expenditure of money for small 
improvements in accuracy, with the research direction dictated by societal goals and 
relevance (e.g. nuclear test detection, and earthquake hazard reduction). The tools are 
standardized accurately calibrated commercial equipment, sophisticated higher order 
theoretical calculations requiring great mathematical skill and large computers, with 
relatively less reliance on intuition. If developing countries are to fully contribute they 
need to participate in multidisciplinary, multinational, mission-oriented projects, and 
need to be willing to invest the relatively large amounts of time and money required. 
Examples of such projects include large seismic arrays such as NORSAR and the 
Gauribidanur array, the worldwide seismograph network, strong motion arrays, 
multidisciplinary and multi-institutional ocean-going surveys, IRIS (Incorporated 
Research Institutions for Seismology) and PASSCAL (Program for Array Seismic 
Studies of the Continental Lithosphere). These projects usually result in papers with 
several co-authors. 

As we are all aware, at this stage politics may play an important role. Large 
international or bi-national projects can impact a nation economically and militarily. 
Research can lead to unexpected results, and politicians are reluctant to commit 
national efforts and funds to projects which they do not understand and fear may 
affect them adversely. As a consequence, important research which could help a nation 
is not carried out. An excellent example might be the seismic hazard research in the 
Himalayas, which could obviously benefit millions of people, but which unfortunately 
is hampered by border security concerns. Because of the uniqueness of the Himalayan 
region, understanding the earth structure, tectonics, earthquake mechanism and wave 
propagation in that region could be of great importance to the international 
community, and in particular to the objectives of the International Decade of Natural 
Hazard Reduction. 

There are strong incentives for government agencies to fund international 
co-operative projects e.g. national, economical and political advancement, national 
hazard reduction, increasing the national educational level, and, last but not least. 
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development of basic research capabilities. In the long run basic research will help all 
of the above objectives, but in the short run is often considered too costly a drain on 
resources for attacking more immediate social problems. 

The career advancement of individual investigators, students, and technicians is one 
of the motivating factors behind co-operative research. Priorities must be balanced 
between, on the one hand, publishing results as soon as possible, and on the other 
hand allowing students and other apprentices to take the initiative on important 
problems in order to learn more effectively. There is the obvious scientific incentive of 
publishing as soon as possible. On the other hand, “hot” research projects are usually 
best for teaching students because they help develop enthusiasm and competence 
which will last into the future. Students and apprentices are less experienced and often 
take longer to work up data and finish research. A wise balance is required in these 
matters in order to maximize the short and long-term benefits to both the nation and 
the international community. 

Of course military and economic interests often interfere with co-operative 
research. I will not discuss this obvious and much analysed problem except to observe 
that the United Nations has been much less effective in helping the situation than I 
and many others had hoped 20 years ago. We may hope that time will slowly improve 
the problem. 

The problems discussed above have general relevance to India. As far as seismology 
in particular is concerned, the two problems I have had some experience with are 
study of the crustal structure in the Bay of Bengal, and attempting to share 
co-operative seismic source mechanism and hazard projects in the Himalayas. Both 
have been effected by such difficulties. Both the Himalayan tectonic belt and the Bengal 
sedimentary basin are unique and the largest examples of their type on our planet, and 
thus attract high interest. They are of obvious military and economic sensitivity. 
Naturally Indian researchers would like to have the funding support to attack these 
problems themselves. And, because of their uniqueness, they are also of great 
importance to the international community and consequently there is strong incentive 
to attempt to initiate co-operative field and sea-going operations. One would hope 
that some way would be found to overcome the difficulties involved with such 
co-operative research in such a way that all parties would benefit. 

A specific problem in seismology is the lack of facilities and incentives to make data, 
especially early data, available to interested researchers. Again two incentives are at 
work. On the one hand there is the incentive to get the research done as soon as 
possible because it would provide answers to scientific problems important to the 
nation. From this point of view one would like to have the data made as widely 
available as possible, to both national and international researchers. On the other 
hand this policy would entail many of the real or perceived national and personal 
risks touched on above. To mention one, if other researchers use the data to make all 
of the exciting discoveries, this will deprive the students and personnel of the 
institution in charge of the data of the opportunity for making the discoveries 
themselves, and thus advancing their own careers. A wise balance would be to make 
sure that the institution in charge of the data has both an incentive for co-operation 
and adequate facilities and funding to co-operate with others on an equal basis. 

As far as academic preparation for research is concerned, I feel that all the 
mechanisms for the high level training, required are potentially in place in India. What 


necessary to be fully effective in the modern age of complex multidisciplinary rest 
The isolated “ivory tower” scientist is obviously still important, however there 
reason that such a person cannot be associated with an applied research group s< 
valuable cross-fertilization of ideas can occur. This cross-fertilization is difficult 
scientist is the only one of his field at a spatially isolated university, even i 
university is of high academic quality. 

Government research groups are obviously critical in modern research, but 
they could be more effective if there was a closer tie with nearby academic institu 
Graduate student research and training, along with the natural energy 
enthusiasm of students is important to stimulating fundamental research at 
centers. 


The future 

As an old saying goes; “The future is upon us”. This is certainly true in seism 
with the advent of exploding computing power, rapid data transmission capab 
greatly improved instruments, modern techniques of analysis of digital dat£ 
rapidly improving transportation and logistical facilities. I believe the way i 
future is toward standardized data collection techniques. A critical loss of effi( 
will result without this standardization. A recent publication of the U.S. Na 
Academy Press (1988) “Geophysical data; Policy issues”, graphically points o 
fact that a tremendous amount of valuable data, collected at considerable cost, 
actually available because of the lack of standardization of formats. This is nati 
the early stages of development of a scientific field, but should be corrected as 
develops. The report recommends that a series of guidelines and policies be ad 
to correct the situation and to make future data more easily accessible and usei 
for many of the issues discussed above, international co-operation will be crit 
this regard. Formats should be standardized as soon as possible so that scient 
developing countries will have easy access to the modern data, and to moderi 
analysis techniques. Data collected at considerable public expense cannot be “1 
up” because only one or two scientists or technicians have the facilities and kno^ 
to use it. especially when this data may relate to such issues as seismic safei 
national economic development. 

India currently has the capability of carrying out state-of-the-art resear 
determination of earth structure, and the primary challenge of the future v 
setting priorities in funding, based on the economic benefit of such si 
Determination of earth structure can be of great importance to mineral explo 
and seismic hazard reduction. Passive continental margins, the Bengal Fan, ai 
Indus Fan have great economic potential, and determination of regional anc 
crustal structure using seismic means may aid local exploration efforts. There n 
some benefit to international co-operative studies in the deeper water parts oi 
structures. 

It is obvious that India, with such a high population concentration i 
Himalayan seismic region, will have to devote a major effort to seismic li 
research in the next decade. The major thrust of this will have to be toward; (i) de 
the actual seismic hazard in India (e.g. through seismicity studies, geologic stuc 
recent fault motion, refined estimates of probability, measurement of Fourier si 


and response spectra of earthquakes, determination of local ground motion 
amplification), and (ii) clarifying the various options for confronting these hazards 
(specifically involving cost-benefit analysis). The importance of such studies is 
highlighted by the controversy over the seismic hazard to Tehri Dam, a construction 
project of tremendous potential benefit to India, and at the same time of tremendous 
potential devastation if it should fail during the great earthquake expected in the 
region. In spite of such tremendous importance to India, knowledge of the tectonics 
and seismology of the region is only primitive, but could be greatly improved within a 
few years with modern instrumentation and geological studies. 

Some specific critical decisions relating to seismology in India will probably need to 
be made in the next few years, or perhaps even in the next few months: 

(i) Should a modern standardized national seismic system, both for strong and weak 
motion recording, be established? 

(ii) How much effort should be made to integrate seismic data systems with those of 
the international community? 

(iii) Should data transmission be primarily by satellite, with event recording, or 
should part of the national array have continuous recording, perhaps with microwave 
transmission? 

(iv) How much emphasis should be placed on PC computer facilities for individual 
scientists, versus emphasis on large centralized facilities for data analysis and 
theoretical computations? 

(v) How much emphasis should be placed on small research grants to individual 
investigators, versus support for large multidisciplinary, multi-institutional projects? 

(vi) How much effort and expense should be put into earthquake resistant building 
design? 

(vii) How much effort should be put into short-term earthquake prediction and 
earthquake warning systems? 

Most of these issues are common to many other nations. In order to expedite the 
decision-making process, I would recommend that a high-level panel or commission 
be set up to make recommendations to the government on the above issues. 

In conclusion, I feel that seismology in India is at the threshold of a most interesting 
and challenging period. The rewards for correct policy decisions can be great, but the 
consequences of error can be disastrous. India has the national capability to 
successfully attack these problems by itself, but the most cost-effective strategy would 
probably be to participate as fully as possible in future international co-operative 
projects such as the International Decade of Natural Hazard Reduction. 

Department of Geological Sciences 
Mackay School of Mines 
University of Reno 
Reno, Nevada 89557, USA 


James N. Brune 
Guest Editor 




A preliminary tectonic fabric chart of the Indian Ocean 
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Abstract. We present a preliminary tectonic chart of the Indian Ocean based on a joint 
compilation of bathymetric data, magnetic anomaly data and Geosat altimetry data. Satellite 
altimeters such as Geosat map the topography of the equipotential sea surface or marine 
geoid. Our interpretation of the GEOSAT data is based on an analysis of the first derivative 
of the geoid profiles (i.e. deflection of the vertical profiles). Because of the high correlation 
between the vertical deflection (at wavelength <200 km) and the seafloor topography, the 
Geosat profiles can be used to delineate accurately numerous tectonic features of the ocean 
floor such as fracture zones, seamounts and spreading ridges. The lineations in the Geosat 
data are compared with bathymetric data and combined with magnetic anomaly identifica¬ 
tions to produce a tectonic fabric chart of the Indian Ocean floor. 

Keywords. Indian Ocean; GEOSAT altimeter; deflection of the vertical; tectonic fabric 
chart. 


1. Introduction 

During the past few years, as a result of development of remote sensing and data 
imaging techniques, better resolution and coverage have permitted the upgrading of 
the tectonic picture of the ocean floor. In the Indian Ocean, one may examine the 
development of this improved resolution by comparing the successive tectonic charts 
produced by Heezen and Tharp (1965) and Udintsev etal (1975), and the bathymetric 
charts of the GEBCO series (1975-1982) with the recent chart of the gravity field by 
Haxby (1985). Satellite altimeters (GEOS-3, SEASAT, GEOSAT) map the topography 
of the equipotential sea surface (marine geoid). Due to the uniform coverage of the 
satellite ground tracks and the excellent correlation between the short wavelengths 
(20-200 km) of the geoid signal with the topographic gradient of the ocean floor, the 
satellite altimeter data bring a wealth of rigorously comparable information bearing 
on the morphology of the seafloor, especially in the remote inhospitable and hence 
poorly charted southern oceans. There, in particular, previously unrecognized tectonic 
elements, specifically fracture zones and seamounts, have been delineated or verified 
(e.g. Lazarewicz and Schwank 1982; Sand well 1984; Haxby 1985; Sand well and 
McAdoo 1988). 

Recently, Gahagan etal (1988) have used the anomalies and lineations in plots of 
the deflection of the vertical from SEASAT data to portray a preliminary fabric chart 
of the ocean floor. In this paper, we expand upon this approach by presenting first, 
plots of the actual deflection of the vertical data along tracks, second, identifying 
lineations in these plots and finally, comparing these lineations with bathymetric 
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features and magnetic anomaly identifications in order to construct a prelim 
tectonic fabric map of the ocean floor. This work is based primarily on the GEC 
data which we combine with the SEASAT data set. This analysis reveals new fea 
on the seafloor at the southerly latitudes that are related to the early phase cvof 
of the Indian Ocean. Recently-refined portrayals of tectonic fabrics of the sea 
provide exceptionally good constraints for plate tectonic reconstructions far be 
those available to the pioneer modellers of the evolution of the Indian Ocean 
Fisher etal \971, McKenzie and Sdater 1971; Sdater and Fisher 1974; Nortoi 
Sdater 1979; Partriat 1987). In fact, problems with previous reconstruction m 
have generally arisen in areas where either the fracture zone lineations or the maj] 
lineations were poorly established through lack of shipborne coverage. Hcnct 
joint compilation of magnetic anomalies and fracture zone lineations from dc 
spaced altimetry traverses has proved to be a very powerful approach as it allo' 
to revise and improve our understanding of the tectonic history of the ocean 
(e.g. Royer etal 1988; Royer and Sandwell 1989; Mayes etal 1989). 


2. A tectonic history of the Indian Ocean 

The Indian Ocean floor (figure 1) is characterized by a system of three active sprc; 
ridges that now separates four major fragments of the former supcrcont 
Gondwana; Africa, India, Australia and Antarctica. The northern branch of this 
system rises in the Gulf of Aden, separating Africa from Arabia, and continues i 
the Carlsberg Ridge and the Central Indian Ridge which separates Africa from 1 
At 25°S, 70°E, the Central Indian Ridge intersects the two other branchef 
Southwest Indian Ridge that extends between Africa and Antarctica toward 
South Atlantic, and the Southeast Indian Ridge that extends towards the South P 
between Antarctica and India, and Antarctica and Australia. Numerous ridge; 
plateaus of intermediate depth, which may be taken to characterize the Indian 0 
appear on either sides of the spreading ridges. The size, number and distributi 
these elevations have raised many questions about their nature and origin as w 
their role in the continuing development of the Indian Ocean. From east to 
they are the South Tasman Rise, the plateaus off western Australia (Natur; 
Cuvier, Wallaby, Exmouth), Broken Ridge and the Kerguelen Plateau, Ninel 
Ridge, Chagos-Laccadive Ridge and the Mascarene Plateau, the Chain and M 
Ridges, Madagascar Plateau, Del Cano Rise, Crozet Plateau, Conrad Rise, Gun 
Ridge, Mozambique Ridge, Astrid Ridge, Agulhas Plateau and finally Maud 
(figure 1). 

According to interpretations of the seafloor magnetic anomaly pattern as ider 
in the Indian Ocean, the major basins lying between the spreading ridgef 
continental margins and the submarine ridges evolved during three main pe 
From the breakup of Gondwana in the Late Jurassic to the mid-Cretaceous, / 
separated from Madagascar and Antarctica, creating respectively the western S( 
Basin (Segoufin and Patriat 1980; Rabinowitz et al 1983; Coffin and Rabinowitz 

Z"' 10 00\ _____1 * ' o* 1 r\i^ f\ • t 
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1978; Larson et al 1979; Veevers et al 1985). From the mid-Cretaceous to the Middle 
Eocene, the major motion resulted in the rapid northward drift of India towards 
Asia. Most of the floor of the Indian Ocean was created during this phase: the basins 
between Africa-Madagascar and Antarctica (Bergh and Norton 1976; Patriat 1979, 
LaBrecque and Hayes 1979; Sclater et al 1981; Fisher and Sclater 1983), the 
Madagascar and the Mascarene Basins, the eastern Somali Basin and the Arabian 
Sea between Africa-Madagascar and India (McKenzie and Sclater 1971; Whilmarsh 
1974; Schlich 1975,1982), the mirrored Central Indian Basin and Crozet Basin between 
India and Antarctica (McKenzie and Sclater 1971; Sclater and Fisher 1974; Schlich 
1975, 1982) and the Wharton Basin between India and Australia (Sclater and Fisher 
1974; Liu et al 1983). In addition, during this period, Australia and Antarctica 
commenced to separate (Cande and Mutter 1982). The latest period in the evolution 
of the Indian Ocean started with a major reorganisation of the spreading centres, 
consequent to the collision of India with Asia in the Middle Eocene. Seafloor spreading 
ceased in the Wharton Basin while major changes in the direction and rate of spreading 
occurred in the Central Indian Basin, the Crozet Basin, the Madagascar Basin, the 
eastern Somali Basin and the Arabian Sea. The new configuration of the spreading 
ridges corresponds to the present-day ridge system. During this phase, the Australian- 
Antarctic Basin (Weissel and Hayes 1972), the southern Central Indian Basin (Sclater 
et al 1976a) and the northern Crozet Basin (Schlich 1975) were created. Also, the 
Mascarene Plateau and the Chagos-Laccadive Ridge separated (Fisher et al 1971; 
McKenzie and Sclater 1971), and the Gulf of Aden opened (Laughton et al 1970). 
Because of the differential motion along the Central Indian Ridge and the Southeast 
Indian Ridge, the Southwest Indian Ridge propagated rapidly towards the east 
(Tapscott et al 1980; Sclater et al 1981; Patriat 1987). 


3. Analysis of the altimetry data: Lineations in the deflection of the vertical charts 

The SEASAT satellite, launched by NASA in 1978, used a radar altimeter to make 
precise measurements (Tapley et al 1982) of the sea-surface or geoid height (~ 1-Om) 
over most of the world’s oceans. Unfortunately, the SEASAT mission failed after a 
few months, and almost no data were recovered from the areas south of 62°S because 
of the ice-coverage around Antarctica at the time that SEASAT operated. In March 
1985, the US Navy launched the GEOS AT altimeter in order to map the marine 
geoid to a high spatial resolution on a global basis. Only the data from the Exact 
Repeat Mission, started in November 1987, which duplicated the 17 days SEASAT 
repeat orbits, have yet been unclassified. However, because the ground track spacing 
(~ 164km at the Equator) decreases at the high latitudes, and because GEOS AT 
operated during the Austral ice-free summer, the information from the high southerly 
latitudes, between 60°S and 72°S, is exceptionally good. In this study, we use mainly 
the GEOSAT data. The accuracy is about 3 times better than for the SEASAT 
(Sandwell and McAdoo 1988) and we improve this accuracy and data recovery along 
tracks by stacking 22 GEOSAT repeat cycles. For wavelengths greater than 20 km, 
the nils noise level of the GEOSAT data is about 1 to 2 mGals. However, in order 
to take advantage of the better density of the SEASAT data north of 60°S, we combine 
both SEASAT and GEOSAT data sets when the amplitudes of the SEASAT signal 
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To enhance the short wavelengths of the marine geoid signal, it is more appropriate 
to use the deflection of the vertical which is4he first derivative of the geoid along the 
track. Deflection of the vertical, expressed in microradians, measures the variation 
of the horizontal gravity field. One microradian (^urad) is roughly equivalent to a 
horizontal variation in gravity of one milligal (mGal). Figures 2 and 3 display the 
deflection of the vertical plotted along respectively the descending and ascending 
satellite ground tracks. Vertical bars along tracks are spaced every 4 data points 
(~14km apart). Amplitudes are positive to the north. The scale of amplitude is 20 
^rad/degree of longitude; however when greater than 30 jurad, the amplitudes are 
reduced through an arctangent function. Long wavelength (>4000 km) and short 
wavelength noise (< 20 km) have been removed from each pass using the procedure 
de.scribed by Sand well and McAdoo (1988). 

Figure 4A through 4J display a series of different interpretations of the deflection 
of the vertical signal. In a tectonic fabric chart, we are mostly interested in the lineated 
features such as fracture zones, spreading ridges, linear troughs, ridges or plateaus, 
and trenches. Depending on their extension, most of these features have a reproducible 
signature on a set of parallel profiles, which make them easy to identify. The signature 
of a fracture zone will vary according to the spreading rate and the age offset (figures 
4A and 4B), Depending on the direction of the satellite and the orientation of the 
tectonic feature relative to the satellite ground track, the signature of a tectonic feature 
can be identical or reversed on the ascending or descending data sets. As a rule of 
thumb, topographic features oriented north-south will have same signatures; in the 
case of a fracture zone (figure 4C), both the ascending and descending profiles cross 
the topographic (age) offset from the same side. East-west oriented features will have 
reversed signatures; in figure 4D, ascending and descending profiles approach the 
fracture zone from opposite sides. Because the age offset reverses at the midpoint of 
the transform segment of a fracture zone, the signature of a fracture zone is symmetric 
relative to the fracture zone axis and of opposite sign on either side of this point 
(figure 4E). Finally, the spacing of parallel fracture zones is also an important factor; 
the combined effects of parallel topographic (age) offsets may produce a large signal 
where the signal produced by each individual fracture zone cannot be identified (e.g. 
Prince Edward FZ area along the Southwest Indian Ridge). The signature of a 
spreading ridge axis depends on the spreading rates as does the topography (figures 
4F and 4G). Seamounts have a typical signature (figure 4H) recognizable at isolated, 
sharp and narrow positive and negative picks of equal amplitude. Trenches (figure 
41) are characterized by a succession of positive and negative picks related to the 
flexure of the lithosphere ahead of the trench and to the trench itself. Continental 
rise (figure 4J) and linear offsets in the oceanic basement are also clearly identifiable. 
Finally, aseismic ridges and plateaus are clearly visible, but, depending on their width 
and height relative to the adjacent basins, they are difficult to chart. The deflection 
of the vertical charts of the Indian Ocean (figures 2 and 3) provides examples for 
each type of feature. We have traced lineations in these charts by correlating respective¬ 
ly peaks and troughs (i.e. steepest horizontal gradient of the gravity) in the deflection 
of the vertical plots, on both the ascending and descending passes. Figure 5 summarizes 
our interpretation of the lineations in the deflection of the vertical plots. We have 
abbreviated the name of this chart to deflection of the vertical lineations. Symbols 
show the actual location of the data points and lines, the lineations. We deliberately 
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Figure 3. GEOSAT ascending profiles of the deflection of the vertical plotted along the satellite ground tracks. Positive is the north. 
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Figure 4. Schematic signatures of the geoid and the deflection of the vertical (horizontal 
derivative of the gravity field) across some tectonic features of the ocean floor: (A), (B) Fracture 
zones. Since the deflection of the vertical is the horizontal derivative of the geoid along the 
satellite ground track, the sign of the signal changes with the direction of the satellite. 
(C), (D) Deflection of the vertical plotted at right angle of the subsatellite tracks. Positive is to 
the north (top of page). Descending and ascending profiles will have an identical signature 
across a north-south oriented fracture zone while they will show an opposite signature across 
an cast-west oriented fracture zone. (E) Variation of the deflection of the vertical signal 
along a fracture zone (same convention as for C and D). The signal reverses at the midpoint 
of the transform segment of the fracture zone as does the age offset. (F), (G) Deflection of the 
vertical signature above spreading ridge axes. (H) Seamount signatures generally show two 
sharp positive and negative picks of equal amplitude. (H) Trench (J) Continental rise. 


of the ocean floor, outlined in figure 5 by bathymetric contours. Such features are 
better represented by direct contouring of the gravity field (e.g. Haxby 1985). 


4. Tectonic fabric chart of the Indian Ocean 

The deflection of the vertical plots amplify the short wavelengths of the geoid. 
Lineations in these charts are in general related to high amplitude short wavelength 
topographic features such as transform faults, fracture zones, seamounts and other 
discrete entities on the seafloor. By combining the lineations of the deflection of the 
vertical charts with transform faults and fracture zones known from topographic and 
magnetic surveys, it is possible to create an improved tectonic chart of the sea floor. 
Because of the greater coverage of the satellite data and the ability of the presentation 
to amplify and extend lineated features both above and below sediment, the resultant 
signal presents a fabric chart of the igneous basement of the seafloor. We have 
abbreviated the name of these combined charts to “tectonic fabric charts”. 

The lineations in the deflection of the vertical charts (figure 5) are documented to 
varying degrees depending on the characteristics of seafloor spreading in the oceanic 
basins: spreading rate, paleo-direction of spreading relative to the satellite ground 
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tracks and the age offsets across the fracture zones. The interpretation of the deflection 
of the vertical charts brings new information mainly on the southeily latitudes around 
Antarctica regarding the early phase of seafloor spreading between Africa and 
Antarctica during the Late Jurassic up to the Late Cretaceous, and between Antarctica 
and Australia during the Middle Eocene. In addition, due to the uniform coverage 
of the satellite profiles, the fracture zones are generally better sampled and their 
pattern better defined on the deflection of the vertical charts than on tectonic charts 
based on shipborne data. This is for instance the case along the Central Indian Ridge, 
or between Broken Ridge and the Kerguelen Plateau, or between Tasmania and 
Antarctica. Except for this latter area, most of the fracture zones in the Indian Ocean 
are oblique or perpendicular to the ascending satellite tracks and parallel or 
sub-parallel to the descending tracks. Therefore, they are more visible on the ascending 
passes (figure 3) than on the descending passes (figure 2), while the opposite occurs 
for the Indian Ocean spreading ridges. However, in some areas such as the 
Australian-Antarctic Discordance where the fracture zones are oriented north-south 
and the ridge axes east-west, both data sets are complementary. 

In regard to the evolution of the Indian Ocean, the best documented period 
corresponds to the latest and most recent phase of seafloor spreading that started 
after the general reorganization of the plate boundaries in the Middle Eocene. Along 
the Sheba, Carlsberg and Central Indian Ridges, the fracture zones strike N50°E to 
N60°E. From the northwest to the southeast, the main fracture zones are the 
Alula-Fartak FZ (13°N,51°E) at the mouth of the Gulf of Aden, the Owen FZ (12°N, 
56°E) separating the Sheba Ridge from the Carlsberg Ridge, the fracture zone at 3°N, 
the Mahabiss FZ (3°S), the Sealark FZ (4°S), the Vitiaz FZ (8°S), the Verna FZ (1TS), 
the Argo FZ (14°S), the Marie Celeste FZ (I7°S) and the Egeria FZ (20°S) which lies 
east of the Rodrigues Ridge. Except for the Alula-Fartak FZ and the fracture zone 
at all the transform offsets are left-lateral. The largest transform offsets are 
located along the Alula-Fartak FZ (200 km), the Owen FZ (300 km), the Verna FZ 
(200 km) and the Marie Celeste FZ (220 km). Along the Southwest Indian Ridge, the 
orientation of the fracture zones varies progressively from N45°E in the vicinity of 
the Bouvet Triple Junction (55°S, 1°W) to north-south, west of the Indian Ocean 
triple junction (25°S,70°E). From west to east, they are Bouvet FZ (2°E), Moshesh 
FZ (5°E), Islas Orcadas FZ (6°E), Shaka FZ associated with the Shaka Ridge (8°E), 
Dingaan FZ (12°E), Mandela FZ (15°E), Du Toit FZ (26°E), the Andrew Bain FZ 
complex (28°E), Prince Edward FZ (35°E), Discovery II FZ (43°E), Indomed FZ 
(46°E), Gallieni FZ (52°E), the Atlantis II FZ (58°E) and the Melville FZ (61°E). The 
transform offsets are right-lateral and generally important (100 to 200 km) for all 
these fracture zones; the major offset (800 km) along the Southwest Indian Ridge axis 
occurs between Bain FZ and Prince Edward FZ (28°E to 35°E). Finally, along the 
Southeast Indian Ridge, the fracture zone trends, oriented N45°E in the vicinity of 
the Indian Ocean triple junction, shift progressively to N15°E-N10°E between 
Australia and Antarctica, and reach N15°W in the vicinity of the Macquarie triple 
junction (61°S, 162°E). The major offsets are located along the Amsterdam and St 
Paul FZ (78°E) and along the George V, Tasman and Balleny FZ. The fracture zone 
pattern between Australia and Antarctica, particularly the prominent fracture zones 
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and Antarctic margins. 

The spacing between satellite profiles is generally too large (160-80km) to allow 
an accurate mapping of the mid-oceanic ridge axes. However in several areas, the 
deflection of the vertical data shows an excellent correlation with the location of the 
ridge axes deduced from bathymetric and magnetic data. The Sheba Ridge is 
characterized by a slow spreading rate (Icm/a, Laughton etal 1970) and an axial 
valley that deepens towards the east. The ridge axis produces a clear signature on 
the deflection of the vertical signal, as illustrated in figure 4G, that defines two continuous 
ridge segments separated by the AIula-Fartak FZ. The Carlsberg Ridge has similar 
characteristics as the Sheba Ridge with spreading rate ranging from 1-2 to 1-3 cm/a 
(Le Pichon and Hcirtzler 1968; McKenzie and Sclater 1971) and shows a linear and 
nearly continuous ridge axis from the Owen FZ to the Equator. The Central Indian 
Ridge, which runs between the Mascarene Plateau and the Chagos-Laccadive Ridge, 
is extremely segmented by closely spaced fracture zones. The first ridge segments 
clearly identifiable out of the geoid signal caused by the fracture zones are located 
southeast of the Rodrigues Ridge. The spreading rates for this mid-oceanic ridge 
range from north to south between 1-8 and 2-4cm/a (Fisher etal 1971; Patriat 1987). 
The Southwest Indian Ridge has a very slow spreading rate (less than 1 cm/a; Schlich 
and Patriat 1971a; Sclater et al 1976b, 1981; Patriat 1987) and is expressed by a deep 
and pronounced inner valley. Several extensive ridge segments are well delineated by 
the deflection of the vertical data between the Mandela and Du Toit FZ which is the 
longest segment, east of the Prince Edward FZ, between the Discovery II, Indomed, 
and Gallieni FZ, and between the Melville FZ and the Indian Ocean triple junction. 
The Southeast Indian Ridge is the most active of the three Indian mid-oceanic ridges 
with medium spreading rates increasing from west to east from 3-0 to 3-5 cm/a (Schlich 
and Patriat 1971b; Weissel and Hayes 1972; Schlich 1975; Patriat 1987). Between 
the Indian Ocean triple junction and the Amsterdam and St Paul FZ, the ridge axis 
is segmented by several fracture zones and is generally well expressed by an axial 
valley, while east of the St Paul FZ the ridge segments are long and continuous and 
expressed by a central rise typical for medium or fast spreading ridges (Weissel 
and Hayes 1972; Royer and Schlich 1988). The area of the Australian-Antarctic 
Discordance presents an anomalous topography and fracture zone pattern (Vogt et al 
1983). The deflection of the vertical signal along the Southeast Indian Ridge varies 
from that of figures 4a and 4b. The trenches along the Java and Sumatra Islands and 
south of New Zealand are the second type of present-day plate boundaries that clearly 
show on the deflection of the vertical data. 

The altimetry data add little information on the fabric of the oceanic basins related 
to the second phase of opening of the Indian Ocean, from the mid-Cretaceous to the 
Middle Eocene. Only the Murray Ridge and the Chain Ridge, on either sides of the 
Owen Fracture Zone, document in the Arabian Sea and eastern Somali Basin the 
drift of India away from the Seychelles and East Africa. Structural directions in the 
Central Indian Basin and the Wharton Basin (McKenzie and Sclater 1971; Sclater 
and Fisher 1974), which record the fast northward drift of India relative to Antarctica 
and Australia, do not show well on the deflection of the vertical data. Their north-south 
orientation is too close from those of the satellite profiles. In the Wharton Basin, in 
addition to the Investigator FZ and the eastern scarp of the Ninety east Ridge, a cross 
grain is apparent, while intraplate deformation within the Central Indian Basin cause 
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east-west undulation of the lithosphere that obliterate and obscure the tectonic fabrics 
(e.g. Weisscl et al 1980; Geller et al 1983; Haxby 1985; McAdoo and Sandwell 1985). 
However, in the conjugate basins, the Madagascar and Crozet Basins, the fracture 
zone pattern is well recognizable. The most important feature is the Kerguelen 
Fracture Zone that offset the magnetic lincations 32. 33 and 34 (Schlich 198^- Patriat 
1987) by more than 800 km. With the GEOS AT data, it becomes clear that th’is feature 
does not extend further south than 62°S, and therefore post-dates the plate boundary 
reorganization in the mid-Cretaceous. In the Mascarene Basin, part of the extinct 
spreading system is visible. This system initiated at the time of the break-up between 
India and Madagascar in the mid-Cretaceous and died in the Paleocene/Early Eocene 
(chron 31 to 27; Schlich 1982) when the spreading ridge jumped to the north to form 
the Carlsberg Ridge and subsecjuently separated India from the Seychelles 3 nd the 
northern Mascarene Plateau. Remnant of the fossil spreading axis can be identified 
southwest of La Reunion Island (La Reunion Trench), north of the Mauritius Trench 
extension that forms the limit of the extinct spreading system with the Madagascar 
Basin. The fossil spreading ridge is offset 400 km to the north by the Mahanoro FZ 
and Wilshaw Ridge complex. Further north, the fracture zone pattern becomes unclear. 
Finally, the most interesting information related to mid-Cretaceous/Middle-Eocene 
time span concerns the relative motion of Africa and Antarctica. A clear set of parallel 
fracture zones can be defined between the Agulhas Plateau and the Maud Rise and 
Astrid Ridge, off Antarctica. At 48°S on the northern flank, and at 58°S, on the 
southern flank, there is a clear disruption of the fracture zone pattern which is shifted, 
by 50 to 100 km, to the east or to the west, respectively. This demonstrates the change 
of motion that occurred in the Paleocene (anomaly 31 to 24; Larson etal 1985; Patriat 
etal 1985; Royer etal 1988). This change of motion is particularly well recorded by 
the large offset northeast of Astrid Ridge (~40Ma age offset). Such a change of 
motion will also match perfectly the extremities of the fracture zones lying immediately 
west of the Madagascar Plateau (northern flank) and Conrad Rise (southern flank). 
The Prince Edward FZ is therefore a recent feature (Middle Eocene, chron 20) and 
is not the continuation of the fracture zone running from Astrid Ridge. 

Finally, new evidence about the early phase of opening of the Indian Ocean, from 
the Late Jurassic to the mid-Cretaceous, are found in the western part of the Enderby 
Basin. East of the fracture zone near Astrid Ridge, several parallel fracture zones are 
delineated in the deflection of the vertical charts. It shows that the direction of motion 


between Africa and Antarctica remained the same for 34 Ma, from chron MO (118 Ma) 
identified east of the Astrid Ridge (Bergh 1977, 1987) to chron 34 (84 Ma) located 
west of Conrad Rise (Bergh and Norton 1976; Patriat etal 1985). To the north, except 
for the Davie Ridge that records the relative motion between Africa and Madagascar, 
no fracture zone shows on the deflection of the vertical data: neither in the Mozambique 
Basin (Segouffin 1978; Simpson et al 1979), nor in the western Somali Basin (Segoufin 
and Patriat 1981; Rabinowitz et al 1983) or off the east African margin (Bunce and 
Molnar 1977; Cochran 1988). This is due to their north south orientation and to the 
very important thickness of sediments in these old basins. Along the western 
Australian margins, only the large offset Investigator FZ, the scarps of the Wallaby 
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produce a typical signature on the horizontal gravity field. Among the most spectacular 
are the conjugate escarpments of Broken Ridge and the Kerguelen Plateau. These 
lines correspond to an offset of the basement along the Kerguelen Plateau (Coffin 
etal 1986) and to a deep linear trough along Broken Ridge (Fisher etal 1982). 
According to Mammerickx and Sandwell (1986), these topographic features are related 
to the demise of seafloor spreading in the Wharton Basin (Liu etal 1983; Geller etal 
1983) and the initiation of spreading between the two plateaus (Mutter and Cande 
1983) at chron 18-19 (~45Ma). The match of these two curves provides a tight 
constraint upon the relative position of the Kerguelen Plateau and Broken Ridge 
prior to their break-up. The continent/ocean boundary mapped along the southern 
margin o^ Australia from magnetic and seismic evidence (Talwani et al 1979; Konig 
1980; Veevers 1986) has a distinct a'nd characteristic deflection of the vertical signature. 
We have used this criterion to identify the conjugate boundary along the Antarctic 
margin where the data are very limited. The deflection of the vertical signature also 
clearly outlines the limit of the continental shelf break (steepest gradient of the 
horizontal gravity field). Generally, this structural limit is well mapped from 
bathymetric data; however, around Antarctica where the data are sparse and scattered 
because of the remoteness and the ice-coverage, the altimery data now permit accurate 
delineation of the limit of the continental plateau. Another clearly visible feature is 
the change of topographic grain between the crust created along the Southwest Indian 
Ridge (at a slow spreading rate) and the crust created along the Southeast Indian 
Ridge and the Central Indian Ridge (at medium spreading rates). Actually, the age 
of this limit decreases towards the east, demonstrating the eastward propagation of 
the Southwest Indian Ridge (Patriat 1979; Tapscott et al 1980; Sclater et al 1981). 
Patriat (1987) mapped independently these two limits as the traces of the Indian 
Ocean triple junction respectively on the African and Antarctic plates. 


5. Conclusion: a revised tectonic summary diagram for the Indian Ocean 

Comparison of bathymetric data, deflection of the vertical lineations, and magnetic 
anomaly data permits the construction of a revised tectonic fabric chart for the Indian 
Ocean (figure 6). This chart is of significant value especially in the far southern ocean 
where new fracture zones have been identified and previously charted fracture zones 
extended. However, it is still only preliminary in nature as we have neither integrated 
the deflection of the vertical data with the actual topographic data to produce 
improved topographic charts (e.g. Driscoll etal., unpublished manuscript 1986) or 
repicked any magnetic anomalies. In addition, the geoid or gravity field of the oceans 
(Haxby 1985) contains very valuable information on the origin and structure of broad 
basic features of the ocean floor. The present study is the first step towards an 
international attempt to integrate all the above information, to provide some useful 
preliminary constraints for studies of the Indian Ocean tectonic history, and to show 
what can be accomplished by such integrated approach. 

For the sake of clarity, we did not represent in figure 6 all the magnetic lineations 
that have been identified in the Indian Ocean; in particular for the Cenozoic, we have 
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Abstract. Recently Brune and Singh (1986) reported evidence for anomalous continent-like 
crustal structure beneath the Bay of Bengal, possibly caused by perturbation in the 
temperature-pressure regime and consequent phase chage, partial melting, or mass transport 
(e.g. convection or underplating). Recent refraction results indicate the existence of an 
anomalous lower crustal or subcruslal layer of P-wave velocity about 7-3 km per second 
along the eastern North America passive margin, possibly a result of underplating of the 
oceanic crust just after initial rifting. We have searched for other evidence of anomalous 
crustal structure. The data suggest some mechanism may cause a general increase in the 
anomalous thickness of the crust with increasing thickness of the accumulated sediments, up 
to a thickness of about 6-7 km. On the other hand, anamolous crustal structure may in fact 
be transitional between oceanic and continental, or may have been modified by aseismic 
ridges, thus requiring no sediment related structure modification mechanism. The 
explanation for all the data may require more than one mechanism, all probably involving 
severe temperature perturbations. The general tendency is for perturbations of normal 
oceanic crust to make it more continent-like, suggesting that normal oceanic crust is an 
un.stable end-member in crustal states. 

Keywords. Anomalous crust; underplating; Bay of Bengal; sedimentary basin. 


1. Introduction 

Recently Brune and Singh (1986) reported evidence for anomalous continent-like 
crustal structure beneath the Bay of Bengal. Surface wave dispersion data could not 
be fit by a simple model of oceanic crust with overlying sediments—the model 
expected from current understanding of the tectonic evolution of the Bay of Bengal 
(Curray et al 1982; Royer et al this volume). The data required a progressively more 
continent-like crustal structure as the sediment load increased northward. This 
suggested that either the sediment deposition had altered the crustal structure, or that 
the structure had not been a simple oceanic type before deposition. If the anomalous 
structure were a result of modification by the sedimentary load it might be caused by 
the perturbation in the temperature-pressure regime and consequent phase change, 
partial melting or mass transport (e.g. convection or underplating), and this might 
have important implications not only for the Bay of Bengal, but also for the classical 
geodynamic problem of basin formation. 

Recent results of Trehu et al (1988) have increased interest in the problem. Using 
ocean bottom seismometers they were able to show the existence of an anomalous 
lower crustal or subcrustal layer of P-wave velocity about 7-3 km per second. They 
suggested that this could be the result of underplating of the oceanic crust just after 
the initial rifting of the North American plate from the African Plate. These results are 
discussed in more detail below. 
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what might have been expected to be normal oceanic crust is so similar to the i 
for the Bay of Bengal that it suggests a common causal factor. The fact that for tl 
of Bengal the perturbation in structure correlates with sediment thickness (incr 
from south to north) suggests that some sediment load effect should be consider 
the Carolina Trough results. On the other hand, whatever mechanism is estab 
for the Carolina Trough results should be considered as a possible explanation 1 
Bay of Bengal results. 

If the original structure beneath the Bay of Bengal was not oceanic, then 
modifications of the current plate tectonic model for the Indian Ocean mi^ 
suggested. However the excellent fit between the coast of India and the co 
Antarctica (Lawver and Scotese 1987) argues against any other structure thai 
oceanic in the majority of the Bay of Bengal. 

An alternate explanation of the anomalous dispersion might be a combinat 
structural perturbations, including effects of a marginal rifting event 
perturbations by the 90 deg east ridge and the 85 deg east ridge (e.g. th 
underplating). 


2. Evidence from other areas 

Brune and Singh (1986) mentioned evidence that anomalous structure unde 
thick sediments on oceanic structure might not be unique to the Bay of Beng 
some anomalous surface wave dispersion evident in figure 8 of Sykes and Oli 
Part II (1964) for the Argentine sedimentary basin. We have searched for 
evidence of the phenomenon, and report the results here. Most of the evidence ii 
published studies. A paper documenting anomalous crustal structure in small 
basins was published by Menard (1967). Although his results may be related 
data considered here, we have not placed heavy emphasis on it because in mosi 
it is not clear that the original structure before sedimentation was oceanic. 

We have not included results from active tectonic areas such as the Bering S( 
Mediterranean, and the Sea of Okotsk, which, although they have consid 
thicknesses of sediments, have not been stable tectonically in the same man 
passive margins and the Bay of Bengal, and in some cases seem to have n 
oceanic crust underlying the sediments, possibly a result of creation of new cr 
back-arc spreading or some other type of active tectonism. 


3. Data 

Selected data correlating the sediment thickness to the thickness of the anon 
crust is plotted in figure 1, and is discussed below. 

3.1 The Bengal Fan 

The data for the Bengal Fan are taken from Brune and Singh (1986). For a disc 
of that data the reader is referred to the publication. The main conclusion frc 
data is that there is an apparent correlation between sediment thickness and thi( 



Thickness of Anomalous (excess) Crust (km) 

Figure 1. Plot of selected data for thickness of sediments versus thickness of anomalous 
crust, i.e. crustal thickness in excess of 6 km. 


of anomalous lower crust (figure 1). This apparent correlation suggests that some 
causal mechanism increases in effect from south to north, rather than parallel to the 
coast as might be expected from a remanent continental margin effect or the initial 
rifting event or from effects of the 85 East and 90 East ridges. However since the width 
of the Bay of Bengal decreases from south to north it would take less of a 
contaminating structure to influence the apparent average crustal thickness. 

3.2 The Arabian Sea 

Singh (1988) has recently published results on surface wave dispersion across the 
Arabian Sea. The Indus Fan sedimentary basin in not as large as the Bengal Fan and 
the location of seismic stations and earthquake sources is not as favourable, with no 
good paths crossing the thick part of the fan. Many of the paths are contaminated by 
the Maidive and Laccadive ridges. The dispersion results suggest the existence of 
anomalous crust for much of the Arabian Sea area studied by Singh (1988). 

3.3 Anomalous crustal structure at passive continental margins 

Recently Trehu et al (submitted 1989) reported new results for the crustai structure 
underlying the sedimentary basin of the Carolina trough. These results were obtained 
using new ocean bottom seismographs which recorded clear deep crust and mantle 
arrivals, arrivals not reliably recorded in earlier studies using surface sonobouys. The 
data provided evidence for a thick (up to 15 km) high-velocity (7-2-7-6 km/s) lower 
crustal or upper mantle layer extending beneath the Carolina trough and the adjacent 
ocean basin. This layer was either very thin or absent landward of the hinge zone 
(continental broader). Trehu et al, following similar results from the Baltimore Canyon 
Trough (EASE Study Group 1986), interpreted this anomalous crust to represent 
material that was underplated to the crust during the late stage of continental drifting 
and early seafloor spreading. 

Both the Carolina Trough and The Baltimore Canyon have very thick sedimentary 
accumulations, similar to the thicknesses of sediments in the Bengal Fan, and thus fit 
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the correlation of anomalous crust with sediment thickness. However results fro 
Hatton Bank, a sediment starved margin, also indicate the presence of an axiom 
layer which Trehu et al (1989) label as “underplated”, but it is more localized 
the continental border. 


3.4 Gravity 

A detailed study of gravity using data collected in recent cruises has not been c: 
out in the Bay of Bengal. Liu et al (1982) presented a free-air gravity map of the B 
Fan. The free-air anomalies tend to be 20 to 30 milligals negative over the w< 
part, and near zero close to the eastern margin near the Ninety-East Ridge. If m; 
sediments were flooded onto a perfectly rigid oceanic crust-mantle system the fr 
anomaly would be positive since more dense sediments would replace less 
water. However in a non-rigid system isostatic compensation will cause the anc 
to move more or less rapidly toward zero. Liu etal (1982) showed that the eff 
flooding sediments around a ridge, in this case the buried Eighty-Five East ridge 
create a local gravity low over the ridge, but we know of no mechanism which \ 
give a broad low over a region where turbidity currents are rapidly depo 
sediments. The negative values are probably related to more global factors, e. 
assumed shape of the geoid. We do not believe that at this stage gravity provi 
reliable constraint on the crust-mantle structure, or on the interpretation c 
surface wave data. 

Rao and Rao (1986) made an approximate gravity interpretation assuming 
sedimentary layers of density 2*3 and 2*5 g/cc, respectively, and a basement dens 
2*84 g/cc overlying a mantle with density 3*4 g/cc. Their interpretation indica 
slight crustal thinning northward in the Bay of Bengal; however the interpret 
depends on unverified density contrasts, does not consider a possible lower c: 
layer, and does not deal with the question of local and regional compensatioi 
thus cannot resolve the questions being dealt with in this paper. 


4. Discussion 

Figure 1 is a plot of sediment thickness versus thickness of the anomalous crustal 
(crustal thickness remaining after subtracting out a “normal” crustal thickn( 
6 km). The Bengal Fan data are connected by a line, starting at the origin (( 
means no sediments and no excess crustal thickness). The data are consistent 
some mechanism causing a general increase in the thickness of the crust 
increasing thickness of the accumulated sediments, but there is considerable scat 
the data and the correlation does not establish any particular mechanism to ej< 
the data. If there is a causal mechanism relating anomalous crustal thickne 
sediment thickness, figure I suggests that it is not linear but has a “corner” 
sediment thickness of about 6-7 km, possibly suggesting a saturation effect sc 
beyond this thickness, further sediments have less of an effect in modifyin 
structure. 

Possible general interpretations of the data include: 


by a mechanism yet to be established e.g. thermal blanketing, phase change, 
underplating; etc. (causal sedimentary effect). 

(ii) Since all the data are from oceanic areas near continental margins, or areas with 
aseismic ridges, the actual underlying crust may not have been (before sedimentation) 
originally normal oceanic as assumed, but in fact transitional between oceanic and 
continental, or may have been modified by aseismic ridges, and thus no sediment 
related structure modification mechanism is required. 

There is no established mechanism for a causal connection between sediment load 
and crustal thickness. The data suggest that the modification in structure occurs near 
the Moho, but there is strong evidence from our recently developed understanding of 
the evolution of oceanic crust that the oceanic Moho is a compositional boundary, 
not a phase change or pressure-temperature boundary. Thus it is difficult to see how 
sediment load could lead to large structural changes near the Moho. Furthermore the 
most obvious candidate for a pressure effect, the basalt-eclogite phase change would 
work the opposite direction and cause the Moho to shallow. This suggests that if there 
is a large structure modification near the Moho, it must involve either temperature or 
mass transport, or both. One such mechanism might be reactivated basalt 
underplating brought about by an increased geothermal gradient, possibly from the 
sedimentary blanketing effect and/or reactivated shallow convection in the mantle. 
Neither of these mechanisms has been quantitatively evaluated. 

Since great sedimentary basins are of necessity near continents (to provide a supply 
of sediment), there is usually some question as to where the exact continent-ocean 
boundary is, and whether or not the structure of the crust underneath the basin is 
truly oceanic. An exact knowledge of plate tectonic history may be required to 
establish the original type of structure unambiguously. In the case of the Bengal Fan 
the location of the remnant boundary between the Antarctic and Indian plates is well 
established from the deflection of the vertical data (Royer et al this volume), and the 
most recent plate reconstructions (Lawver and Scotese 1987), and thus we expect the 
Bengal Fan to be underlain by oceanic crust. 

Even if it is accepted that the Bay of Bengal is underlain by oceanic crust, there is 
the possibility of structural perturbations not only by continental margin effects but 
also by anomalous structures such as aseismic oceanic ridges. The nature of the 
Ninety-East Ridge is not known for certain, and conceivably could exert some 
influence on the surrounding structure. The same is true for the buried Eighty-Five 
Degree East Ridge. 


5. Conclusion 

Evidence for anomalously thick crust in areas expected from plate tectonic 
reconstructions to have normal oceanic crust suggest that various factors may be 
acting to perturb the crust-mantle structure. Evidence from the Bay of Bengal suggests 
a strong correlation with sediment thickness. Evidence from eastern U.S. passive 
continental margin sedimentary basins suggests a similar anomalous crustal structure 
related to the original rifting event. If there is a strong effect related to sediment 
blanketing it could result from temperature perturbations and related mass transport 
phenomena, e.g. underplating. The explanation for all the data may require more than 
one mechanism, e.g. mechanical mechanisms related to stretching and rifting. 



jiicuiiauisms rciiiicu lo iriciiimi cucuis oi scuiiiicxiLciij uiauN.cuiig, auu znet 
related to emplacement of aseismic ridges. All of these probably involv 
temperature perturbations. Whatever mechanisms are operating, the general t 
for perturbations of normal oceanic crust is to alter the structure in such a v 
make it more continent-like. Perhaps this indicates that normal oceanic cn 
unstable initial crustal state, and various mechanisms can restructure it toward 
continental structure. 
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High velocity anomaly beneath the Deccan volcanic province: 
Evidence from seismic tomography 
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Abstract. Analysis of teleseismic P-wave residuals observed at 15 seismograph stations 
operated in the Deccan volcanic province (DVP) in west central India points to the existence 
of a large, deep anomalous region in the upper mantle where the velocity is a few per cent 
higher than in the surrounding region. The seismic stations were operated in three 
deployments together with a reference station on precambrian granite at Hyderabad and 
another common station at Poona. The first group of stations lay along a west-northwesterly 
profile from Hyderabad through Poona to Bliatsa. The second group roughly formed an 
L-shaped profile from Poona to Hyderabad through Dharwar and Hospet. The third group 
of stations lay along a northwesterly profile from Hyderabad to Dhule through Aurangabad 
and Latur. Relative residuals computed with respect to Hyderabad at all the stations showed 
two basic features: a large almost linear variation from approximately -I- 1 s for telcseisms 
from the north to — 1 s for those from the southeast at the western stations, and persistance 
of the pattern with diminishing magnitudes towards thp east. Preliminary ray-plotting and 
three-dimensional inversion of the P-wave residual data delineate the pre.sence of a 600 km 
long approximately N-S trending anomalous region of high velocity (1-4% contrast) from 
a depth of about 100 km in the upper mantle encompassing almost the whole width of the 
DVP. Inversion of P-wave relative residuals reveal the existence of two prominent features 
beneath the DVP. The first is a thick high velocity zone (1-4% faster) extending from a 
depth of about 100 km directly beneath most of the DVP. The second feature is a prominent 
low velocity region which coincides with the westernmost part of the DVP. A possible 
explanation for the observed coherent high velocity anomaly is that it forms the root of the 
lithosphere which coherently translates with the continents during plate motions, an 
architecture characteristic of precambrian shields. The low velocity zone appears to be 
related to the rift systems (anomaly 28,65 Ma) which provided the channel for the outpouring 
of Deccan basalts at the close of the Cretaceous period. 

Keywords. Deccan volcanic province; seismic tomography; deep structure; high velocity 
anomaly; teleseismic residuals; three-dimensional inversion. 


1. Introduction 

The Deccan volcanic province (DVP) encompasses an area of about a half million 
km^ in west-central India with a possible extension of an additional 1 million km^ 
beneath the Arabian Sea to the west. The Deccan basalts are primarily of tholeiitic 
composition with a few pockets of acidic volcanism. The volcanic flows occur as a 
relatively thin veneer at the surface with a thickness of 1-1-5 km in the western third 
of the province, thinning to about 0-4 km in the middle section and to a few meters 
in the eastern section (Kaila 1982). The bulk of the Deccan volcanism is now known 
to have occurred around• 65 Ma (Mahoney 1988). About 200 Ma ago, India was 
wedged in the southern latitudes between Africa and Antarctica as part of the 


migrating northwards as an independent continent eventually colliding with Ei 
about 45 Ma ago. Morgan (1972) proposed that the Deccan basalts erupted c 
this northward migration as the sub-continent traversed over a mantle | 
associated with the Reunion hot spot. 


1.1 Gravity 

The Deccan Traps lie in a broad gravity low which encompasses the whole 
south Indian shield. Within this broad low are shorter wavelength “lows” and “h 
A few of these short wavelength features have been interpreted variously: upli 
subsidence due to cycles of volcanic activity (Kailasam 1975); crustal intrusi 
mafic rock (Takin 1966) and linear channels of magma flow (Qureshy 1981). W 
show later that preliminary two-dimensional modelling of the gravity data ind 
a lower-than-normal density in the upper mantle beneath the central part of the 

1.2 Deep seismic sounding 

Deep seismic sounding was carried out along two profiles over the Koyna reg 
the Deccan Traps and several others in the Saurashtra Peninsula and acro: 
Narmada-Son lineament (Kaila et al 1981, 1985; Kaila 1982). Though the ve 
models obtained from these profiles reflect the thickness of the surficial basa 
to 2 km, the crustal structure beneath the DVP is not significantly different fror 
of the precambrian south Indian shield model elsewhere as revealed by data 
the Kavali-Udipi profile (Kaila 1982). Moho depth is found to vary from 39-^ 
west of Koyna to 36-38 km eastwards. There is no evidence of drastic changes 
depth of major discontinuities like the Conrad discontinuity across the profi' 
is there evidence of unusually high upper or lower crustal velocities whic 
characteristic of rift zones. 

1.3 Heat-flow 

Except for the Cambay graben where the heat-flow is somewhat higher-than-nc 
other parts of the Deccan Traps show normal heat flow typical of stable contii 
regions (Gupta and Gaur 1984). Except for implying that the volcanic ej 
responsible for the Deccan volcanic flows did not leave a detectable thermal rema 
in the lithosphere, the heat-flow data do not provide any information on the 
structure beneath the DVP. 

In summary, the available geophysical data are inadequate to constrain a s 
model for the Deccan volcanism. For example, if the magmatism was causec 
mantle plume, extensive melting, differentiation and variations in physica 
chemical properties of the deep structure of the lithosphere might have taken 
But the available heat-flow and deep seismic sounding data are inadequate to i 
the structure of the lithosphue and asthenosphere beneath the DVP which ’ 
expose alterations in structure associated with that volcanism. With this 
ground, we designed a tomographic experiment to model in three dimensioi 
compressional wave velocity structure of the upper mantle beneath the DVP. R 
from this experiment reveal a vast anomalous region in the upper mantle 


DVP where the P-wave velocity is a few per cent higher than in the surrounding 
region. This is the first effort to map the deep structure beneath the DVP. 


2. Tomography using teleseismic residuals 

Teleseismic residuals provide an excellent approach to model the two- and three- 
dimensional velocity structures of the earth’s crust and mantle. Numerous publications 
describe the use of the technique to delineate heterogeneous earth structure from 
global and continental scales to local scale. Recent introduction of inversion methods 
for teleseismic residuals (also called seismic tomography) has considerably enhanced 
the power of the teleseismic residual technique (Aki et al 1977; Humphreys and 
Clayton 1988). Tomographic modelling of small and large heterogeneities in the 
earth’s crust and mantle have been carried out in diverse geologic settings (Thurber 
and Aki 1987; Iyer and Hitchcock 1988). The present study is the first application of 
this technique in India. 


3. Estimation of teleseismic residuals 

Teleseisms (earthquakes occurring at distances greater than approximately 20 degrees 
from the recording station) are recorded by a dense array of seismic stations sited 
over the region of interest to provide the data for tomography. The average instrument 
spacing in the array limits the smallest size of the velocity anomaly that can be 
resolved. Two-dimensional arrays yield three-dimensional models, while linear arrays 
yield two-dimensional models. Arrival times of the first seismic phases (P or PKIKP, 
PKP) are read from the seismograms and travel time from the source to the recording 
stations are estimated from a knowledge of the earthquake origin times. Knowing 
the hypocentral location of the earthquake, theoretical travel times are computed 
using a standard seismological table. In the present study, we used the P-wave travel 
time tables of Herrin (1968a). The travel-time residuals are simply the differences 
between the observed and theoretical travel times. In the mathematical notation we 
write, 

^ ij ^ ij 

where, is the absolute residual in seconds, and 7^ and TjJ are the observed and 
theoretical travel times respectively at station i for event j. In addition to the 
information on the deep structure beneath the seismic array, Rij contains effects due 
to errors associated with the location and origin times of the earthquakes travel-time 
tables, and anomalous structure at the source and along the seismic wave paths. In 
order to isolate the local structures from extraneous effects, relative residuals are 
computed independently for each event by subtracting the residual at a reference 
station which hopefully is located outside the anomaly under investigation from the 
absolute residuals at each station. Alternatively, a weighted mean of absolute residuals 
for each event can be used as the reference residual. Thus, the relative residual RRij 
for event j can be defined as, 

RRij = Rij — R^j, 



where the subscript r represents the reference station, or as 


RRij = Rij'Rr 

where R, represents the weighted mean of Rij using all readings available for the 
event j. In this study all weights were chosen to be equal to unity. Consequently, Rj 
is a simple average of absolute residuals for the whole array. 

RRij contains the information needed to model the three-dimensional relative- 
velocity structure beneath the seismic array. The magnitudes of the relative residuals 
indicate the change in velocity within the anomaly; the sign of the residuals being 
negative for high-velocity anomalies and vice versa. The azimuthal pattern of RRij 
gives an idea of the depth of the anomaly. As a rule of thumb, it may be remembered 
that a 10-km long anomalous ray-path in which the seismic velocity differs from the 
normal by 10% produces a residual of about 0-2 s. 


4, Three-dimensional inversion 

Wc have used the three-dimensional inversion technique developed by Aki et al (1977) 
{now popularly referred to as the ACH technique) to model the teleseismic residuals. 
Since the ACH technique is described in detail in various papers (Aki etal 1977; 
Ellsworth and Koyanagi 1977; Iyer et al 1981; Achauer et al 1986) and also its 
application (Aki 1982; Thurber and Aki 1987; Iyer and Hitchcock 1988), only a brief 
outline of the theory will be given here along the lines presented by Acha"er et al 
(1986). In this technique, the volume under the array is divided into layers and each 
layer into a grid of rectangular blocks. An initial T-wave velocity from a known or 
assumed model is assigned to each layer for ray-tracing purposes. It turns out that 
the final results are not very sensitive to the choice of this initial velocity model. The 
incident wave at the bottom of the deepest layer is assumed to be plane. As it travels 
to the seismograph, it encounters an unknown velocity perturbation in each block. 
The resultant travel time delays and advances accumulated along the ray path equal 
the observed residual at the surface. In this treatment, the refraction of the seismic 
ray due to the velocity perturbation is assumed to be negligible. Thus for each ray, 
the following linear equation can be generated in a matrix form. 

Am = d, 

where ^4 is a matrix with calculated unperturbed travel times of ray segments, m is 
a vector containing the unknown fractional slowness perturbations (slowness is the 
reciprocal of velocity), and d is a vector containing the travel time residuals. Since 
A'^ A is near singular, Aki cr a/(1977) used the damped least-squares method to obtain 
the model estimate m as, 

m = {A^A + Q^iy '■A^d, 

where 6^ — is the damping parameter and cr| and are the estimated variances 
of the data and the (unknown) true model m (Achauer etal 1986). It should be 
remembered that the estimated model contains only relative velocities, as the inputs 
are relative residuals. 

:_•__ i _i _ / t. • r •. 


Lt it pfovides only a ‘fuzzy’ picture of the real model. The degree of‘fuzziness’ 
^ tj-ie density of instruments deployed in the seismic network, the number 
thal distribution of teleseisms used and the block size. Fortunately, the 
^g|._g<q[uares inversion also produces a resolution matrix and standard errors 
t>le to estimate the figure of merit for a computed model. The ‘true’ 
5 related to the ‘fuzzy’ model m by the relation 


i tlie resolution matrix. Each element of the symmetric singular matrix R 
le t>etween — 1 and +1 with positive values on the diagonal. Therefore, 
Df the matrix can be considered as an average kernel relating all of m, the 
I, to one block of the estimated ‘fuzzy’ model. Thus large diagonal elements 
olu-tion matrix imply that m is a good representation of m. In general, 
is good towards the middle of each layer and is rather poor near the edges, 
ly related to the number of ray-hits in the block. The standard errors, 
:>m the diagonal of the covariance matrix, are indicators of the accuracy 
elocity perturbation. As a rule of thumb, the accuracy of individual 
on estimates are about (±) 2 times the standard error for 0-8 resolution 
—4 times for 0-6 resolution (Evans 1982). When small blocks are used, the 
the inversion can be assessed by looking at the whole resolution matrix 
blocks. 


periment 

ted fifteen vertical-component seismic stations (Teledyne protacorders with 
s, seismometers) in and around the southern half of the DVP (figure 1, 
1 three consecutive deployments with Hyderabad (HYB) and Poona (POO) 
imon to the whole experiment. A total of 168 teleseisms mainly from north, 
east and southeast, and in the distance range of 19 to 95 degrees were 
by this network (table 2). The recorders were operated at their maximum 
ed of 240mm/min and internal clocks synchronized with radio time 
oadcast by the National Physical Laboratory, New Delhi. Considerable 
taken in reading the arrival times of the first phases, thereby achieving a 
racy of 0 05 s. Visual correlation of seismic waveforms and multiple readings 
k, peak, zero-crossing, trough etc) were used to ensure a high relative 
of picks (figure 2). Despite these precautions substaptial signal changes 
; array led us to discard about 10% of the readings (figure 3). This as well 
'' amplitudes of the signals also contributed to the scatter in data. Even so, 
magnitude of the observed relative-residuals permits us to draw significant 
as. 


nalysis 

tion of residuals with azimuth and distance 

"7 show the relative residuals, with resnp.r.t to HYTt at 1"^ statinn« nlntt#»r1 



Figure 1. Map of west-central India showing the locations of seismic stations (three-letter 
names) used in this study. The approximate outline demarking the extension of the DVP. 

Table I. Coordinates of seismic stations. 


Station name 

Latitude 

(N) 

Longitude 

(E) 

Elevation 

(meters) 

Duration of 
operation 

Hyderabad (HYB) 

17“25-03' 

78°33-18' 

5100 

Sept. 1985-Dec. 1987 

Poona (POO) 

18°32-84' 

73“49-07' 

579-3 

Sept. 1985-Dec. 1987 

Bhatsa (BHA) 

19°3M2' 

73°24'06' 

185-9 

Dec. 1985-Apr. 1986 

Ujani (UJN) 

18“04-00' 

75°06'87' 

496-9 

Jan. 1986-Apr. 1986 

Solapur (SOP) 

17“38-90' 

75°54'93' 

475-6 

Jan. 1986-May. 1986 

Gulbarga (GUL) 

17019.57 

76“53-25' 

470-0 

Jan. 1986-Apr. 1986 

Satara (SAT) 

17“40-17' 

74°00-19' 

670-5 

Jun. 1986-Feb. 1987 

Kolhapur (KOP) 

I6°40-20' 

74° 15-10' 

563-8 

Jun. 1986-Feb. 1987 

Dharwar (DWR) 

15‘’25-27 

74°55-79' 

731-5 

Jun. 1986-Jan. 1987 

Hospet (HPT) 

15°15-31' 

76°20-58' 

478-5 

Jun. 1986-Jan. 1987 

Srisailam (SLM) 

16°06 00 ' 

78°53-50' 

438-9 

Jan. 1987-Mar. 1987 

Aurangabad (AUB) 

I9“54-38' 

75°18-66' 

569-0 

May. 1987-Dec. 1987 

Bir(BIR) 

I8°58-37' 

75°46-28' 

560-0 

May. 1987-Dec. 1987 

Latur (LTR) 

18°21-25' 

76°33-75' 

610-0 

May. 1987-Dec. 1987 

Dhule (DHL) 

20°55-49' 

74°47-0r 

265-0 

Jun. 1987-Dec. 1987 


Table 2. List of teleseisms used in this study which provided readings in at least 3 stations. 


SI. 

No. 

Source region 

Origin time 

Depth 

(km) 

Coordinates 

1 . 

MOLUCAzc(01/03/86) 

09-43-28-6 

39 

00:58-32S 

126;52-08E 

2 . 

AFGHAN(01/14/86) 

03-03-37-4 

245 

36:20-46N 

71:01-44E 

3. 

HONSHUzc(01/16/86) 

08-34-43-7 

436 

29:48-12N 

138:39-72E 


(continued) 
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Table 2. (contiinied) 


SI. 

No. 

Source region 

Origin time 

Depth 

(km) 

Coordinates 

4. 

MINDROzc(01/16/86) 

15-45-06-7 

223 

13.-41-10N 

120.-49-32E 

5. 

ALUTINtr(01/18/86) 

01-5901-6 

33 

51:33-18N 

173:06-54W 

6 . 

MOLUCAtr(01/22/86) 

14-57-13-0 

59 

00:28-028 

124:21-96E 

7. 

HOKADOtr(01/31/86) 

17-48-04-4 

69 

42:13-68N 

142:59-76E 

8 . 

CELBEStr(02/01/86) 

15-00-35-8 

392 

02.-54-36N 

124.-04-50E 

9. 

VANATUzc(02/02/86) 

01-44-05-4 

31 

13:37-568 

166.-41-58E 

10 . 

BONINtr(02/03/86) 

20-47-35-3 

508 

27:47-46N 

139:33-12E 

11 . 

HINDKUzc(02/04/86) 

1943-56-6 

203 

36:25-56N 

70:41-76E 

12 . 

AFGHANpk(02/ll/86) 

07-14-24-2 

119 

36.-22-32N 

70:54-60E 

13. 

HONSHUpk(02/12/86) 

02-59-30-4 

30 

36:23-04N 

14l:07-68E 

14. 

PAPNGUzc(02/12/86) 

11-27-45-4 

34 

06:32-168 

147:25-98E 

15. 

HALMHApk(02/12/86) 

2302-12-8 

135 

01:36-24N 

127:18-30E 

16. 

CARSBR2c{02/15/86) 

19-56-35-9 

10 

04:24-72N 

62:42-24E 

17. 

KURILEpk(02/19/86) 

10-54-46-2 

115 

48:34-74N 

153:24-90E 

18. 

PHLPNEzc(02/19/86) 

11-40-27-5 

77 

18:56-70N 

121:18-24E 

19. 

AUSTLAzc(02/24/86) 

01-26-58-2 

32 

12:45-248 

114.-36-54E 

20 . 

HALMHApk(02/24/86) 

02-31 •26-7 

118 

01:43-74N 

127:21-24E 

21 . 

BANDCApk(03/01/86) 

16-41-40-6 

80 

06:18008 

130:56-94E 

22 . 

KAMCHK(03/02/86) 

03-14-41-8 

118 

51:40-68N 

156:56-16E 

23. 

SOLMONzc(03/06/86) 

12-31-24-0 

72 

07:00-008 

155:46-86E 

24. 

NEWBRTfb(03/07/86) 

02-46-52-0 

116 

04:59-408 

151:42-60E 

25. 

HANDUKzc(03/ll/86) 

23-07-38-2 

206 

36:28-92N 

70:39-60E 

26. 

BONINzc(03/17/86) 

09-18-25-2 

476 

27:25-26N 

139:51-96E 

27. 

FLOREStr(03/21/86) 

21-35-35-8 

614 

07:27-308 

120:38-76E 

28. 

HONSHUtr(03/25/86) 

04-00-43-3 

144 

37:15-72N 

139:25-38E 

29. 

B'RAZILzc(03/26/86) 

22-06-57-6 

609 

07:07-508 

71:38-28W 

30. 

NBRTANzc(04/03/86) 

07-39-53-9 

30 

06:20-408 

151:42-90E 

31. 

KURILEzc(04/05/86) 

22-59-09-4 

80 

44-27-96N 

147:51-42E 

32. 

MINDNOtr(04/09/86) 

21-56-19-6 

56 

09:5 8-14N 

126:09-24E 

33. 

MOLUCAzc(04/10/86) 

02-21-13-3 

33 

00:58-148 

126:50-28E 

34. 

MARANAzc(04/13/86) 

03-00-20-9 

304 

17:11-28N 

145:37-26E 

35. 

VANATUzc(04/14/86) 

00-25-12-4 

33 

13:56-768 

166:58-62E 

36. 

KURILEzc(04/16/86) 

12-52-19-7 

52 

43:52-98N 

147:33-96E 

37. 

TANMBRzc(04/18/86) 

08-08-39-2 

44 

06:00-428 

131:32-82E 

38. 

WIRIANzc(04/20/86) 

07-03-30-8 

33 

02:22-508 

139:19-08E 

39. 

CRETE tr(04/27/86) 

09-27-02-3 

10 

34:38-46N 

.23:24-42E 

40. 

SUMTRAzc(04/29/86) 

13-59-21-8 

47 

04:2400N 

95:00-24E 

41. 

RYUKYUtr(04/30/86) 

23-14-43-7 

44 

28:42-66N 

130:02-82E 

42. 

PAPUAGzc(06/24/86) 

03-11-30-9 

102 

04:26-888 

143:56-58E 

43. 

FLOREStr(06/27/86) 

03-09-50-8 

229 

07:56-828 

122:47-46E 

44. 

FLORESpk(06/27/86) 

03-09-50-8 

229 

07:56-828 

122:47-46E 

45. 

TAGYKApk(06/29/86) 

21-47-59-6 

20 

05:20-168 

20:32-34E 

46. 

MEXICOpk(07/05/86) 

22-09-36-9 

112 

15:28-02N 

92:34-92W 

47. 

CARBRGfb(07/07/86) 

16-26-56-6 

8 

10:23-34N 

56:49-92E 

48. 

MINAHAzc(07/08/86) 

04-27-34-6 

245 

01:59-82N 

124:18-30E 

49. 

MOLUCApk(07/09/86) 

23-10-53-1 

28 

01:54-24N 

126:31-50E 

50. 

SIRANfb{07/12/86) 

07-54-26-8 

10 

29:57-72N 

51:34-92E 



Table 2. (conliimecl) 


SI. 

No. 

Source region 

Origin time 

Depth 

(km) 

Coordinates 

51. 

AFUSSRzc(07/17/86) 

15-46-37-0 

47 

36:40-08N 

71:14-82E 

52. 

FOXIS zc(07/19/86) 

04-31-55-9 

33 

53:21-12N 

165:52-92W 

53. 

KURILEzc(07/19/86) 

05-59-36-2 

141 

47:15-84N 

151:07-62E 

54. 

BURMAtr(07/26/86) 

20-24-47-8 

27 

23:45-18N 

94:10-62E 

55. 

EUSSRpk(07/26/86) 

14-46-18-9 

310 

45;25-32N 

137:04-20E 

56. 

CHINApk(08/06/86) 

19-55-15-6 

34 

29:20-64N 

100:54-90W 

57. 

PHILPNzc(08/09/86) 

00-53-12-6 

80 

14;05-94N 

120:20-1OE 

58. 

HALMRAzc{08/10/86) 

04-40-49-7 

104 

01:59-10N 

128:I6-26E 

59. 

MOLUCAtr(08/17/86) 

15-27-41-2 

31 

02:16-56N 

126:57-54E 

60. 

BANDASzc(08y 19/86) 

12’5241-1 

33 

04:07:32S 

129:19-44E 

61. 

PHILPNzc{08/19/86) 

22-41 •35-8 

98 

12:30-84N 

124:28-68E 

62. 

SIOCENtr(08/20/86) 

21-15-46-4 

10 

01-.58-38S 

87:05-28E 

63. 

ROMANAzc(08/30/86) 

21-28-35-4 

132 

45:32.82N 

26:18-86E 

64. 

VOLCNOpk(08/31/86) 

09-23-04-7 

46 

23:02-52N 

144:07-20E 

65. 

KURILEtr(09/01/86) 

21-47-35-0 

148 

46;37-86N 

150:00-60E 

66 . 

HALMRAzc(09/09/86) 

15-58-22-5 

33 

00:55-62N 

127:56-64E 

67. 

NBRTANtr(09/ll/86) 

00-18-25-5 

51 

05:11-22S 

152:26-52E 

68 . 

MOLUCAtr(09/ll/86) 

17-54-03-0 

33 

00:28-44N 

125:51-72E 

69. 

HrNDKUtr(09/13/86) 

14-14-54-3 

200 

36:26-52N 

70:45-96E 

70. 

SOLMONzc(09/I4/86) 

20-58-23-1 

59 

06:28-20S 

154:54-96E 

71. 

ANDREOtr(09/l5/86) 

06-29-35-8 

33 

51:33-66N 

177:05-low 

72. 

MARINAtr(09/16/86) 

18-20-17-7 

48 

19:22-56N 

146:18-06E 

73. 

AFUSSRzc(09/17/86) 

12-08-09-4 

120 

37:17-40N 

71:43-80E 

74. 

AFUSSRtr(09/17/86) 

12-08-09-4 

120 

37:17-40N 

71:43-80E 

75, 

CARBRGnD(09/17/86) 

21-25-15-0 

10 

10:29-82N 

56:59-98E 

76. 

HONSHUzc(09/20/86) 

03-04-58-7 

75 

36:27-24N 

140:31-20E 

77. 

BANDAzc(09/22/86) 

16-15-05-8 

123 

06-.42-96S 

130-.25-62E 

78. 

BONINzc(09/27/86) 

19-38-42-1 

30 

27-.52-02N 

142:46-62E 

79. 

BANDASzc(10/01/86) 

13-00-06-3 

343 

05:34-74S 

128:37-02E 

80. 

MEDRNSpk(10/02/86) 

10-12-45-6 

46 

34:50-76N 

28:18-84W 

81. 

H A LMRAtr( 10/04/86) 

02-00-08-0 

106 

02:58-92N 

128:02-16E 

82. 

ANDREOzc( 10/06/86) 

04-21-46-7 

43 

51:51-78N 

176:15-42W 

83. 

SHONSUtr(l0/07/86) 

11-40-55-0 

400 

31;54-42N 

137:40-56E 

84. 

JAVARzc(10/10/86) 

17-48-24-6 

82 

07:29-88S 

107:13-98E 

85. 

HINDKUtr(10/13/86) 

16-11-40-4 

117 

36:04-02N 

70:50-88E 

86 . 

HINDKUzc(10/I3/86) 

16-11-40-4 

117 

36:04-02N 

70:50-88E 

87. 

HONSHUtr{10/13/86) 

21-17-50-9 

65 

37:06-12N 

141:00-66E 

88 . 

BANDASzc(10/17/86) 

07-32-51-3 

67 

05:16-32S 

131:25-92E 

89. 

ANDREOpk(10/18/86) 

01-02-52-1 

33 

51:43-80N 

175:17-10W 

90. 

JAVAzc(10/18/86) 

22-09-31-7 

643 

05:37-86S 

109:59-82E 

91. 

HALMRApk(10/21/86) 

21-57-59-4 

142 

01:32-46N 

127:25-44E 

92. 

SACRUZtr( 10/22/86) 

08-59-28-8 

165 

10:34-148 

166:02-40E 

93. 

EPGNEAzc{10/23/86) 

03-54-20-8 

127 

06:05-828 

146:18-36E 

94. 

NCHILEpk{10/24/86) 

02-42-51-6 

51 

25:19-148 

70:10-56W 

95. 

NCHILEzc( 10/24/86) 

02-42-51-6 

51 

25:19-148 

70:10-56W 

96. 

AIRISEpk(10/28/86) 

15-11-23-3 

10 

30:28-268 

60:10-92E 

97. 

MINDNOfb(10/29/86) 

20-11-39-7 

72 

05:43-32N 

125:19-86E 



Table 2. ( continued) 


SI. 

No. 

Source region 

Origin time 

Depth 

(km) 

Coordinates 

98. 

MINDNOzc(10/29/86) 

20-11-39-7 

72 

05;43-32N 

125:19-86E 

99. 

BURMAzc(11/01/86) 

0502-42-4 

26 

26:54-12N 

96:25-50E 

100 . 

MINDNOzc(l 1/06/86) 

02-48-23-6 

62 

09;00-42N 

126:14-40E 

101 . 

ANDREOzc(l 1/06/86) 

18-27001 

33 

51:28-14N 

176:41 •04W 

102 . 

MOLUCAfb(l 1/11/86) 

0002-320 

79 

02;24-84N 

126:48-66E 

103. 

MOLUCAzc(l 1/11/86) 

00 02-32 0 

79 

02:24-84N 

126:48-66E 

104. 

TAKXJGfb( 11/12/86) 

100615-3 

115 

38.-27-30N 

73.-16-68E 

105. 

TAKXJGzc(l 1/12/86) 

1006-15-3 

115 

38:27-30N 

73:16-68E 

106. 

TAIWANtr(l 1/14/86) 

21-2010-5 

34 

23:54-06N 

121:34-44E 

107. 

TAIWANzc(l 1/14/86) 

23-04-370 

33 

23:51-9bN 

121:42-66E 

108. 

PRUECQfb(l 1/23/86) 

01-39-23-9 

106 

03;20-52S 

77:24-66W 

109. 

PRUECQzc( 11/23/86) 

01-39-23-9 

106 

03:20-52S 

77:24-66W 

no. 

HJAPANzc(l 1/28/86) 

22-29-35 1 

41 

36:20-70N 

141:10-44E 

Ill. 

H PA JANpk( 11/30/86) 

20I5-30-3 

37 

38:51-18N 

141:56-88E 

112 . 

NIBRISpk(12/07/86) 

05-40-39-2 

215 

06:46-98N 

95:07-20E 

113. 

SOLMrSfb(12/ll/86) 

I9-56125 

61 

10:29-28S 

I60:42-90E 

114. 

HJAPANfb(12/14/86) 

0011-31-2 

352 

32:43-62N 

137:43-38E 

115. 

AFUSSRzc(12/17/86) 

08-31-30-3 

225 

36:32-34N 

71:07-50E 

116. 

AFUSSRzc( 12/18/86) 

07-23-30-0 

182 

36.-34-08N 

71.-21-18E 

117. 

ANDREO(b(12/19/86) 

13-5010-3 

33 

51:31-02N 

176:58-62W 

118. 

ANDREOtr(l 2/19/86) 

13-50-10-3 

33 

51:31-02N 

176:58-62W 

119. 

TAIWANzc( 12/21/86) 

02-42-49-0 

275 

25:31-02N 

121:30-42E 

120 . 

MARINAtr(01/01/87) 

16-25-341 

74 

02:42-908 

138:21-78E 

121 . 

MARINAzc(01/01/87) 

I6-25-341 

74 

02:42-90S 

138:21-78E 

122 . 

. CHINAzc(01/05/87) 

22-52-46-5 

17 

41:57-84N 

81:19-14E 

123. 

TAIWANfb(01/06/87) 

0507-48-1 

38 

23:58-56N 

121:43-74E 

124. 

TAIWANzc(0I/06/87) 

0507-48 1 

.38 

23:58-56N 

121:43-74E 

125. 

NIRELRfb(01/08/87) 

19-48-55-4 

44 

04:44-16S 

153:06-30E 

126. 

NIRELRzc(01/08/87) 

19-48-554 

44 

04:44-16S 

153:06-30E 

127. 

HONSHUzc(01/09/87) 

06-14-44-8 

68 

39:53-70N 

141:40-62E 

128. 

HOKKAIzc(0l/14/87) 

11 03-48-7 

102 

42:33-90N 

142:51-OOE 

129. 

KAMCHTzc{01/19/87) 

06-47-43-0 

42 

54:44-64N 

163:16-62E 

130. 

NEPALzc(01/19/87) 

07-46-24-4 

33 

28:23-10N 

83:40-92E 

131. 

NEPALzc(01/19/87) 

08-12-05-8 

33 

28:14-58N 

83:34-32E 

132. 

MARINAzc(01/21/87) 

11-26-36-6 

118 

20:36-48N 

144:53-lOE 

133. 

MOLUCAfb(01/23/87) 

17-51 •09-2 

72 

01:38-76N 

126:31-86E 

134. 

MOLUCAzc(01/23/87) 

17-5109-2 

72 

01:38-76N 

I26:31-86E 

135. 

HINDKUtr(05/05/87) 

15-40-47-5 

208 

36:28-80N 

70:40-38E 

136. 

KAZAKHzc(05/06/87) 

04-02 05-6 

1 

49:49-80N 

78:07-50E 

137. 

ANDREOzc(05/06/87) 

0406-14-1 

20 

51:16-32N 

179:53-88W 

138. 

NBRITNzc(05/06/87) 

12-39-49-1 

20 

05:42-90S 

152:39-36E 

139. 

EUSSRzc(05/07/87) 

03 05-49-1 

430 

46:44-16N 

139:13-92E 

140. 

BALIStr(05/10/87) 

00-37-100 

42 

07:44-22S 

115:59-82E 

141. 

TALAUDzc(05/ll/87) 

09-59-34-1 

94. 

04:28-20N 

127:42-54E 

142. 

MINDNOzc(05/12/87) 

01-30-25-0 

25 

07:05-40N 

126:42-06E 

143. 

ICELNDzc(05/25/87) 

11-31-54-3 

8 

63:51-00N 

19:43-68W 

144. 

BANDASzc(05/30/87) 

16-54-04-7 

138 

06:03-848 

130:31-08E 


Table 2. {continued) 


SI. 

No. 

Source region 

Origin time 

Depth 

(km) 

Coordinates 

145. 

KURILE2c{05/30/87) 

17-19-00-1 

53 

44:40-14N 

150:17-34E 

146. 

MINHSAzc{05/31/87) 

18-32-17-0 

79 

00:44-88N 

121:56-46E 

147. 

ANDREOzc(06/01/87) 

00'15-14'3 

33 

51;32-40N 

177:30-54Y 

148. 

CHINAzc(06/05/87) 

04'59-58-3 

1 

41:35-04N 

88.-44-22E 

149. 

P.HILIPtr(06/05/87) 

21-25-IE2 

45 

05:22-86N 

127;32-04E 

150. 

PHILIPtr{06/05/87) 

22'00-03-l 

47 

05:15-90N 

127:30.78E 

151. 

PHILIP2c(06/07/87) 

05-49-43-6 

15 

20:25-74N 

121;21-96E 

152. 

KURILEfb(06/13/87) 

14-00-39-3 

42 

44:40-26N 

150;23-52E 

153. 

TALAUDtr(06/15/87) 

06-31 •44-7 

136 

03:54-94N 

125.-56-94E 

154. 

KAZAKHzc(06/20/87) 

00-53-04-8 

1 

49:54-78N 

78:44-10E 

155. 

IRIANtr(06/27/87) 

00-17-04-6 

21 

02:09-848 

138:10-20E 

156. 

KYUSHUzc(07/03/87) 

10-10-43-7 

168 

31:11-76N 

130:19-32E 

157. 

CHAGOStr{07/03/87) 

18-03-59-7 

26 

06.-47-16S 

72:14-82E 

158. 

ALEUTNpk(07/05/87) 

09-23 OO-O 

33 

51:29-16N 

174:39-60E 

159. 

ALEUTNpk(07/06/87) 

00-03-25-6 

33 

51:30-48N 

174:43-26E 

160. 

TALAUDpk(07/07/87) 

10-36-03-7 

70 

05:24-18N 

125:01-68E 

161. 

KURILEzc(07/08/87) 

22-56-02-7 

152 

46;26-22N 

149;33-48E 

162. 

TALAUDzc(07/14/87) 

04-31-23-0 

33 

03:50-04N 

126:32-52E 

163. 

OKHTSKzc(07/14/87) 

23-46-03-5 

591 

49;37-86N 

147;49-68E 

164. 

HONSHUtr(07/16/87) 

05-46-29-5 

310 

33:03-54N 

138;05-76E 

165. 

KAZAKHpk(07/17/87) 

01-17-07-0 

1 

49:47-94N 

78:06-60E 

166. 

TALAUDtr(07/22/87) 

08-03-17-2 

51 

04:03-06N 

125:34-92E 

167. 

MINHASn3(07/31/87) 

00-27-32-0 

167 

00:10-50N 

123:36-30E 

168. 

KAZAKHzc(08/02/87) 

00-58-06-8 

1 

49:52-80N 

78:55-02E 


Note; Code after event name indicates pick: fb, first break; zc, zero crossing; pk pea 
tr, trough. There are a few multiple readings in the data set. 



SIGNAL PICKS FOR ANALYSIS 


Figure 2. Seismograms showing how arrival-time are made. 







SIGNAL ATTENUATION ACROSS THE 
ARRAY 


Figure 3. Signal attenuation and signal-shape changes observed at some of the stations. 


Station because of its location on Precambrian granite outside the DVP. A comparison 
of residuals at HYB with those at Srisailam—a station also on Precambrian granite 
about 150 km south-southeast of HYB—supports our assumptions that stations on 
the Precambrian granite are “normal” compared to those over the Deccan basalts. 

For each station, the data were grouped into four categories according to epicentral 
distances: all distances, closer than 30°, 30° to 70° and beyond 70°. We first examine 
the data from Poona (POO) where the largest data set is available (figure 4). Note 
that the range of variation of residuals is about 2 s varying from 1 s in the north to 
— 1 s in the east (figure 4a). The relative residuals also depend on the epicentral distance 
and which determines the angle of emergence at the seismic station. Observe that for 
distances <30° the relative residuals are mostly positive (about +1 s) for all azimuths 
(figure 4b), with strong azimuthal variation occurring in the 30° to 70° range 
(figure 4c) and a less pronounced but still clear variation occurring for distances > 70° 
(figure 4d). Comparison of the relative residuals along the east-southeasterly profile 
from BHA, the westernmost station, to GUL on the eastern margin of DVP shows 
the pattern of residuals at BHA to be quite similar to that at POO but with larger 
variations from -)-1-5 s to — 1-8 s (figure 5). Relative residuals at UJN, SOP and GUL 
show a similar azimuthal variation, but with a progressive decrease in magnitude 
towards HYB (figure 5). GUL, closest to HYB, shows the smallest residuals and is a 
key station to substantiate our argument that the strong azimuthal variation observed 
at other stations are associated with some anomaly in the deep structure beneath the 
DVP. The relative residuals along the southeasterly profile from DHL to AUB, BIR 
and LTR (figure 6) also show a decreasing azimuthal variation towards the eastern 
margin of the DVP, though LTR has five strong negative residuals in the 40-100° 
azimuth range and 30-70° distance range, cautioning against simplified inferences. 
Note that the residuals at these stations for distances < 30° are generally negative, 






AZIMUTH (DEGREES) AZIMUTH (DEGREES) 

Figure 4. Relative residuals at POO with respect to HYB plotted as a function of azi 
of arrival, (a) All distances; hollow squares represent events at distances < 30°; solid dian 
represent events in the distance range of 30-70°; hollow diamonds represent ever 
distances >70°. (b) Events at distances <30°. (c) Events in the distance range 3C 
(d) Events at distances >70°. 

unlike those observed at other stations. The southerly profile, comprised of stat 
SAT, KOP and DWR (figure 7) shows a behaviour similar to BHA and POO (fig 
4 and 5). One surprise is at the southernmost station HPT (figure 7), south of C 
(figure 5) but outside the DVP. This station behaves similar to the western stat 
implying that the velocity anomaly responsible for the observed azimuthal variat 
extends even outside the DVP to the southeast. 

Note the abrupt transition of residuals from negative values to positive near 
100° azimuth range observed at all stations in the experimental area. In order to 
whether this feature in the azimuthal pattern of residuals could be due to a pos 
residual bias at HYB, we examined the relative residual patterns using the ave 
event residual as reference and found that no such bias exists. The relative resid 
thus calculated retained this feature at all these stations including Hyderabad, impl 
that it could not have arisen from any fortuitous cause related to the choic 
Hyderabad as a reference station. We feel that this sharp transition in the resi 
patterns reflects a sharp velocity transition in the upper mantle outside the re 
being modelled. 
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Figure 7. Relative residuals at SAT, KOP, DWR and HPT. 


7. Analysis using the simple ray plotting technique 

In this section we use simple ray-plotting techniques to show that the azimutha 
spatial variation of the relative residuals discussed above can be explained oi 
terms of velocity perturbations in the upper mantle of the DVP region. But fii 
evaluate the possible contributions that surficial and crustal layers can make 1 
residuals. 

Variation in the thickness of basalt cover which attains a maximum of about 
in the western part, can given rise to a maximum relative residual of +0-04 s asst 
the P-wave velocity of the Deccan basalts and the basement rock to be 4-8 
and 6-Okm/s respectively (Kaila 1981; Kaila 1982). Further, even a 1-0 km 
low-velocity sedimentary layer with a velocity of 4-Okm/s underlying the D 
Traps, such as that found by Kaila et al (1981) in the Saurashtra Peninsuh 
contribute a residual of less than 0-1 s. This is even lower than the noise level i 
experiment. 

7.1 Effect of the variation of crustal thickness 

The sub-trappean crustal model in the Koyna area of Deccan basalts has 
approximated by three layers of 10, 15 and 15 km thickness with average vek 
of 5-7, 6-5 and 7-Okm/s respectively (Kaila etal 1982). Variation in the thickn 
any of these layers by as much as 5 km will give rise to a total residual of only 
OT s. Crustal models of the Koyna region show that its thickness varies from 31 
near the coast to 39 km near Koyna (Kaila 1982). Even if we take the velocity co 



across the Moho to be bout 'l*Okm/s this order of variation in crustal thickness will 
only cause the arrivals to be 0-2 s earlier at the western stations compared with those 
at the eastern stations. The effect is approximately the same if the thinning of the 
crust is distributed over all the three layers. 

It is thus clear that possible variations in the thicknesses of basalts, and sub-trappean 
sedimentary layers, if present, as well as those of the upper or of the whole crust can 
at the most account for only + 0-2 s of the observed spatial variation of the relative 
residuals. The azimuthal variation cannot be accounted for by such shallow features 
at all. We therefore conclude that most of the observed residual pattern can only be 
accounted for by deep upper mantle features. 

7.2 Ray plotting analysis 

To delineate qualitatively the zones of anomalous velocity structure responsible for 
the large azimuthally varying observed residuals, we first used a simple ray-tracing 
procedure similar to that used by Iyer et n/ (1981). From a knowledge of the emergent 
angles of seismic rays at various stations, inferred from the depth and epicentral 
distances of the respective teleseisms, rays arc projected backwards from the various 
stations for each teleseism. The corresponding relative residuals are then projected on 
flat layers at three different depths (100, 200 and 300 km). If the observed residuals are 
due to the presence of an anomalous body rather than a random velocity perturbation, 
the whole set of projected residuals would be expected to show coherent signs and 
magnitudes so that they can be contoured when the depth chosen for the projection 
is correct. After considerable experimentation using data corresponding to A > 30°, 
we found that the maximum coherence was achieved for projection at a depth of 
200 km. Also, some form of subjective averaging was needed to reduce too many 
overlapping numbers. Figure 8 shows the result of this projection to 200 km depth 
using the above constraints. The main inference from this figure is that it outlines, 
albeit broadly, a 200-km wide north-south trending region of negative residuals of 
magnitude >0-5s extending from AUB in the north and probably continuing south 
beyond HPT. Since better constraints on the region of anomalous high velocity is 
provided by inversion of residuals which is presented below, we do not want to 
overinterpret this rather crude modelling of our data. However, we would like to 
state that this analysis, though not providing adequate constraints to delineate the 
anomalous structure beneath the Deccan Traps, suggests that the observed residuals 
are best interpreted in terms of a heterogeneous structure in the upper mantle 
consisting of higher-than-normal velocities beneath the south-central part of the DVP 
and lower-than-normal velocities in the west-central part. 


8. Tomographic results 

Ideally, to obtain well-constrained tomographic models, we need a set of teleseismic 
arrivals from events that are fairly well-distributed in various azimuths and distance 
ranges. Further, the arrival-time readings should be of high quality with minimal 
scatter, and the data set should be homogeneous and simultaneous at all stations 
in the seismic network. Though our experiment does not satisfy any of these criteria 
strictly, the similarity of inverted models for varying block and layer configurations. 




Figure 8. Selected averaged relative residuals projected at 200-km depth. 


suggests that they correctly resolve first-order features of the anomalous three- 
dimensional P-velocity structure beneath the DVP. Also, these models agree with 
the qualitative inference derived from the azimuthal patterns of residuals using the 
ray-plotting approach discussed above. We performed numerous ACH inversions for 
a variety of layered models. The dimensions of the block were 200 km in the east-west 
direction and 150 km in the north-south direction. We use mantle velocities from 
Herrin’s (1968b) earth model and the crustal velocity of 645km/s (Kaila 1982). For 
each inversion the respective resolution matrix and standard errors were studied with 
a view to assessing the figure of merit of the solution. Note that ACH inversions use 
relative residuals based on the average of all available readings over the network for 
each event. Hence no specific reference stations were used. 



9. Models 


Although a large number of models were used for inversion, description of a few of 
these models together with the corresponding data variance improvements are given 
in table 3. The main feature of the models, that is an anomalous high-velocity zone, 
appears to be significant only at sub-crustal depths. In order to constrain the upper 
boundary of the anomalous zone, we first analyse a series of two-layer models using all 
available data. These are shown in figure 9. The numbers indicate velocity perturb¬ 
ations, positive numbers denoting higher-than-normal velocities and negative numbers 


Table 3. Description of models. 


No. of Layers 

Layer thicknesses 
(km) 

Variance improvement 
(%) 

2 

40, 60 

16-0 

2 

40, 160 

27-0 

2 

40 260 

300 

2 

40, 360 

31-0 

3 

too, 300, 200 

53-0 

3 

150, 250, 200 

63-0 

4 

40, 160, 200, 200 

46-0 




AOO Km 



MODEL A : iO- ^00 Km 



Figure 9. Second layer of 2-layer tomographic models derived using data at all available 
distances. Numbers indicate percentage velocity perturbations. 
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indicating lower-than-normal velocities. Note that all these models invariably show 
the presence of a high velocity region within the DVP surrounded by a relatively low 
velocity region. The magnitude and lateral extent of the velocity perturbations vary 
with the model. Although the anomalies are reflected in all the models presented 
below, the first one with a layer thickness of 40-100 km shows only a small variance 
improvement (16%) compared with others whose variance improvement ranges from 
26-32%, suggesting that the high velocity anomaly (upto 4%) is confined to a depth 
range of 100-400 km (table 3). To further constrain the anomalous domain we then 
analysed a large set of three-layer models. Two of these (with variance improvements 
53% and 63% respectively) are shown in figure 10, which show that the high velocity 
anomaly decreases both in magnitude and areal extent for the layer 0-100 km 
compared with that for the layer 0-150 km. Significantly, the anomaly remains 
unaltered for depths of 400 km. In order to further test that the anomalies exposed 
by the above models are not an artefact of the choice of the number of layers, we 
also tested a suite of 4 layered models for inverting the travel time residuals. 
Figure 11 and table 4 show the model with various layers lying between 0 and 40 km, 40 
and 200 km, 200 and 400 km and 400 and 600 km (variance improvement 46%). These 
models show that the high velocity anomaly is confined to layers 2 and 3 in the depth 
range of 40-400 km. As the depth to the top of the anomaly is reasonably 
well-constrained by the models (figure 10) described earlier we conclude that the high 
velocity anomaly coinciding with the southern two-thirds of the DVP exists in the 
depth range of 100-400 km. The 400-600 km layer is characterized by the presence 
of a very diffused low velocity and is distinct from the overlying layers. Table 4 shows 
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Figure 10b. Velocity perturbations in the three layers (0-150, 150-400, 400-600 km). 
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Table 4. Full three-dimensional inversion for the four-layer model. 
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Table 4. (continued) 
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layer model. 


I. Interpretation 

icistence of large-scale velocity anomalies in the upper mantle no longer evokes any 
irprise. Regional studies of one-dimensional velocity structure using travel-times, 
nthetic seismograms, and surface-wave dispersion, show that heterogeneous 
ructure of the upper mantle is a rule rather than an exception in almost any region 
here observations are available. Teleseismic-residual studies show that velocity 
lomalies with dimensions of a few tens to a few hundreds of kilometers exist in the 
ust and upper mantle in North America and also in many other regions of the 
orld (Aki 1982; Iyer 1984, 1988a, b; Thurber and Aki 1987; Iyer and Hitchcock 
)88). In general, low-velocity anomalies are easier to interpret than high-velocity 
lomalies as the former seem to be generally associated with regions of young 
)lcanism or recent orogeny. Some typical examples are: the Yellowstone hot spot 
yer et al 1981); Rio Grande Rift (Parker et al 1984); Imperial Valley Spreading Centre 
lumphreys etal 1984); Eastern U.S. (Taylor and Toksoz 1979); the Rhinegraben 
Laikes and Bonjer 1983); and the Tarbela Area of Lesser Himalaya (Menke 1977). 
igher-than-normal velocities, invariably found in active subduction zones and 
ileo-subduction zones are also reasonably well understood as they delineate 
ibducting plates (for e.g., Japan: Hirahara 1981; Carapathian Arc: Romania: Oncescu 
al 1984). Other high-velocity anomalies, not so well understood, are found in the 
Dper mantle beneath the Baltic Shield (Husebye and Hovland 1982); Eastern Senegal, 
^est Africa (Dorbath etal 1983); the Alps (Baer 1980); the North China Basin 
ihedlock and Roecker 1987); the Timber Mountain Caldera Complex, Nevada (Evans 
id Oliver 1987); and the Southern Sierra Nevada, California (Humphreys etal 1984; 
umphreys 1987), to name a few. Humphreys (1987) interprets the vertical, tongue- 
laped, high-velocity anomaly extending to depths of 300-40C km in Sierra Nevada, 
i representing a descending sub-crustal lithosphere caused by small-scale convection, 
vailable case histories on large-scale high-velocity anomalies, though similar in size 
id velocity contrast, are simply too few to allow any detailed comparisons with the 
lomaly which we have delineated beneath the DVP. 

A fundamental question to be addressed regarding the DVP anomaly is the 
illowing. Is the anomaly at all related to the Deccan volcanism or is it simply an 
itegral part of the architecture of the upper mantle of the Indian shield? There is 
lOugh evidence from global and continental scale geophysical data to show that old 
intinents are indeed far from being homogeneous in their deep crust and upper 
lantle structure. Our modelling results lead to startling revelations indicating the 
resence of an anomalously thick high-velocity region beneath the DVP. The presence 
F such an anomaly at depths much greater than the conventionally accepted plate 
lickness (of about 100-200 km) could be best explained in the light of the concept 
f continental tectosphere (Jordan 1975,1979a, b) and stabilization of cratons (Pollack 
?86). The existence of lithospheric roots to depths as great as 400 km has been 
iferred earlier in Western North America, Superior Province of Canada and South 
frica based on petrological, geothermal and seismological investigations (Grand and 
[elmberger 1984; Pollack 1986; Lam and Jordan 1987; Silver and Chan 1988). The 


(west to east and north to south) taken from the 5 mgal contour Bouguer gravity 
anomaly map of India. These indicate the presence of a long wavelength, 35 mgal 
gravity low present over the region of high-velocity anomaly. A preliminary inversion 
of the data delineates the lower density region roughly lying within the high velocity 
anomalous region in the depth range of 100-300 km, with a density contrast of about 
— 0-006 to —0-013 grams per cc (D Srinagesh, personal communication). 

The other line of evidence is from S-wave residuals at POO relative to HYB taken 
from ISC catalogues (figure 12). In spite of the poor quality of data and consequent 
high scatter, it is quite obvious that we see an azimuthal pattern broadly similar to 
that shown by P-wave residuals at POO. The S-residuals are about three times as 
large as the corresponding P-wave residuals shown alongside for comparison. Shear 
wave structure beneath the DVP delineated by Montagner (1986) using surface wave 
data also reveals the existence of a high-velocity anomalous zone in the depth range 
of 85-250 km. However, the continuation of this anomaly beyond 250 km depth is 
not observed due to data constraints in the surface wave study. 


11. Source of Deccan basalts 

Let us now consider the question of the source of Deccan basalts in the light of the 
existence of a deep continental root under this region, which must have acquired a 
coherent stable structure long before the Gondwanaland break-up considering that 
the craton came into existence earlier than at least 2-5 billion years. The formation 
of tectosphere due to the depletion of volatiles which elevates the solidus, renders the 
upper mantle region more resistant to subsequent melting. Such stable regions are 
observed to have normal shield heat flow values (Pollack 1986). This is in conformity 
with the observed normal heat flow over DVP (Gupta and Gaur 1984). The very 
presence of the tectosphere beneath the DVP would therefore inhibit the development 
of large scale perforations which are believed to have been responsible for the extrusion 
of such extensive flood basalts, implying the absence of their source region directly 
beneath the DVP. On the other hand, it is significant to note that in all the models, 
the per cent velocity anomaly beneath the westernmost part of the DVP in the depth 
range of 40-200 km is relatively less than that observed in other parts. This is also 
reflected in the results of surface wave study (Nataf etal 1984) showing a regional 
low sheaf wave velocity enveloping the entire Arabian Sea extending from the east 
coast of Africa to the west coat of India. However, it is significant to observe that 
the above low velocity zone is confined to the west of our westernmost profile 2. 

Also, recent studies of major and trace element abundances and inter-element ratios 
as well as Sr isotopic ratios in selected samples combined with extensive field 
checks in the western parts of the DVP show that the flow sections do not correspond 
to a classical layer cake-sequence. Instead, they appear to be comprised of a sequence 
of thicker and older formations in the north which become thinner in the south and 
are superceded by southward thickening of younger ones that taper out further 
southwards. This interlocking sequence strongly suggests that the flows are derived 
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Nd and Sr isotopic values of least contaminated Deccan Thoieiites from different 
areas of the DVP fall in a restricted field of high which closely approach those 
of old and present-day Central Indian ridge basalts. This strongly suggests that the 
Deccan basalts were derived from the same large region of isotopically homogeneous 
oceanic mantle (Mahoney 1988) which feeds the Central Indian ridge. Considering 
further that the bulk of Deccan basalts were formed about 65 Ma ago shortly after 
the rifting between India and Madagascar (80 Ma) when the ridge axis jumped towards 
India (Norton and Sclater 1979; Schlich 1982), a more natural channel for their 
outpouring would appear to be associated with the rift systems that eventually 
fashioned the Indian west coast passive continental margin. If the above visualization 
of the formation of Deccan basalts is correct, we should expect to see a low velocity 
zone (LVZ) on the western periphery of the DVP representing the source of the 
Deccan basalts perhaps characterized by high heat flow values but no thermal anomaly 
over the DVP. While the heat flow over the LVZ still remains to be measured, the 
other two predictions from the model are indeed confirmed by observational evidence. 


12. Conclusions 

Using teleseismic P-wave residual tomography, we have identified a large region in 
the upper mantle beneath the Deccan volcanic province in which the P-wave velocity 
is higher than that in the surrounding region by 1 to 4%. This high velocity anomaly 
exists coherently in the depth range of 100-400 km. This implies that a coherent, 
colder and comparatively rigid lithospheric root extends to a depth of about 400 km 
beneath the DVP. The anomaly in the depth range of 400-600 km is distinctly weaker 
and appears to be decoupled from the overlying layers. This is as we expect owing 
to diminished ability of the asthenosphere to retain heterogeneities over extended 
periods of geologic times. Persistence of the anomaly in the upper mantle is however 
maintained through a stable geochemical reordering and formation of a tectosphere 
beneath the craton which together translate coherently over the underlying weaker 
mantle. 
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Abstract. The four major earthquakes that have occurred in the Himalayan region since 
1897 seem to have ruptured as little as about 15% or 20% to perhaps as much as 45% of the 
thrust zone separating the underthrusting Indian Shield and the overthrusting Himalayan 
crystalline nappes. Because of various difllculties in estimating the rupture zones for each of 
these earthquakes, we cannot place a tight constraint on the fraction of the Himalayan belt 
for which the risk of an imminent great earthquake is high. If a slip between the Indian Shield 
and the Himalayan crystalline nappes occurs largely by slip associated with major 
earthquakes, then recurrence intervals of such earthquakes are likely to be between 200 and 
500 years, with a likely value of 300 years. 
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1. Introduction 

Four great earthquakes have occurred in the Himalayan region in the last 100 years. 
The 1897 Assam earthquake seems to have ruptured a gently dipping plane beneath 
the Shillong Plateau (Oldham 1899; Seeber and Armbruster 1981; Molnar 1987b), and 
the rupture may have extended north beneath the Himalaya. The earthquakes of 1905, 
1934 and 1950 occurred within the Himalaya, and were responsible for extensive 
damage along the Himalaya in zones up to 300 km in length. Earthquakes of this type 
almost surely constitute the most serious seismic hazard to the residents of the 
Himalaya and the neighbourhood. 

Several aspects of the structure and tectonics suggest that major earthquakes in the 
Himalaya occur as the Indian shield is thrust beneath the Himalaya (Seeber and 
Armbruster 1981). Fault plane solutions of nearly all moderate earthquakes 
(5-5 ^ M ^ 7) in the Himalaya indicate thrust faulting with one nodal plane dipping 
gently north or northeast and with the other nodal plane striking parallel to the range 
and dipping steeply south or southwest (figure 1) (e.g. Fitch 1970; Rastogi 1974; 
Molnar et al 1973, 1977; Chandra 1978; Baranowski et al 1984; Ni and Barazangi 
1984). The gently northward or northeastward dipping planes almost surely are the 
fault planes. Focal depths of about 15 km for these earthquakes are consistent with 
their occurring on the top surface of the Indian shield as it underthrusts the Himalaya 
(e.g. Molnar and Chen 1982; Baranowski et al 1984; Ni and Barazangi 1984). Gravity 
anomalies over the Indo-Gangetic Plains and the Himalaya are consistent with the 
Indian plate extending intact northward from the Indo-Gangetic Plains 



Figure I. Rupture zone of great Himalayan earthquakes and fault plane solution of moderate earthquakes. Shaded areas show the approximate 
areas underlain by the rupture zones of the four major earthquakes in the Himalaya since 1897. Lower hemisphere projections of the focal 
spheres, with blackened quadrants showing those with compressional first motions, were taken from Molnar et al (1977), Baranowski e£ al (1984) 
and Ni and Barazangi (1984). 



Himalaya, with an abrupt increase in elevation along the southern margin of the 
range, the mature landscape of the Lesser Himalaya, and the sharp incision of the 
Greater Himalaya, is consistent with the Indian plate underthrusting the Lesser 
Himalaya along a gentle plane that becomes steeper beneath the Greater Himalaya 
(Molnar 1987c). Thus, several observations suggest that the Indian shield is being 
underthrust coherently beneath the Lesser Himalaya along a gently dipping fault that 
steepens northward. 

None of the four great earthquakes seems to be associated with primary surface 
faulting, and therefore the types of faults that ruptured cannot be proven definitively. 
Seeber and Armbruster (1981) suggested that the great Himalayan earthquakes 
occurred on a gently dipping thrust fault that underlies the entire Himalaya and part 
of the Indo-Gangetic Plains. We proceed with the assumption that this contention is 
basically correct, but we disagree with their estimates of the dimensions of the rupture 
zones, and in particular, with their inference that the ruptures extend beneath the 
Indo-Gangetic Plains. 

We review briefly the constraints on the rupture zones of the four great Himalayan 
earthquakes of 1897, 1905, 1934 and 1950, and we also report some tentative 
implications for recurrence intervals of great earthquakes in the Himalayan region. 

2. Rupture zones of great Himalayan earthquakes 

2.1 The 1897 Assam earthquake 

Virtually everything that is known about the effects of this earthquake is given in 
Oldham’s (1899) report. Later, others have extracted factual material to draw 
conclusions about the style and extent of faulting associated with that earthquake (e.g. 
Richter 1958; Seeber and Armbruster 1981; Molnar 1987b). 

Very high accelerations and particularly great destruction were noted at numerous 
localities on the Shillong Plateau, and there seems little doubt that the plateau is close 
to the rupture zone. Indisputable evidence of surface faulting was found only on the 
Shillong Plateau, but the short segments and varying orientations of such faults 
suggest that this faulting was secondary, not primary. Warping and tilting of the 
earth’s surface, as reflected by changes in drainage indicated important strains 
associated with the earthquake. Finally, a nearly continuous belt of landslides along 
the southern edge of the plateau suggested high accelerations there. Oldham (1899) 
inferred that a slip occurred on a thrust fault that dips gently northward from the 
southern edge of the Shillong Plateau, an inference that he later abandoned, but that 
seems reasonable (e.g. Seeber and Armbruster 1981; Molnar 1987b). 

The likelihood of such a thrust fault cannot be denied, but its proof is lacking. The 
east-west striking Dauki fault, which bounds the southern edge of the Shillong 
Plateau, marks a vertical separation of the Pre-Tertiary basement of more than 10 km 
(e.g. Evans 1964; Sengupta 1966). Although Evans (1964) inferred a large component 
of strike-slip faulting, thrust faulting can account for most of the observations that he 
presented. The very large isostatic gravity anomalies of the Shillong Plateau (e.g. 
Verma et al 1976; Verma 1985) imply that the strength of the Indian lithosphere must 
support the plateau, and therefore underthrusting of the strong Indian plate beneath 
the Shillong Plateau seems likely. Neither the locations of microearthquakes (Khattri 
etal 1983; Kayal 1987) nor fault plane solutions and locations of moderate 



earthquakes (Chen and Molnar 1989) give a hint oi sucn a lauit. me live moaerate 
earthquakes studied by Chen and Molnar (1989) occurred at depths between 25 and 
50 km, too deep to have occurred on a gently dipping plane that crops out as the 
Dauki fault, and four of the fault plane solutions show faulting on relatively steep 
places (30°-60°). Thus, although Oldham’s (1899) inference of a gently dipping plane 
seems reasonable, it is not proven. 

The extent of the rupture remains in dispute. Damage associated with the 1897 
earthquake was not confined to the Shillong Plateau. Oldham (1899) noted that much 
of the damage, particularly west and south of the plateau, could be associated with 
fissuring, and he specially noted that the soil conditions associated with and without 
fissures were different. Thus, he considered this damage to be due to superficial 
deformation and not to faulting of the underlying basement. One of us, Molnar 
(1987b), found his arguments convincing, but Seeber and Armbruster (1981) inferred 
that the rupture extended 100 km south and 200 km west of the plateau. The 
numerous reports of aftershocks felt on the plateau and just north of it and their 
paucity from areas to the south or west (Oldham 1899) also support the contention 
that the fault that ruptured underlies the plateau but not the idea that this rupture 
zone extends south and west of the plateau. 

The eastern limit of the rupture is poorly defined, largely because Oldham and his 
deputies were unable to gather much data from that area. From the few observations 
Oldham (1899) reported, Molnar (1987b) inferred an east-west extent of the rupture of 
200 ± 40 km, from the western margin of the Shillong Plateau, considerably less than 
the 550 km suggested by Seeber and Armbruster (1981). 

The northern limit of the rupture is also uncertain. Reports of landslides and 
damage in Bhutan (Oldham 1899) allow the possibility that the rupture reached as far 
north as the Himalaya, as Seeber and Armbruster (1981) suggested. In any case, for 
evaluating earthquake hazards in the Himalaya, knowing the northward extent of the 
rupture is probably less important than determining the degree to which the 1897 
earthquake relieved strain caused by the convergence of the Indian plate with the 
Himalaya. If the 1897 earthquake did relieve such strain, then obviously its occurrence 
must be considered in evaluating the earthquake hazard of the Himalaya. If, instead, 
the 1897 earthquake did not relieve such strain, and if the convergence of the Shillong 
Plateau with the Himalaya occurs essentially as rapidly as that between the rest of the 
Indian subcontinent and the Himalaya, then clearly the 1897 earthquake could be 
ignored in evaluating the earthquake hazard in the Himalaya. The available evidence 
seems to allow either possibility. 

2.2 The 1905 Kangra earthquake 

Middlemiss’s (1910) report contained virtually all the published factual material 
concerning the macroseismic effects of this earthquake. The highest intensity, X on the 
Rossi-Forel scale, was felt near the towns of Kangra and Dharamsala. The intensity 
decreased sharply to the west, with the intensity VIII and X isoseismals separated by 
only 15 km. To the east, the decrease in intensity was more gradual. The intensity VIII, 
IX and X isoseismals are all aligned northwest-southwest, parallel to the trend of the 
range and surrounding the Main Boundary Fault, the thrust fault separating Tertiary 
from pre-Tertiary rock at the foot of the Himalaya. The limited number of 
observations northeast and immediately southeast of the I = VIII isoseismal make 



occurred by slip on a gently northeastward dipping plane, but there is no proof of this 
and little evidence to constrain the northeast limit of the rupture. 

Damage associated with the 1905 earthquake was greater in DehraDun and 
Mussoori, about 250 km southeast of Kangra, than in most of the area southwest and 
northeast of these two towns. Accordingly, Middlemiss (1910) drew an isolated 
contour of intensity VIII surrounding Dehra Dun and Mussoori. This contour must 
close southeast of Dehra Dun, because very little damage was reported at Rishikesh, 
Hardwar, and neighbouring towns. Very few observations, however, define the closure 
to the northwest and probably are insufficient to rule out the possibility that the two 
separate / = VIII isoseismals’ought to be connected. 

Relevelling in 1906 along a line across the Siwalik Hills to Dehra Dun, which had 
been surveyed in 1862, and its continuation to Mussoori, which had been surveyed in 
1904, revealed 4 cm of uplift. Middlemiss (1910) and Seeber and Armbruster (1981) 
associated this uplift with the 1905 earthquake, but we remain unconvinced of this 
association. The amount of only 4 cm seems small to be associated with the rupture of 
a major earthquake. Moreover, the uplifted part lies in a topographic basin flanked by 
ranges that clearly have risen in late Cenozoic time. Thus ascribing the 4 cm of uplift 
to the 1905 earthquake seems to us to be risky. 

The extent of the rupture is poorly defined, and several different possible rupture 
zones seem to be allowed (Molnar 1987a). One is that rupture occurred only beneath 
the area delimited by the I — VIII isoseismal surrounding Kangra, Dharamsala, a 
zone about 100 km in length. A second possibility is that the two separate areas 
encompassed by the time / = VIII isoseismals ruptured, for a total of about 200 km of 
the Himalayan front. Third, perhaps the entire zone from northwest of Kangra to 
southeast of Dehra Dun, 280 km in length, ruptured, as Seeber and Armbruster (1981) 
inferred. Included within the latter, two plausible rupture zones is the possibility that 
the slip on the southeastern part of the zone was less than that in the northwest, where 
the maximum intensity was greater. In evaluating the earthquake hazard, it seems 
reasonable to assume that only the segment approximately 100 km in length 
surrounding Kangra ruptured in 1905, but also to consider the possibility that the 
entire zone 280 km in length ruptured. 

2.3 The 1935 Bihar-Nepal earthquake 

Although the 1934 earthquake was recorded by more than 100 seismograph stations, 
to the best of our knowledge, there are not enough consistent P wave first motions to 
determine a fault plane solution. Singh and Gupta (1980) inferred from surface waves 
that strike-slip faulting had occurred, but they did not show their data to be 
inconsistent with thrust faulting. Chen and Molnar (1977) assumed thrust faulting and 
used the amplitude spectra of surface waves to infer a seismic moment. Because the 
inferred value of the moment is inversely proportional to the assumed dip for slip on 
gently dipping planes, however, their value is very uncertain. Most of what we know of 
this event derives from macroseismic observations. 

Two important sources of macroseismic data were published independently and 
separately from one another; Dunn et al (1939) gave extensive information on the 
effects of the earthquake in India and along three traverses into Nepal, and Rana 
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(1936) presented extensive data from the interior of Nepal. Neither Rana nor Dunn 
and his co-workers seem to have been aware of the others’ work. Thus the agreement 
of the two reports concerning areas discussed by both gives credibility to the areas 
discussed by only one of them (Pandey and Moinar 1988). 

Dunn et al (1939) reported three areas of very extensive damage to which they 
assigned an intensity of / = IX and X on the Mercalli scale: a belt along the 
India-Nepal border, a narrow zone along the Ganga river through Monghyr, and the 
Kathmandu valley in Nepal. The damage in all three areas seems to be due, at least in 
large part, to local site conditions. The extensive damage along the India-Nepal 
border is primarily due to slumping (the “slump belt”). Evidence for faulting in the 
basement is absent, and several observations imply that the slumping was superficial 
(Dunn etal 1939; Pandey and Moinar 1988). Localized high intensities had been 
noted at Monghyr in association with other important Himalayan earthquakes (Baird 
Smith 1843; Dunn etal 1939). Similarly, the contrast in damage within and 
surrounding the Kathmandu valley, both in 1934 and in 1833, when a smaller 
earthquake caused extensive damage, implies that local conditions are largely 
responsible for this localized damage. Moreover, a recent study of variations in earth 
noise in the Kathmandu valley shows marked variations that correlate with variations 
in damage in 1934 (Pandey, unpublished data, 1986). Thus, the areas assigned the 
highest intensity by Dunn et al (1939) do not seem to place useful constraints on the 
extent of the rupture. 

Following the 1934 earthquake, JB Auden made three excursions into Nepal to 
evaluate the damage: one to Kathmandu, a second into the frontal range of the 
Himalaya in eastern Nepal, and a third, a traverse across the Lesser Himalaya of 
easternmost Nepal and the frontier with India to Darjeeling. Except for the 
Kathmandu valley, Auden saw little evidence for destruction as severe as that in the 
slump belt. These observations, together with an apparent mislocation of the epicenter 
of mainshock, led to the inference that the rupture zone underlies the slump belt in the 
Indo-Gangetic Plains and not the Himalaya (Dunn et al 1939; Richter 1958; Singh 
and Gupta 1980; Seeber and Armbruster 1981). 

Ironically, Auden did not visit most of the areas of Nepal where destruction appears 
to have been most severe, and where it cannot be easily ascribed to local conditions. 
Rana (1936) reported almost complete destruction of the area north of Auden’s central 
route and west of his eastern route, an area that includes Chen and Molnar’s (1977) 
relocated epicenter of the 1934 earthquake. The severity of the destruction decreased 
eastward; Darjeeling was not affected much. To the west the level of destruction 
diminished more gradually to the area just south of Kathmandu, but west of 
Kathmandu, there was almost no destruction. Thus, the rupture zone probably 
underlies the Lesser Himalaya of eastern Nepal. The region of most severe destruction 
spans a segment parallel to the Himalaya approximately 100 km in length, but 
destruction was not negligible in the area from the eastern frontier of Nepal to just 
west of Kathmandu, about 300 km in length. Thus the length of the rupture zone is 
best described as 200+ 100 km. 

2.4 The 1950 Assam (or Chayu) earthquake 

This earthquake occurred at the northeast end of the Himalaya; although the 



epicenter lies in China, numerous aftershocks beneath the Himalaya in eastern Assam 
[e.g. Chen and Molnar 1977) imply that at least part of the rupture zone underlies the 
Himalaya. 

Tandon (1955) reported a fault plane solution with numerous compressional first 
motions at stations north, northwest, and northeast of the earthquake. He inferred 
normal faulting, either on a steeply northward or on a gently southward dipping 
plane. Ben-Menahem et al (1974) revised and added data and inferred a large 
component of right-lateral strike-slip faulting on a north-northwest trending steeply 
east-dipping plane. They used surface waves and the existence of scattered aftershocks 
south-southeast of the epicenter of the mainshock to corroborate this solution. Chen 
and Molnar (1977) found the P-wdve first motions used by Ben-Menahem et al (1974) 
to be consistent with thrust faulting on either a gently north-northwest or a steeply 
south-southeast dipping plane, similar to those of moderate earthquakes (figure 1). 
Because of the extensive aftershock activity beneath the Himalaya and the similar 
solutions for nearby moderate earthquakes (figure 1), they inferred thrust faulting on 
the gently dipping planes. This solution, however, does not agree with Ben-Menahem 
et fl/’s (1974) measured ratios of Love to Rayleigh wave amplitudes, and the possibility 
that two separate planes ruptured should not be overlooked (Molnar et al 1977; 
Molnar and Deng 1984). Finally, a re-examination of the aftershocks of the 1950 
earthquake shows that if only those earthquakes recorded by more than 50 stations 
are considered, nearly all of the relocated aftershocks lie beneath the Himalaya in a 
zone extending roughly 250 ± 50 km west of the epicenter of the mainshock and about 
100 km wide (Molnar, unpublished data, 1985). Thus, we presume that rupture 
occurred primarily on a gently north-northeast dipping thrust fault. 

The remoteness of the region afiected by the 1950 earthquake prohibited a 
comprehensive investigation of the damage associated with it. Apparently no special 
investigation of the damage in the epicentral region was made. Ramachandra Rao 
(1953) showed several maps of intensity distributions for India, and all indicate a 
maximum intensity, / = XI (or XII) on the Modified Mercalli scale, along the 
Himalaya to about 250 km west of the epicenter. The absence of data within the 
Himalaya limits the northward extent of this isoseismal. A Chinese map show, these 
isoseismals closing just east of the epicenter (Anonymous 1979), consistent with the 
rupture underlying the Himalaya, Thus, given both the aftershock distribution and 
the isoseismals, we infer that dimensions of the rupture zone are about 250 km in its 
east-west dimension of about 100 km north-south. 


3. Summary and implications 

The lengths of the rupture zones parallel to the Himalaya for all four earthquakes 
exceed 100 km; those for the events of 1897, 1934 and 1950 probably are at least 
200 km; and all may reach 300 km, but probably not more than that. Thus, if the slip 
during the 1897 Assam earthquake relieved the strain in the Himalaya and 
accommodated some of the convergence between India and the Himalaya, the strain 
along nearly half of the 2500 km of the Himalaya might have been relieved in the last 
100 years. Alternatively, if the 1897 Assam earthquake did not relieve such strain, 
perhaps strain along as little as only 15-20% of the Himalaya chain was relieved. If we 
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choose conservative, but probable, lengths for the ruptures of 200 km (1897), 100 km 
(1905), 200 km (1934), and 250 km (1950), then 30-35% of the Himalayan chain 
ruptured in the last 100 years. This large range of uncertainty (30% + 15%) for the 
fraction of the chain along which strain seems to have been relieved since 1897 
underscores our ignorance of the seismic hazard in the Himalaya. 

If (1) the convergence of the Indian plate with respect to the Himalaya occurs 
largely by slip associated with great earthquakes, if (2) these four earthquakes relieved 
all of the strain that had accumulated since the last such earthquakes in the same 
areas, and if (3) the rate of 4 great earthquakes per century were representative of 
earthquake recurrence in the Himalayan region, then recurrence intervals should be 
between about 200 and 500 years, with a most probable value of 300 years. The areas 
most likely to rupture next are those that did not rupture in 1897, 1905,1934, or 1950 
(e.g. Khattri and Tyagi 1983). 

Clearly none of the three assumptions listed above are proven, and there is reason 
to doubt each of them. Several earthquakes with M = 7-5 have occurred this century. 
They, as well as fault creep, may contribute a non-negligible fraction of the 
convergence. The possibility that different amounts of slip occurred on different parts 
of the 1905 rupture zone, and the differences in the measured seismic moments, if both 
very uncertain, for the 1934 and 1950 earthquakes imply that the amounts of slip 
associated with ruptures in different parts of the Himalaya have been different. 
Consequently, recurrence intervals for different segments also should be very different 
from one another. Finally, although major earthquakes did occur during the last 
century before 1897, reports of them (e.g. Baird Smith 1843; Oldham 1883; Dunn etal 
1939; Seeber and Armbruster 1981) do not suggest that any was as large as those in 
1897, 1905, 1934, or 1950. Thus, the last 100 years may not be representative of the 
seismic history of the Himalaya for the last few thousand years. 

Studies of the great Himalayan earthquakes clearly are vital for evaluating the 
earthquake hazard, but enough remains unknown about each event that alone they 
do not place a tight constraint on the earthquake hazard. The value of these studies 
could be enhanced by a thorough study of archived historic documents that give 
information about earthquakes prior to 1897. 

Because the geologic record is not likely to contain a complete record of 
earthquakes that did not cause surface faulting, techniques of Quaternary geology may 
not provide tight constraints on the average recurrence intervals of earthquakes along 
the Himalaya. Obtaining such constraints probably will require a synthesis of 
geologic, geophysical, and historical data of various kinds. Among such data are those 
that constrain average Quaternary, Holocene, or more recent (geodetically measured) 
slip rates on the main active thrust fault in the Himalaya would be very important in 
future estimates of recurrence intervals. 
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Active tectonics of the Himalaya 


JAMES F NI 

Departmeht of Physics, New Mexico State University, Las Cruces, NM 88{X)3, USA 

Abstract. The Himalayan mountains are a product of the collision between India and 
Eurasia which began in the Eocene. In the early stage of continental collision the 
development of a suture zone between two colliding plates took place. The continued 
convergence is accommodated along the suture zone and in the back-arc region. Further 
convergence results in intracrustal megathrust within the leading edge of the advancing 
Indian plate. In the Himalaya this stage is characterized by the intense uplift of the High 
Himalaya, the development of the Tibetan Plateau and the breaking-up of the central and 
eastern Asian continent. Although numerous models for the evolution of the Himalaya have 
been proposed, the available geological and geophysical data are consistent with an 
underthrusting model in which the Indian continental lithosphere underthrusts beneath the 
Himalaya and southern Tibet. Reflection profiles across the entire Himalaya and Tibet are 
needed to prove the existence of such underthrusting. Geodetic surveys across the High 
Himalaya are needed to determine the present state of the MCT as well as the rate of uplift 
and shortening within the Himalaya. Paleoseismicity studies are necessary to resolve the 
temporal and spatial patterns of major earthquake faulting along the segmented Himalayan 
mountains. 

Keywords. Himalayan mountains; continental collision; earthquakes; active tectonics; 
tectonic setting; Main boundary thrust; Tethyan slab. 


1. Introduction 

The Himalaya is the largest mountain chain in the world. Accordingly, understanding 
the mechanics responsible for its formation and for its support must be two of the 
most important geological problems. Moreover, study of the active tectonics of the 
Himalayan mountain belt is important because such a region contains information 
that is unavailable in other, older mountain belts, but necessary in order to 
understand how mountain building takes place. 

Early episodes in orogenic research (last quarter of the 18th and first half of the 19th 
centuries) were characterized by the interpretation that mountain belts are a result of 
vertical uplift. It was not until the late 19th century that the role of horizontal 
movements became more widely appreciated. The development of geosynclinal theory 
in the early 20th century to explain the prehistory of mountain belts and the 
subsequent horizontal shortening in orogeny as expressed by fold and nappe 
structures delayed a more accurate understanding of mountain building. As the 
concept of plate tectonics was developed during the 1960s, numerous models for the 
evolution of the Himalaya were proposed. The past decade has included a transition 
from a focus on testing of the Himalayan collision models to research based more on 


This review focuses on the deep crustal and uppermost mantle structures of the 
Himalaya, because these features reflect the mechanical processes that result in 
mountain building. The review also will discuss recent work on the Pakistan Himalaya. 


2. Collision models 

Early ideas about the architecture of the Himalaya assumed that the Himalayan 
Mountains rest on a rigid substratum that is similar to the rocks beneath the Indian 
Shield. However, geodetic measurements (plumb line) of Pratt (1855) showed that the 
observed deflection was much smaller than expected for a rigid model. Shortly 
following this experiment the concept of isostasy was born and a floating Himalayan 
range with a “crustal root” became the dominant geological thinking for the early 
twentieth century. With the improvement of gravity meters, numerous studies have 
shown that gravity measurements deviate from local isostatic equilibrium. This 
discovery suggested that the Himalaya is not completely supported by the buoyancy 
force; rather it is held up, at least in part, by an upward force. Searching for the origin 
of this force is essential to understanding the dynamics of mountain building. 

The concept of underthrusting of one continental lithosphere beneath another as 
the cause of mountain building and plateau uplift was first proposed by Argand 
(1924). With the advent of the new global tectonics (e.g. Isacks et al 1968) a number of 
different mechanisms for the uplift of the Himalaya and Tibet have been developed. 
Dewey and Bird (1970) suggested that the Himalaya was formed solely by the collision 
of an Atlantic-type continental margin with the Asian continent bordered by a 
marginal trench. In their model they assumed that the buoyancy of the continental 
lithosphere prevents further subduction and they overlooked the possible significance 
of intracontinental thrusts. Powell and Conaghan (1973) pointed out that the 
geological history of the Himalaya does not agree with the simple continent- 
continent collision model of Dewey and Bird (1970), because the collision itself did 
not directly produce the present-day Himalaya. Alternatively, Powell and Conaghan 
(1973, 1975) and Powell (1979) proposed an “evolutionary” model (figure 1) invoking 
two phases of orogeny. In their model, the first phase is related to active Mesozoic- 
early Tertiary subduction zones along the present-day Indus-Tsangpo suture, 
analogous to the North American Cordilleras. This phase ended in the Eocene (e.g. 
Molnar and Tapponnier 1975; Tapponnier and Molnar 1977) when India collided 
with Asia. The Indus-Tsangpo suture is a direct result of this collision. The second 
phase involves formation of intracontinental fractures within the northern Indian 
margin during the Middle Tertiary and underthrusting of the Indian continent along 
this fracture from Miocene time to the present. 

Since 1973, increasing amounts of data on the Himalaya and Tibet have become 
available and, consequently, many plate tectonic models have been proposed for the 
evolution of the Himalaya (e.g. LeFort 1975; Molnar et al 1977; Bird 1978; Seeber et al 
1981; Chen and Molnar 1981; Barazangi and Ni 1982; Molnar and Chen 1982; Ni and 
Barazangi 1984; Molnar 1984). Among these models, much of the debate about the 
tectonics of the Himalaya is concerned with the nature of the Main Boundary Thrust 
(MBT) and Main Central Thrust (MCT). According to LeFort (1975), the Himalaya 
was formed by the imbrication of slices of India’s northern margin. In his evolutionary 
model (figure 1), the MBT and MCT are assumed to be similar, but successive, 



A. La Fort 1975 


B. Powell 1979 


MBT MCT 



Himalayan 

Schuppen 

Zone 


. Tethyan . 
Himalaya 


Indus ^ 
Suture 



^_ Sedimentary _ 

Wedge Tethyan Himalaya 




[^Molasse [||] Cretaceous ^ Ophiolite Metamorphic 
Plutons Foliation 

Figure 1. Comparison of tectonic models of the Himalaya. A. After La Fort (1975). The 
Himalaya was formed by imbrication of slices of India’s northern margin. B. After Powell 
(1979). The MCT and MBT are fundamentally different structures in this model. C. After 
Seeber and Arrnbruster (1981). All the elements in this model are active simultaneously and 
subduction proceeds without fundamental changes in the structure. This model was named 
the “steady-state” model. D. After Burg and Chen (1984). In their model, decoupling and 
thrusting, which separate lower and upper crustal layers, occur contemporaneously with each 
other. 

tectonic thrusts. The southern thrust zone (i.e. the MBT) becomes the new boundary 
of the continental convergent zone, while the older thrust zone (i.e. the MCT) becomes 
a less active or a dormant feature. 

In contrast, Seeber and Arrnbruster (1981) proposed a “steady-state” model 
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(figure 1) for the evolution of the Himalaya. This model requires that the MBT and 
MCT are contemporaneous features. Central to this model is a detachment 
(decollement) surface underlying the entire Himalaya. The detachment represents the 
upper surface of the underthrusting Indian plate. Seeber et al (1981) also postulated 
the presence of a Basement Thrust (BT) where the northerly and relatively steeply 
dipping MCT merges with the shallow dipping detachment surface at depth (see 
figure 11 of Seeber et al 1981). The BT zone represents the front of the overriding 
crystalline basement (analogous to a bulldozer blade), and it is along this front where 
the sedimentary rocks above the Indian Shield are scraped off. 

A more complex model for the Himalaya has been proposed by Him et al (1984a) 
and Burg and Chen (1984) (figure 1). Their model is mostly based on the interpretation 
of recent reflection data in which a variation in the level of the Moho has been 
suggested. In their model, decoupling and thrusting, which separate lower and upper 
crustal layers, occur contemporaneously with each other. However, polarities of 
thrusting are opposite to each other (figure 1). 


3. Tectonic setting 

The hypotheses of the evolution of the Himalaya discussed above depend largely on 
understanding the evolution and structure of major tectonic features and their 
behaviour at depth. The following is a brief summary of the main structural 
characteristics of these features (figure 2). 

3.1 Indus-Tsangpo suture 

The Indus-Tsangpo suture (e.g. Gansser 1964,1980; LeFort 1975) marks the northern 
limit of the Indian subcontinent following the late Cretaceous-early Tertiary closure 
of the Neo-Tethys. The suture is made up of imbricated melanges of flysch sediments, 
radiolarite, pillow lavas, volcanics, and basic and ultrabasic rocks that are cut by steep 
thrust faults (e.g. Stocklin 1980; Allegre et al 1984). 

3.2 Main central thrust (MCT) 

The MCT, at the base of the crystalline zone, (e.g. Gansser 1964) dips 30° to 45° 
northward, separating the Higher Himalaya from the Lesser Himalaya. The thrust is 
characterized by a zone of intense shearing; no master thrust surface is apparent (e.g. 
Valdiya 1980; Sinha-Roy 1982). An abrupt change in the style of structures, especially 
folding, and in the grade of metamorphism is the main manifestation of the presence of 
the MCT. The MCT appears to have developed since mid-Tertiary time, and there are 
some geological indications of minor recent movements along the MCT (e.g. Valdiya 
1980; Stocklin 1980; Gansser 1982). 


3.3 Main boundary thrust (MCT) 




Figure 2. Map of Himalayan arc, Tibetan plateau, and surrounding regions showing major structural features (after Gansser 1964; Molnar and 1 
1975). MCT, main central thrust; MBT, main boundary thrust; MFT, main frontal thrust. Neogene granites appear to be related to MCT and 
to the Central Himalaya. 
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planes of the MBT at the surface gradually flatten out at depth (Valdiya 1980). The 
MBT has developed since Pliocene time, and abundant geological evidence shows 
that the thrust zone was active through the Pleistocene (e.g. Mathur and Evans 1964; 
LeFort 1975; Valdiya 1981). 

3.4 Tethyanslab 

The region between the MCT and the Indus-Tsangpo suture is composed of the 
Tethyan sedimentary zone in the north and the central crystalline zone in the south, 
and is often referred to as the Tethyan slab. The richly fossiliferous, platform-type 
sediments of this region (Paleozoic-Mesozoic age) are clearly of Gondwana affinities 
and represent the former Indian margin. Burg and Chen (1984) showed that the 
Tethyan slab consists of several thrust sheets with distinct stratigraphies, states of 
deformation and metamorphic grades. Crustal shortening of the northern Tethyan 
slab is produced by decoupling and thrusting of both upper and lower crustal layers 
(Burg and Chen 1984). As the result of continental collision the Tethyan sediments are 
also intruded by Tertiary leucogranite. The leucogranite, dated between 10-20 Ma, 
(LeFort 1975, 1981; Vidal 1978; Wang et al 1981) appears to have a certain 
stratigraphic and structural relationship when compared with the position of the 
MCT. First, the metamorphic isogrades run roughly parallel to the MCT for great 
distances (e.g. LeFort 1975). Second, the metamorphic isogrades are inverted in the 
upper level of the Lesser Himalaya. A high ®'^Sr/®^Sr ratio (e.g. LeFort 1981; Wang 
et al 1981) of the leucogranite clearly indicates that it was derived from a crustal 
source, however, the source of heat remains an enigma. 

3.5 Normal faulting in the High Himalaya 

Burg et al (1984) indicate the existence of an east-west striking, gently north-dipping 
normal fault of regional extent in the High Himalaya. This normal fault zone 
separates the highly metamorphosed crystalline rock and a transitional zone from the 
overlying late Precambrian to Cambrian sedimentary sequence. A similar type of 
normal faulting has been observed in the Zanskar region of the western Himalaya 
(Herren 1987). Movement along the shallow-dipping normal fault is believed to be 
about 20 km and limited to the deformation of upper crust. Analysis of two-mica 
granite shows that normal faulting is coeval with the south-directed overthrusting 
along the MCT (Copeland et al 1987). Burchfiel and Royden (1985) suggest that the 
tensional stress induced by the difference in topography between the Indian foreland 
and the southern edge of the Tibetan Plateau is a possible cause for the extensional 
faulting parallel to the Himalayan ranges. At present there are little data on these 
normal faults. It appears that the normal faulting can be considered as an important 
effect superimposed on the north-south regional compression. 

3.6 Strike-slip faulting 

A secondary effect of the north-south shortening of the Himalaya is the pervasive 
strike-slip faulting throughout most of the Lesser Himalaya. Valdiya (1976, 1981) 
indicated that the strike-slip faults have right-laterally displaced the MBT by as much 

nc 1 9 Vm Oftipr PYiam r\\ PC r\f otriVf»_o1ir*k foiiltinrr in Mien'll 


I iLartnquaKes ana jauit plane solutions 

The Himalaya is seismically active. During the last decade not only have four great 
earthquakes (M > 8) occurred along the Himalaya, but also a large number of 
moderate-sized earthquakes. Local seismic networks are restricted to the western 
Himalaya and the Assam region. Therefore, earthquake hypocenters for the 
Himalayan events are based on teleseismic locations. Routine location procedures 
usually yield reasonably good locations with uncertainties not exceeding 20 km. 
However, the focal depths of these events can have errors of as much as 50 km. 
Recently developed techniques of synthetic body wave modelling of P and SH waves 
not only provide a better constraint on the orientation of nodal planes but also 
provide good depth control for shallow events to ± 3 km. Using this technique 
Baranowski et al (1984) and Ni and Barazangi (1984) showed that the moderate sized 
thrust-type events in the central and eastern Himalaya lie between 10 and 15 km 
depth. Based on a qualitative selection criteria, Ni and Barazangi (1984) showed that 
the better located epicenters are concentrated in a narrow zone, about 50 km wide, 
lying just south of the MCT (figure 3). Focal mechanisms obtained from forward 
waveform modelling indicate thrust faulting with one nodal plane dipping gently 
(< 30°) north or northeast beneath the Himalaya and the other nodal plane dipping 
steeply southward (figure 3). It is generally accepted that the shallow north-dipping 
nodal planes are the fault planes (figures 3,4). 

Projection of all reliable thrust-type events along the Himalayan arc into a single 
cross-section reveals that these events define a simple planar zone, having an apparent 
northward dip of about 15° (figure 4). This shallow northward dipping zone has been 
interpreted as a part of the master detachment of the Himalaya (Ni and Barazangi 
1984). 

Baranowski etal (1984) suggested that the orientation of slip vectors is 
approximately perpendicular to the Himalayan range. Their method of constraining 
the direction of slips was primarily based on forward modelling of P waves (SH 
waves were used for only a few earthquakes). The orientation of slip vectors is 
approximately perpendicular to the Himalayan range. Thus, the gravitational 
spreading of the overriding Himalaya and Tibet contributes greatly to the 
deformation of the Himalaya. Teleseismic earthquake data indicate that the MBT, 
rather than the MCT, is currently the active feature of the Himalayan arc. It is most 
probable that the great Himalayan events occur along the same detachment surface as 
defined by the moderate-sized thrust-type events (figure 5). 

4.2 P„ and S„ wave propagation beneath the Himalaya 

Variations in P„ and S„ velocities can be used to indicate regional differences in 
the physical properties of the uppermost mantle. Examples of high P„ velocities 
(8-2-8-4km/s) and S„ velocities (4-7km/s) have been reported for a variety of 
continental shield regions (Brune and Dorman 1963; Bath 1966; Huestis et al 1973), 
whereas low P„ velocities (7-7-7-8 km/s) and S„ velocities (44-4-5 km/s) characterize, 
for example, most of the Basin and Range province in which the lithosphere is thin 
and heat flow values are high (Eaton 1980). 

P„ and S„ travel times from earthquakes in the Himalaya to the World Wide 
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Figure 4. A composite cross-section using the MCT and the 4-km elevation contour as a 
reference for projection. Bottom section shows ISC hypocenters. The larger the open circle 
the better the epicenter location. Notice that the majority of earthquakes are located south of 
the MCT In the top section the projected low-angle nodal planes for thrust-type focal 
mechanisms are shown as line segments drawn across the accurately located depths of 
hypocenters. From Ni and Barazangi (1984). 


Standard Seismograph Network (WWSSN) stations at northern Afghanistan, India, 
and to the local network at northern Pakistan, have been used to resolve the 
uppermost mantle P and S wave velocities. For paths along the Himalaya to the 
Tarbela Dam network in northern Pakistan, Menke and Jacob (1976) obtained P„ and 
S„ velocities of 8-5 ±3-5km/s and 5-(X)±2-Okm/s respectively. With WWSSN data, Ni 
and Barazangi (1983) obtained 8-45 ± 0-08 km/s and 4-74 + 0-5 km/s (figure 6), which are 
similar to the corresponding P„ and S„ velocities for the uppermost mantle beneath 
the Indian Shield and Southern Tibet. The crustal and uppermost mantle velocity 
structure of the Himalaya derived from body waves (Ni and Barazangi 1983) is 
consistent with that obtained from surface wave dispersion (Gupta and Narain 1967). 
The high P„ and S„ velocities indicate that the uppermost mantle structure beneath the 
Himalaya is not significantly different from that of India and southern Tibet and is 
consistent with a model in which the Indian continental lithosphere is being 
underthrust beneath the Himalaya and southern Tibet. 







Figure 5. A schematic cross-section of the present geometries of converging plates i 
Himalaya-Tibetan region. Events with normal-type focal rnechanisms located in the sh 
crust and uppermost mantle of Tibet are shown as open circles with crosses. All Hima 
thrust-type events can be interpreted most simply to define a part of the detachment be 
the underthrusting Indian plate and the upper Himalayan crustal blocks. 



Figure 6. Reduced travel-time plots of P„ and S„ data recorded at station NDI, KB 
SHL. At least 50% of the propagation paths are within the Tibetan Plateau. The bottoi 
shows ISC hypocentral depth vs epicentral distance for events used in the above travt 
plots. 



^3 Refraction profiles and wide-angle reflections beneath the Himalaya 

lignificant information on the velocity structure of the Himalaya comes from: (i) an 
inreversed profile that crosses the western Himalaya (Kaila et al 1978; Beloussov et al 
980); and (ii) wide-angle reflection profiles recorded in Nepal (Him and Sapin 1984; 
firn et al 1984a, b). 

Two spreads, one 6 km long and the other 11 km long, located near the Srinagar 
egion of the Kashmir Himalaya recorded shots fired in Pakistan and USSR. These 
iata recorded clear reflection from the Moho. Kaila et al (1978) showed that the 
4 oho is at a depth of about 45 km near Sopur, dips to 54 km in the region of Wular 
.ake and continues to deepen further toward the north. Such geometry indicates a 
Qoderate, 15° to 20° northwesterly dip of the Moho. The dip of the Moho could be an 
iverestimate because the profile is not reversed. 

Wide-angle reflection data in the Tethyan and High Himalaya were collected by 
firn and his colleagues. Seismograms recorded at Arun valley and Kathmandu 
alley from shots in southern Tibet showed two clear phases separated by about 2 s. 
"he later phases were interpreted as a P^P (Moho) reflection from 75km depth 
Lepine et al 1984). The earlier phase was attributed to reflection from the top surface 
>f the underthrusting Indian Plate (Him and Sapin 1984). Although Lepine et al 
1984) suggested that a step in the Moho existed somewhere just north of the High 
limalaya, Molnar (1988) indicated that six seismograms are insufficient to constrain 
he geometry of the Moho. 

Most of the evidence supporting a southward dipping Moho beneath the 
limalaya is from a fan profile constructed with seismograms at distances of 170 to 
00 km recorded in Nepal and southern Tibet (figure 7). The arrival times from the 
rst strong signals recorded by seven stations located north of the High Himalaya 
igure 3 of Him et al 1984a) show a southward dipping reflector from 45 to 50 km 
epth. Him et al (1984a) interpreted this reflector as representing a southward dipping 
doho. Molnar (1988) argued that this reflector can be associated with reflection from 
n interface within the crust. Moreover, he pointed out that reflections in the fan 
Tofile are consistent with a Moho increasing smoothly from 50 to 55 km beneath the 
jreat Himalaya to about 70 km in the southern Tethyan Himalaya. 

.4 Gravity anomalies across the Himalaya 

jravity anomalies in the Himalaya place constraints on the deep structures and the 
ynamics of mountain building. In the Himalaya large Bouguer anomalies become 
icreasingly negative as the elevation increases, their values reach about — 400mGals 
»ver the Great Himalaya (Marussi 1964; Quershy 1969; Quershy and Warsi 1980; 
'''erma and Prasad 1987). This northward decrease in Bouguer anomalies is 
nterpreted to indicate some sort of isostatic compensation in the Himalaya. However, 
arge negative isostatic anomalies (— lOOmGals) over the Ganga Basin and 
-f SOmGals) in the High Himalaya indicate deviation from local isostatic compen- 
ation (e.g. Warsi and Molnar 1977; Lyon-Caen and Molnar 1983). The apparent causes 
)f this deviation from isostatic equilibrium are the focus of recent studies (Lyon- 
2aen and Molnar 1983, 1985). Negative isostatic anomaly over the Ganga Basin can be 
ixplained, at least in part, by the several kilometers of low density upper Cenozoic 
ediments. Positive isostatic anomaly over large parts of the High Himalaya could be 
i result of thickening of a basaltic layer in the lower crust (Quershy and Warsi 1980). 
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Figure 7. A north-south fan-profile across the Greater Himalaya, from Him et al (1984a). 
Reflection times are corrected to a constant offset of200 km. Reflections from about 70 km on 
the north and from about 55 km on the south are quite clear. Northward shallowing 
reflectors to depths of 50 to 45 km are seen on the middle traces. Note strong signals that 
follow these mid-crust reflectors from 65 km (9th and 10th traces) to about 75 km on the 14th 
and 15th traces. Thus, the wide angle reflection data are consistent with a northward dipping 
Moho from 55 km on the southern to 75 km on the northern Greater Himalaya. 


Deficiency and excessive distribution of mass in the Himalaya must be supported by 
stress differences within the lithosphere and mantle. Thus, gravity anomalies in the 
Himalaya can be used to examine how the mountains are supported. In a simple 
elastic flexure model where the Indian lithosphere extends to the Indus-Tsangpo 
suture zone, the overburden is too great to depress the Ganga Basin deeper than 
observed, and the calculated gravity anomalies are much too negative (Lyon-Caen 
and Molnar 1983). Therefore, a simple flexure model cannot account for the gravity 
anomalies across the Himalaya. 

The Indian lithosphere underthrusts beneath the Ganga Basin and Lesser 
Himalaya at a shallower angle than beneath the High Himalaya implying the Moho 
also has different dips. Recognizing this, Lyon-Caen and Molnar (1985) computed 
Bouguer gravity anomalies across the Himalaya using an elastic plate model with 
different flexure rigidity. Their calculations indicate that the load of the mountain is 
still too great to yield gravity anomalies over the Lesser Himalaya and Ganga Basin 
that are comparable with the observed values, A moment must be applied to the end 
of the elastic plate to properly support the Himalayan mountains. This moment can 
be related to either a horizontal force or a vertical force. Lyon-Caen and Molnar 
(1985) ruled out the horizontal force because it must have an excess horizontal 
compressive stress of at least 330 MPa. On the other hand, the vertical force can be 
related to phase change of the lower crustal material in the underthrusted Indian 
lithosphere and is more plausible. 

In our opinion, analysis of gravity data explains the broad features of the Moho 
configuration that are consistent with a smoothly increasing Moho from 38 km 



eneath the Ganga Basin to about 55 km beneath the High Himalaya. Mountain 
uilding involves processes not only confined in the lithosphere but also in the upper 
lantle. The nature of the forces that support the excess weight of the mountain 
;main to be discovered. 


. Pakistan Himalaya 

^formation about the geometry of the underthrusting Indian plate within the 
akistan Himalayan foreland and fold-and-thrust belt has been reported by Lillie 
; al (1987). Based on the interpretation of oil company reflection data they showed 
lat beneath a salt detachment, the Indian crystalline basement dips gently northward 
om about 1° beneath the Salt range to about 4° beneath the northern Potwar 
lateau (figure 8). Their results confirm earlier interpretations that detachment lies 
ithin an Eocambrian evaporite sequence. The extent of this detachment and the 
eometry of the Indian plate farther north beneath the Peshawar basin are not clear. 

Based on an analysis of local seismicity data Seeber and Armbruster (1979) showed 
VO northwesterly trending basement features, the Hazara Lower Seismic Zone 
iLSZ) and the Indus-Kohistan Seismic Zone (IKSZ) (figure 9). They interpreted 
lese features as buried extensions of the northwesterly trending Himalayan thrust 
Lults beyond their surface termination at the western Himalayan Syntaxis. The IKSZ, 
le more active of the two seismic zones, is associated with underthrusting toward the 
ortheast (Seeber and Armbruster 1979). The HLSZ defines a steeply dipping 
asement fault which may merge southward with a shallow dipping surface of 
etachment between the basement and the sedimentary rocks of the Salt Range, 
roving such a relation, however, requires that one knows the deep structures beneath 
lis region. 
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Figure 8. Structural cross-section and Bouguer gravity profile of northern Jhelum Plain, 
Salt Range and central Potwar Plateau, from Lillie et a/ (1987). The gradient in the Bouguer 
gravity in this region is about 1 mgal/km and is similar to that of the Indo-Gangetic Plain. 





Figure 9. Simplified tectonic map of the western Himalayan syntaxis region. MBT, main 
boundary thrust; MMT, main mantle thrust; HTZ, Hazara thrust zone; PT, Panyal thrust; 
SRT, salt range thrust; HLSZ, Hazara lower seismic zone; IKSZ, Indus-Kohistan seismic 
zone. 

A recent study (Ibenbrahim et al 1987) of the three-dimensional P wave velocity 
structure beneath the northern Hazara thrust zone and the eastern Peshawar basin 
from inversion of P wave travel time data produced by the Tarbela Dam array 
indicates that the 6-2km/s velocity surface dips northerly from about 13 km beneath 
the Hazara thrust zone (equivalent of MBT of the central Himalaya) to 17 km beneath 
the northern Peshawar basin. This surface has an apparent dip of about 8° and is 
interpreted as the top surface of the underthrusting Indian Plate. This result 
incorporated with observations from the foreland (e.g. Farah et al 1977; Yeats and 
Hussain 1987; Lille et al 1987) is used to illustrate the geometry of the underthrusting 
Indian lithosphere in the Pakistan Himalaya. Beneath the Salt Range and Potwar 
Plateau the crystalline basement of the Indian Shield dips smoothly northward from 
about 1° to about 4°. The detachment extends farther north beneath the Peshawar 
basin with an increased slope of 5°-10° (figure 10). The relocated hypocenters of 
shallow crustal events cluster at 11-14 km depth aad are associated with the activity 
in the interplate seismic zone beneath the Peshawar basin. A steeper dip of the 
detachment is consistent with the mechanical thrusting model in which shale or other 
incompetent materials rather than evaporite form the zone of decollement (Davis 
et al 1983). Thus, an important implication of the steeper underthrusting of the Indian 
plate beneath the Peshawar basin is that this portion of the detachment appears to 
involve no salt and is capable of producing great earthquakes. 





Figure 10. Structure cross-section (AA' in figure 9) shows that the gently dipping n lan 
plate beneath the Peshawar Basin correlates with the 6-2km/s ?-wave velocity contour. 
From Ibenbrahim et al (1987). 


S. Rates of uplift of the Himalaya 

[t has been suggested by many authors that the uplift “fXdla 

recent geological time, especially during Quaternary and Holocene time (Wadta 1979 
Zeitler et al 1982; Zeitler 1985; Coward et al 1987). A repeated levelling pro i e acros 
the Himalayan frontal zone indicates as much as 6mm/yr of uplift in a 12-year [Kno 
of thHr^t with respect to the Ganga Basin (Chugh 1974; Molnar et al 1977). 
Although the levelling data suggest very rapid uplift, lower rates have been given w en 

Uplift rateTusinrfisSn track methods imply movements of about 1 mm/yr or less 
in Ae t tern Himalaya. However, higher rates of 7 mm/yr have been reported for the 
S Pfirhflt-Haramosh region (Zeitler 1985). The increased uplift around the 
NanS Parbat-Haramosh massif represents the culmination of Nanga Parbat Gneiss^ 
Si miUtaets per year uplift rates in the Higher Himalap have been tnfermd 
from pollen J-ymbedded m ^ tmen s 

that recent uplift is at an accelerating rate. 
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Table 1. Uplift rates of the Himalayan mountains. 


Uplift rates (mm) 

Method of determination 

Source 

0-1 to 4-5 

Fission-track dating 

Zeitler (1985) 

0-6 to 0-8 

Fission-track dating 

Mehta (1980) 

1-0 to 1-4 

River profiles 

Seeber and Gornitz (1983) 

2 

Bengal fan deposition 

Curray and Moore (1971) 

0-1 to 1 

River terraces 

Iwata et al (1984) 

1-2 to 1-5 

Levelling of river terraces 

Iwata (1987) 

6 

Second order levelling 

Chugh (1974) 

5 to 12 

Flexural modelling 

Lyon-Caen and Molnar (1983) 


7. Summary 

The Himalaya is a product of the collision between India and Eurasia which began in 
the Eocene. In the earlier stage of continental collision the development of a suture 
zone between two colliding plates took place. A modern analogue of this stage may be 
the collision between the northern Australian continental margin and the Banda arc 
in the Timor island region. The continued convergence is accommodated along the 
suture zone and in the back-arc region. The Zagros represents an example of this 
stage. Further convergence results in intracrustal megathrust within the leading edge 
of the advancing Indian plate. In the Himalaya this stage is characterized by intense 
uplift of the High Himalaya and the development of the Tibetan Plateau in the back- 
arc region. The tectonics of the modern Himalaya can be best explained by the model 
of Ni and Barazangi (1984). In this model (figure 5), the southern edge of the High 
Himalaya is likened to a curved bulldozer blade that scrapes most sediments from the 
Indian crust, and underthrusting proceeds without fundamental changes in the major 
structure. 

Geological and geophysical constraints presented in this review are still insufficient 
to tightly constrain models outlined in the earlier section. Reflection profiles across 
the entire Himalaya and Tibet are needed to better resolve the deep crustal structures. 
Dense seismic and geodetic networks in the High Himalaya are needed to resolve the 
seismogenic nature of the MCT and the detailed structures. Initiation of 
paleoseimicity studies is necessary to resolve the temporal and spatial patterns of 
major earthquake faulting along the Himalaya. 
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New seismological results on the tectonics of the Garhwal Himalaya 
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Abstract. This paper reports data pertaining to 90 local earthquakes recorded during 
1984-86 using seismographs in arrays of 5-7 stations deployed near the Main Central Thrust 
between Bhagirathi and Alakhananda valleys. The results which are also compared with 162 
earthquakes recorded in 1979-80 provide a local view that refines and complements 
information recorded at distant seismic stations. 
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1. Introduction 

Broadly speaking, the Garhwal Himalaya, drained by the Yamuna, the Bhagirathi 
and the Alakhnanda rivers constitute the middle part of Burrard’s and Gansser’s 
Kumaun Himalaya extending from the Sutlej to the Kali river. The meizoseismal area 
of the 1905 Kangra earthquake (magnitude 8-6) extended into the western part of the 
Garhwal Himalaya (Molnar, this volume). The eartern part of the region marks the 
western end of the 700-km-long seismic gap extending up to the rupture zone of the 
1934 Bihar-Nepal earthquake (magnitude 8-4). The region has not experienced a 
major earthquake (magnitude > 8 0) in recorded history and unless it happens to be 
an exceptional segment of the Himalaya collision zone distinguished only by strain 
releases through small and intermediate magnitude earthquakes, the ambient stresses 
here must be fairly high. 

In this study, data for 90 local earthquakes recorded during 1984-86, using 
portable seismographs in arrays of 5 to 7 stations deployed in the vicinity of the Main 
Central Thrust (MCT) between Bhagirathi and Alakhnanda valleys, are reported and 
considered along with the previously published results for 162 earthquakes recorded 
in 1979-80 in the vicinity of the MCT between the Yamuna and Bhagirathi valleys. 
The epicentres of these 252 earthquakes define a 140-km-long seismic belt stretching 
from northwest of the Yamuna valley southeastward upto the Alakhnanda valley. 
The width of the belt varies between 30 and 50 km. The estimated focal depths of 
225 earthquakes were less than 13 km, while those of the remaining 27 were between 
13 and 23 km. The first motion data for 50 earthquakes provided two composite focal 
mechanism solutions. Data for 29 earthquakes defined a strike-slip solution. Twenty 
of these earthquakes occurred NW of the Bhagirathi river. Data for the other 21 
earthquakes defined a thrust solution. Nineteen of these earthquakes occurred SE of 
the Bhagirathi river. The maximum compressive stress direction is very nearly 
identical in both cases, being oriented NE-SW. The thrust solution yields a modest 
plunge of 15° to the NE for the axes of maximum compressive stress. 
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or reverse faulting appears to exist in the overriding plate from the ground surface 
down to its lower boundary where it is in contact with the under-thrusting Indian 
plate. The strike-slip earthquakes may be explained by postulating a modification of 
this stress-state through tensile stresses due to a flexure of the overriding plate over a 
ridge-like feature on the top surface of the subducting plate. The ridge would be a 
northeastward salient of the Pre-Cambrian Aravalli ranges of Rajasthan whose 
outcrops can be seen as far north as New Delhi and which probably continue 
northeastward under the Sindhu-Ganga alluvium and even the Lesser and High 
Himalaya. The crest of the ridge may be approximately along 78°E. The tensile 
stresses generated on this account would attenuate with distance perpendicular to the 
ridge crest. The ridge which is transverse to the Himalaya perhaps decouples the 
rupture zone of the 1905 Kangra earthquake from the seismic gap that lies east of it, 
thereby facilitating the strain release on these segments at their separate rates. 

The seismic belt associated with the Himalayan collision boundary has so far been 
delineated mainly on the basis of teleseismic and regional data. Here we present 
information about earthquake locations and source mechanisms for a part of the 
Garhwal Himalaya (figure 1) segment of the seismic belt as gleaned from seismograms 
recorded locally, using portable seismographs deployed in arrays of small areal extent. 
Combined with data and results from Gaur et al (1985) (hereafter referred to as 
paper I) this provides a unique local view that refines and complements the 
information recorded at distant seismic stations. 


1.1 Nomenclature 

The administrative districts of Pauri, Chamoli, Uttarkashi, Tehri and Dehradun 
constitute the region of the Garhwal Himalaya. This region forms a somewhat bigger 
middle third of Burrard’s (see Wadia 151, p. 9 and Gansser 1964, p. 8) “Kumaon 
Himalaya” which extended from Sutlej to Kali river. Alternatively, the region forms 
the western, bigger part of Valdiya’s (1980, p. 1) “Kumaun Himalaya” which run from 
the state boundary between Himachal Pradesh and Uttar Pradesh to the 
international boundary between India and Nepal along the Kali river. The Yamuna, 
the Bhagirathi and the Alakhnanda are the main rivers of the Garhwal Himalaya. 

We shall have occasion to use the terms Lesser and High Himalaya. Originally 
(Wadia 1951, p. 10) these terms denoted physiographic provinces running along the 
entire length of the Himalaya. “The High (or Great) Himalaya constitute the 
innermost line of high ranges rising above the limit of perpetual snow. Their average 
height extends to 6-7 km. The Lesser Himalaya are a series of ranges closely related to 
the former but of lower elevation, seldom rising much above 4 to 5 km” (Wadia 1951). 
Gansser (1964) used terms Lower and Higher Himalaya with the connotation that the 
boundary between the Lesser and Higher Himalaya is marked by the Main Central 
Thrust (MCT), the other boundary of the Lesser Himalaya being the Main Boundary 
Thrust (MBT). According to Valdiya (1980, p. 10), the highest point in Garhwal 
Himalaya is Nag Tibba (elev. 3022 m), the general elevations being in the range of 
1500 to 2500 m. 

We shall call the seismic belt in the Garhwal Himalaya as GLHSB, short for 
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Figure 1. Showing the location of the area of interest in the Himalaya. Abbreviations for 
river names are: AR, Alakhnanda river; BH, Bhilangana river; BR, Bhagirathi river; GR, 
Ganga river; SR, Sutlej river and YR, Yamuna river. Abbreviations for place names are: BD, 
Badrinath; DD, Dehradun and ND, New Delhi. The thrusts shown are the Main Boundary 
Thrust (MBT) and the Main Central Thrust (MCT). 


hwal Lesser Himalaya Seismic Belt. This appellation appears to be strictly true 
i respect to the earlier definition of the Lesser Himalaya. It takes only minor 
rties with Gansser’s Lower Himalaya concept because at places in Garhwal the 
Tiic belt appears to lie some distance NE of the MCT. 


Relevance of the study 

r large earthquakes (magnitudes 8-7, 8-6, 8-4 and 8*7, Richter 1958) have occurred 
iifferent parts of the Himalaya in the last nine decades. But the last of them 
irred 37 years ago. The meizoseismal area of the 1905 Kangra earthquake 
gnitude 8-6) extended well into Garhwal Himalaya from the west. The portion of 
hwal Himalaya immediately east of this is a part of the nearly 700-km-long 
Tiic gap recognized by Seeber and Armbruster (1981), Khattri and Tyagi (1983) 
Khattri (1987). The region is not known to have experienced a magnitude 8 or 
Lter earthquake in recorded history and unless it happens to be an exceptional 
nent of the Himalaya collision zone in which most of the stress is released through 




Further, as noted above, the Garhwal Himalaya are drained by three mighty 
perennial rivers with numerous tributaries. The region thus possesses a vast potential 
for irrigation and generation of hydroelectric power. A large number of multipurpose 
river valley schemes have been conceived for the region and many of them are under 
execution or in advanced stages of design. The 260 m high earth and rockfill Tehri 
Dam on the Bhagirathi river is in the early stages of execution. Detailed seismicity 
information would be vital for planning seismic risk mitigation measures for these 
projects. Scientifically too, detailed information about earthquake processes in the 
region would be useful for evolving earthquake prediction and seismotectonic models. 


1.3 Local recording 

Sporadic recording of local earthquakes using one to three seismograph stations has 
been carried out in the Bhagirathi valley between 1974 and 1977 (Agarwal and Kumar 
1982). Four permanent photographically recording seismographs have also been 
present for some time and are still in operation in the region (Srivastava 1986). But the 
instruments have low magnification and the array configuration is unsuitable for 
detecting and locating many local earthquakes. Results of recording local earthquakes 
with a five-station array (largest dimension 45 km) during 1979-80 were reported in 
paper I. The array stations were located in the Yamuna and the Bhagirathi valleys in 
the vicinity of the MCT. We report here analyses of local earthquake data recorded 
with 5 to 7 element seismograph arrays during 1984-86 in the region immediately SE 
of that considered in paper I. This region lies between Bhagirathi and Alakhnanda 
valleys in the vicinity of MCT. 

1.4 Review 

A cardinal fact about the Himalayan seismic belt is that its overall definition so far 
rests on regional and teleseismic observations of intermediate magnitude earthquakes. 
The information has been supplemented in recent years by data from a few local 
networks (Seeber et al 1981; Gaur et al 1985; Srivastava 1986; Khattri 1987). Hence 
most of the articles about seismicity of Himalaya refer to the four large and numerous 
intermediate magnitude earthquakes. We recall here the available models of 
Himalaya against which we compare our results and conclusions. 

The Himalaya is considered to be the result of continent-continent collision of the 
Indian and Eurasian plates (e.g. LeFort 1975; Valdiya 1980; Seeber et al 1981; Khattri 
1987). The two continents have sutured along the Indus Suture Thrust. Portions of the 
northern edge of the Indian plate have been thrusted southwards as a result of 
continued convergence between India and Eurasia. Two major intraplate thrusts are 
the MCT and MBT. According to one model, the so-called evolutionary model (e.g. 
LeFort 1975), the MCT and MBT are similar in nature and along these crustal slices 
of the Indian plate have been thrusted southwards. In the alternative steady-state 
model (e g. Seeber et al 1981; Seeber and Armbruster 1981; Ni and Barazangi 1984) the 
MBT and MCT are different in character. The MCT is a thrust fault that involves the 
upper crust. The region south of the MCT constitutes a sedimentary wedge consisting 
of the Ganga foredeep, the other (Siwalik) Himalaya, and the Lesser Himalaya. This 



wedge IS thought to be underlain by a quasi-horizontal detachment surface coinciding 
with the upper surface of the subducting Indian lithosphere. Earthquakes of 
intermediate magnitudes occur near the northern margin of the detachment whereas 
the great earthquakes further to the south are caused by slip on this detachment. 

The faultplane solutions and focal depths determined by Baranowski et al (1984) 
confirm the above model, though they, as well as Lyon-Caen and Molnar (1983), 
favour a model in which the upper portion of the Indian slab under the Greater 
Himalaya has been sliced and thrusted southward. 

On a more local scale, for the Garhwal Himalaya specially, Agarwal and Kumar 
(1982) reported their tripartite observations of seismicity in the Bhagirathi valley. 
They noted sporadic activity along the MCT. 

A belt of local earthquake epicentres was defined over a distance of about 70 km 
from the analyses reported in paper I. The best data pertained to the region between 
the Yamuna and the Bhagirathi valley which was enclosed within the recording array. 
In this region the seismic belt was observed to lie distinctly south of the MCT which 
was thus considered to be seismically inactive in the region. A sinistral offset along on 
E-W line was observed in the trend of the belt. A composite focal mechanism solution 
indicated strike-slip faulting. The strikes of the two nodal planes were N-S and E-W. 
The latter plane was inferred to be the fault plane because the sinistral relative motion 
observed across it coincided with the sinistral E-W offset in the seismic belt. The focal 
depth of the majority of the earthquakes was less than 10 km and a case was made 
that the seismicity may be confined only to the top 20 km of the crust. 

The new results presented in this study enable the belt of local earthquakes to be 
defined over a further distance of 70 km between the Bhagirathi and Alakhnanda 
valleys. The estimated focal depths of earthquakes are again predominantly shallow. 
A composite ‘thrust’ solution has been observed for earthquakes SE of the Bhagirathi 
valley. The seismicity in the Garhwal Himalaya appears to be confined primarily to 
the overriding slice of the plate which may be composed of the sedimentary wedge and 
the Tethyan slab according to Seeber and Armbruster (1981), of the upper portion of 
the Indian plate itself (Baranowski et al 1984). 


2. Seismicity 

The sites in the Bhagirathi and Alakhnanda valley occupied with portable 
seismographs during parts of 1984-86 are shown in figure 2. The details of the station 
networks operated are given in table 1. Earthquakes occurring locally as well as at 
regional and teleseismic distances were recorded. But we present data for local 


Table 1. Details about seismograph arrays operated during 1984-86. 

List of stations in the array 

Period of operation 

AKM, CHA, LMG 

AKM, CHA, GHU, LMG, DHO 

AKM, CHA, GHU, DHO, LAT 

AKM, CHA, GHU, DDA, DAB, TIL, UKH 

November-December 1984 
January-March 1985 
April-June 1985 
October 1985-June 1986 


Note: The station codes match with those shown in figure 2 





Layer 1 
Layer 2 


17 

Infinite 


5- 2 

6 - 0 


2- 97 

3- 43 


earthquakes only because by doing so at least the estimates of epicentral locations can 
be relied upon for earthquakes o.ccurring within or close to a seismograph station 
array. 

2.1 Data analysed 

The data for only those earthquakes were analysed which provided a minimum of 
five-wave arrival times in the following combinations. 

(i) P readings at 5 stations; 

(ii) P readings at 4 stations and S reading at one of the 4 stations; 

(iii) P readings at 3 stations and S readings at two of these stations. 

A restriction was also imposed that the difference in arrival times of P and S waves 
should be less than 10 s. The aim was to improve the chances that the first arriving P 
wave would be a direct wave. A simple application of the head wave theory (Dobrin 
1976) shows that in the context of the wave speed model of table 2, this would be the 
case for earthquakes with focal depths less than 13 km. Still, in every case, an attempt 
was made to use as much data for an earthquake as possible and available within 
these restrictions. 

2.2 Hypocentral parameter estimation procedure 

The same computer program was employed for this purpose as discussed and used in 
paper I. 

2.3 Crustal model 

Estimation of the four hypocentral parameters (viz. epicentral latitude and longitude, 
hypocentral depth and origin time) from arrival times of P and S waves requires 
information about the distribution of P and S wave speeds between the earthquake 
focus and the stations. We are unaware that any study elucidating this information 
has been carried out in the region of interest, namely the region along the GLHSB 
between the Bhagirathi and the Alakhnanda rivers. We adopted for the purpose a 
uniform half-space model with a P wave speed of 5-2 km/s and P to 5 wave speed ratio 
of 1-75, close to 1*73 for a Poisson solid. This P wave speed is the same as estimated by 
Chander et al (1986) for the region in the vicinity of the MCT between Yamuna and 
Bhagirathi valleys immediately to the NW of the region of interest here. We offer the 
following remarks in defence of this simple wave speed model. Many investigators (e.g. 
Kaila et al 1968; Tandon and Dubey 1973; Verma 1974) have attempted to estimate 





Lve also oecn usea lo esumaie seismic wave speeas in me iirsi auu sccoiiu oiuaiai 
yers. The above mentioned result (5-2 km/s) for the first crustal layer applies to the 
gion between Yamuna and Bhagirathi valleys near the MCT, while the result 
0 km/s) for the second crustal layer (Kumar et al 1987) applies strictly to the region 
(mediately NW of that. Both these values have been adopted for the region of the 
LHSB southeast of the Bhagirathi valley as being the values estimated for the region 
)sest to the region of interest. The depth to the interface between the two layers has 
>t been estimated in these studies. We assume that the interface in question is the 
rust surface between the colliding lithospheric plates in the Lesser and High 
imalaya. The surface may be delineated approximately by the foci of intermediate 
ignitude thrust earthquakes (Ni and Barazangi 1984) in the Himalaya seismic belt, 
le depth to the interface beneath the MCT is estimated by us to be 17 km from Ni 
d Barazangi’s figure 11 (middle part). The interface dip of 15° has been estimated by 
and Barazangi but may be ignored when only the order of magnitude of the travel 
ries of head waves travelling predominantly parallel to the strike of the interface is 
ing considered. Since only eight out of 92 earthquakes had estimated focal depths 
low 13 km in the context of the model of table 2, these were the only earthquakes for 
lich there was a need to consider at least a two-layered half space model. This 
itifies the use of the simple half-space model on grounds of convenience and its 
atively minor influence on the results obtained. 


I- Seismicity results 

e stringent conditions on phase data discussed in the foregoing were met for only 41 
rthquakes recorded during the 1984-85 season and 49 events during the 1985-86 
ison. The data for these 90 events are considered along with 162 events considered 
paper I which met the same requirements. However, we recall that the hypocentral 
rameters of the latter events had been estimated using a P wave speed of 5 km/s. 
The 252 epicentres are plotted in figure 2. The seismic belt delineated by these local 
5nts is compared with that defined by regionally and teleseismically located events 
1684 to 1985 in figure 3. Two depth sections along the perpendicular to the trend 
the seismic belt (local earthquakes only) are shown in figures 4 and 5. The coda 
Ignitude of these events is in the range of 2 to 4. 


» Belt of seismicity 

c 252 epicentres of the local earthquakes recorded during 1979-80 and 1984-86 
fine a definite seismic belt extending over a distance of about 140 km from NW of 
imuna river to Alakhnanda river. The width of the belt appears to vary between 30 
d 50 km. Around 78-5°E longitude the belt appears to have a 50 km long, N-S 
nding offshoot to the south. The concentration of epicentres between Yamuna and 
lagirathi valleys clearly defines a sinistral E-W oriented offset in the belt (figure 2). 


I Focal depth 

frequency/km table (table 3) and a histogram (figure 6) were prepared for the 



Figure 2. A combined plot epicentres of local earthquakes recorded in the Garhwal Himalaya during 1984-86 (90 events) and 1979-80 (1 
represent epicentres and squares stations. The stations operated during 1984-86 and are shown here as AKM, LAM, DHO, LAT, BN A, 
UKH, DHO. The schedule indicating the stations operating at a given time during 1984-86 is given in table 1. 




The epicentre of the 1905 Kangra earthquake (M = 8-6) is also included. The location of this earthquake has been taken from Molnar 
Delhi-Hardwar Ridge. 
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Table 3. Distribution of estimated focal depths 
in 1 km intervals. 


Depth interval Number of Cumulative 

(km) earthquakes Number 

in the interval 


-1-0 

1 

1 

0-1 

20 

21 

1-2 

23 

44 

2-3 

25 

69 

3-4 

19 

88 

4-5 

20 

108 

5-6 

22 

130 

6-7 

25 

155 

7-8 

11 

166 

8-9 

16 

182 

9-10 

11 

193 

10-11 

11 

204 

11-12 

12 

216 


13-14 

9 

234 

14-15 

3 

237 

15-16 

4 

241 

16-17 

2 

243 

17-18 

3 

246 

18-19 

2 

248 

19-20 

0 

248 

20-21 

2 

250 

21-22 

1 

251 

22-23 

1 

252 



Figure 6. Histogram showing the numbe 
estimated 252 earthquakes. 








estimated focal depths of the 252 earthquakes. The estimated focal depth did not 
exceed 23 km in any case. It was less than 10 km for 193 earthquakes and between 10 
and 13 km for another 32 earthquakes. In all 27 (or about 10%) earthquakes were 
assigned estimated focal depths in the range of 13 to 23 km. The preponderance of 
very shallow earthquakes among the local earthquakes considered is manifest. 

2.7 Activity with time 

The frequency of earthquake occurrence in time is not uniform along the trend of the 
GLHSB. The region with the highest level of activity in this regard is the part of the 
belt between the Yamuna and Bhagirathi valleys. The next most active part of this 
stretch of the belt is considered on available data to be between the Mandakini and 
Alakhnanda rivers. This stretch of the belt appears to be most active, however, when 
the number of earthquakes located from regional and teleseismic observation in the 
period 1684-1985 is considered. 


3. Focal mechanisms 

Focal machanisms of earthquakes have provided crucial evidence for sea floor 
spreading and plate tectonic hypotheses. Fitch (1970) presented the first focal 
mechanism solutions based on teleseismic and regional data for earthquakes in the 
Himalayan region. The results for the four earthquakes which occurred just east of 
Garhwal and were considered by Fitch showed support for the view that the Indian 
lithospheric plate was underthrusting the Himalaya. Subsequently many other 
investigators have presented such solutions for the Himalaya, but none of them 
referred particularly to the Garhwal Himalaya. Among the Himalaya earthquakes 
with available focal mechanism solutions, six occurring in the India-Nepal Border 
region immediately east of the Garhwal Himalaya showed thrust-faulting and had 
reliably determined focal depths in the range of 12 to 18 km (Ni and Barazangi 1984; 
Baranowski et al 1984). 

The seismograms recorded during 1984-86 were also examined for the sense of 
ground motion in the first P waves. Reliable readings were obtained at one or more 
stations for 28 earthquakes. The data were consistent with two distinct composite 
focal mechanism solutions as follows. 

Data from 20 earthquakes were consistent enough to yield the composite solution 
of figure 7. Five of the earthquakes occurred at focal depths of 15 to 20 km. We take 
the NE dipping nodal plane with reverse fault motion as the fault plane. This is 
consistent with the view held by earlier workers (e.g. Fitch 1970; Chandra 1978) for 
Himalayan earthquake to take the northerly dipping plane as the fault plane. Ni and 
Barazangi (1984) and Baranowski et al (1984) reported thrust faulting in the 12 to 
18 km range. The NE dipping plane of figure 6 has a dip of 60° and it should strictly be 
termed a reverse fault. Since it refers mainly to shallower earthquakes, it is in 
conformity with the geological observations and opinion (e.g. Valdiya 1986) that the 
thrust faults in the Himalaya are listric with steeper dips at shallow depths and gentler 
dips at greater depths. Accordingly, we also call our reverse fault solution as a ‘thrust’ 
solution. Numerical data used for this solution is given in table 4. 

Data from eight earthquakes were found to be consistent with the same strike-slip 



N 



Figure 7. Composite focal mechanism solution in which only data from 1984-86 were 
involved. Equal area projection of the upper half of the focal sphere is shown. P and T mark 
the positions of the pressure and tension axes. The north-easterly dipping nodal plane is 
identified as the fault plane. 


Table 4. Numerical data for the thrust solution of 
figure 6. 



Strike 

Dip 

Nodal plane 1 

N48°W-S48°E 

60° along N42°E 

Nodal plane 2 

N48°W-S48°E 

30° along S42°W 


Trend 

Plunge 

Pressure axis 

N42°E-S42°W 

15° along N42°E 

Tension axis 

N42°E-S42°W 

75° along S42°W 

b axis 

N48“W-S48°E 

Nil 


solution as obtained in paper I from 1979-80 recordings. Seven of the eight 
earthquakes had their epicentres in the SE part of the region where earthquakes used 
in solution of paper I had occurred. One of these seven earthquakes had an estimated 
focal depth of 17-7 km while the rest had focal depths shallover than 15 km. The eight 
earthquakes occurred well inside the region of earthquakes contributing to the above 
‘thrust’ solution. The data on first motion for these earthquakes are shown in figure 8 
plotted on the data and strike-slip solution of paper I. The new data 1984-86 provide 






Figure 8. Composite focal mechanism solution in which 1984-86 data are combined with 
the 1979-80 data. Equal area projection of the upper half of the focal sphere is shown. P and 
T mark the positions of the pressure and tension axes. 


no reason to change the earlier conclusion regarding E-W nodal plane with sinistral 
relative motion as the fault plane, because it matches with a sinistral E-W offset in the 
epicentral belt. 

The epicentres of all earthquakes whose data were used in obtaining the composite 
solutions are shown in figure 9 where T and S mark the epicentral positions according 
to whether the corresponding first motion data appear in figures 6 or 7 respectively. 

A remarkable feature of the composite solutions is that the direction of compressive 
stress is very nearly the same NE-SW in both cases (figures 6-8). 


4. Discussion 

It is evident from a review of literature and the foregoing analyses that investigations 
so far made, of Himalayan earthquakes based on magnitude, may be divided into 
three categories. Investigations of large (magnitude > 8) earthquakes are mainly 
based on a visual inspection of the effects, in the four known cases of the last nine 
decades. Only for the 1934 Bihar-Nepal and 1950 Assam earthquakes are some 
instrumental investigations available (apart from estimation of hypocentral 
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Figure 9. Epicentres of the earthquakes whose first motion data were used in the composite 
plots. T and S mark the epicentre positions according as the data for the corresponding 
earthquake were utilized in the thrust solution of figure 6 or strike-slip of figure 7. The large 
arrows indicate continuity in the maximum compressive stress direction across the 140 km 
length of the belt. 


ameters and magnitudes). The intermediate magnitude earthquakes have been 
orded at regional and teleseismic stations and have yielded data on hypocentral 
•ameters, magnitudes, focal mechanisms, space-time seismicity studies, etc. In the 
t few decades local networks have started to yield information on small magnitude 
thquakes in the Himalaya. The networks operated by our group in Garhwal 
nalaya during 1979-80 and 1984-86 have yielded data on only one earthquake for 
ich NOAA also has provided hypocentral data. In the modelling of seismicity of 
nalaya attention has only been paid (e.g. Seeber and Armbruster 1981; Khattri and 
agi 1983; Baranowski et al 1984; Ni and Barazangi 1984, 1987) to large and 
irmediate magnitude earthquakes. 

rhe belt of local earthquakes epicentres investigated by us in the Garhwal 
nalaya appears to agree with that deduced from regional and teleseismic data 
ure 3). The epicentre of the earthquake of 28 December 1979 was located using 
al data and by NOAA using regional and teleseismic data. The estimates of the 
) epicentres differed by 23 km of map distance. This implies that the epicentral 
>rdinates and focal depths of intermediate magnitude earthquakes are subject to 
certainties of this order. We thus take the view that the agreement between small 
1 intermediate magnitude earthquake epicentre belts may be regarded as tentative 
I subject to further verification. 

9epthwise, the small magnitude earthquakes in the Garhwal Himalaya, according 
our estimates, have shallow foci from ground surface down to a depth of about 
cm in a verv laree number of cases. On the other hand, the accurate focal depth 



depth wise region oi low seismicity in the 12 to 10 km range has been speculated upon 
by Chander et al (1985) but bears further investigation. Also, we still hold the view put 
forward in paper I on the basis of numerical experiments with observed and synthetic 
data and multilayered wave speed models that the focal depths in the range of 30 to 
70 km reported by ISC and NEIS for Himalayan earthquakes are probably unreliable. 
This is borne out by the data presented in many figures by Ni and Barazangi (1984) 
and Baranowski et al (1984). 

The relation of the seismic belt to the MCT is of interest. Seeber and Armbruster 
(1981) regarded the seismic belt to be down dip of the MCT. Baranowski et al 
(1984) and Ni and Barazangi (1984) argued that the belt lies between MBT and MCT 
mostly close to the latter. The same inference was also drawn in paper I. However, the 
surface location of the MCT is a matter of debate among investigators. The trace of 
the MCT shown in the figures of the present article is taken from the tectonic map of 
Gansser (1964) because he was one of the propounders of the MCT concept and his 
map is widely available and frequently referred. The results of seismicity obtained in 
this article on being combined with those of paper I show that over a distance of about 
140 km, the seismic belt follows a straighter course than the surface trace of the MCT. 
In fact the seismic belt and the MCT criss-cross each other in this part of the Garhwal 
Himalaya. 

The status of the MCT and the MBT in regard to the current seismotectonics of the 
Himalaya is a matter of debate among the earth scientists. For example, LeFort (1975) 
considered MCT and MBT as successive tectonic boundary thrusts. MBT is the 
current boundary of the continental convergence zone while the MCT is a less active 
or dormant feature. Seeber and Armbruster (1981) argued that MBT and MCT are 
contemporaneous and that MCT is currently active. Ni and Barazangi (1984) state 
that it is the MBT and nearby surface and blind thrusts rather than the MCT that are 
presently the most active structures in the Himalaya. Chander (1987, in press) found 
that the coseismic ground elevation changes observed during the 1905 Kangra 
earthquake can be explained only by assuming that slip did not take place on the 
MBT. Lastly Krishnaswamy et al (1970) observed that coseismic slip did not occur on 
MBT or any other thrust in the Himalaya. We thus take the view that the current 
seismotectonics of the Himalaya is associated with buried (or blind) thrust(s) 
corresponding to the upper surface of the subducting Indian plate. 

The main cause of the stresses responsible for the small magnitude earthquakes 
recorded by us is probably the same as that responsible for relative displacements 
between the overriding and the subducting plates in the Himalaya. They are likely to 
be conducive to thrust earthquakes. Locally the elastic (plastic) response of the 
overriding plate to topographic relief in the underthrusting slab may lead to 
perturbations of this stress regime. Chandra (1978) put forward the agrument that the 
normal faulting during the 1975 Kinnaur earthquake, approximately 150 km north of 
MCT near 78°E, was in response to stresses induced by a north-eastward extending 
salient of the northeast trending Precambrian Aravalli ranges of Rajasthan under the 
Sindhu-Ganga (Indo-Gangetic) alluvium and the Lesser and High Himalaya. We 
recapitulate in figure 10 how the principal stresses for a thrust fault or earthquake can 
be modified with a superimposed tensile stress to produce principal stresses suitable 
for strike slip faulting and earthquakes. We thus argue (figure 11) that in the region of 
GLHSB between the Yamuna and Bhagirathi valleys, the flexing of the overriding 
rock material in response to the presence of the same Aravalli salient on the under- 




THRUST FAULT STRESS TENSILE STRESS STRIKE SUP FAULT 

REGIME DUE TO FLEXURE STRESS REGIME 

Figure 10. Schematic drawing to show the alteration of a thrust-fault stress regime to a 
strike-slip fault stress regime due to action of stresses arising from flexure in the lithospheric 
plate. 



Figure 11. A cartoon to show the origin of flexure in the overriding plate due to 
topographic relief in the subducting plate. DHR is the Delhi-Hardwar Ridge, a postulated 
buried northeastward salient of the Pre-Cambrian Aravalii ranges of Rajasthan under the 
Sindhu Ganga (Indo-Gangetic) Plains (SGP) and Outer (OH), Lesser (LH), Bhagirathi (BR), 
Ganga (GR) and Yamuna (YR). 

thrusting Indian slab gives rise to tensile stresses that modify the configuration of 
stresses that would have been conducive to thrust earthquakes in the absence of the 
salient into a stress configuration conducive to strike-slip earthquakes. South east of 
the Bhagirathi valley, the influence of this superimposed tensile stress regime wanes 
sufficiently with the result ‘thrust’ faulting is observed even at shallow levels under the 
Lesser Himalaya. We however note that in a setting as complex as the Himalaya other 
processes can give rise to strike-slip faulting. 


5. Conclusions 

(i) Our observations with seismograph arrays of small areal extent in the Garhwal 
Himalaya have so far shed light only on the small magnitude local earthquakes. 

(ii) The belt-of local earthquakes has been defined over a distance of 140 km parallel 
to the general strike of the mountains in the Garhwal Himalaya. 



1 ne Dell oi lOcai earinquaKes appears lo comcioe wim me oeii oi miermeuiaie 
magnitude earthquakes. However, as uncertainties in the locations of the latter events 
are greater, this agreement can only be regarded as being tentative. 

(iv) The seismic belt follows a straighter course than the outcrop of the MCT in the 
Garhwal Himalaya and the two criss—cross each other. 

(v) The vast majority of earthquakes have focal depths of less than 13 km. 

(vi) The first motion data of this paper on being combined with that of paper I 
indicate strike-slip and thrust mechanisms in operation along the 140 km length of the 
GLHSB, investigated. Both solutions indicated a NE-SW oriented compressive 
maximum principal stress. 

(vii) The composite solution indicating thrust-faulting establishes that the overall 
stress regime in the region is conducive to thrust or reverse fault earthquakes. 

(viii) The strike slip earthquakes in the region occupied by the valleys of Yamuna to 
Bhagirathi, may be the result of an aberration of this stress regime due to tensile 
stresses induced in the overriding slab, in response to topographic relief on the 
subducting slab surface. 

(ix) The seismic risk in eastern Garhwal Himalaya is high because it lies in the area of 
a seismic gap along the Himalayan collision zone. 
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1 intermediate-term prediction of strong earthquakes in the 
malayan arc region using pattern recognition algorithm M8 
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Abstract. Seismicity of the Himalayan arc lying within the limits shown in figure 1 and 
covering the period 1964 to 1987 was scanned using M8 algorithm with a view to identifying 
the times of increased probabilities (TIPs) of the occurrence of earthquakes of magnitude 
greater than or equal to 7-0, during the period 1970 to 1987. In this period, TIPs occupy 
18% of the space time considered. One of these precedes the only earthquake in this magnitude 
range which occurred during the period. Two numerical parameters used in the algorithm, 
namely the magnitude thresholds, had to be altered for the present study owing to incomplete 
data. Further monitoring of TIPs is however warranted, both for testing the predictive 
capability of this algorithm in the Himalayan region and for creating a base for the search 
of short-term precursors. 

Keywords. Himalaya; intermediate term prediction; M8 algorithm; pattern recognition 
algorithm; times of increased probabilities; seismicity. 


Introduction 

any scheme of pattern recognition, distinctive patterns of attributes describing the 
its of events (objects) of a particular class are first established in the course of a 
rning process. The ensemble of objects which are thus unequivocally grouped into 
3 or more classes and the attributes describing them, constitute the learning material 
m which recognition rules are abstracted. These rules are then applied to recognition 
lew objects, that is assigning them to the appropriate class, which is the ultimate 
ective. For predicting earthquakes, the objective is to recognize instants of times 
ich presage the approach of a strong earthquake. This can be attempted by 
ablishing premonitory patterns based on the analysis of past data and designing 
appropriate pattern recognition algorithm. 

earthquake sequences constitute so far the longest and the least incomplete data 
relevant to earthquake prediction in many parts of the world. This is also largely 
B for the Himalayan region. Here we analyse the earthquake sequences in the 
nalayan region using the M8 algorithm designed by Keilis-Borok and Kosobokov 
84) to test the possibility of diagnosing the approach of a strong earthquake. The 
monitory patterns expressed by this algorithm appear to be rather similar in widely 
erent regions of the world (Gabrielov et al 1986). 

t is of course recognized, at the very outset, that the symptoms of an approaching 
thquake can be quite different from one case to another. This may happen because 
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the distinctive seismotectonic environment of different seismic regions or due to 
lomogeneous distribution of rock materials and stresses within the same region, 
en earthquakes occurring in the same area may show different premonitory 
;ociations if the stress regime and rock properties preparatory to fracture have been 
iatly altered by earlier or neighbouring earthquakes. The nature of earthquake 
[uences is therefore so diverse (Gabrielov et al 1987), as to defy classification and 
strate attempts to abstract from them a consistent set of premonitory signals. In 
ier to circumvent this problem, Allen et al (1984) and Gabrielov et al (1986) devised 
cheme of integral representation of earthquake sequences which smoothens the 
ersity of seismicity patterns without, however, losing any premonitory signals that 
y be contained in the data. This is realized in the following way. 

First, the earthquake flow within a certain large territory is carefully scanned, 
^eral integral functions (attributes) describing the nature (traits) of the earthquake 
[uence in the region are then evaluated within a sliding time window, each assigned 
the instant of time that terminates the window. In all, 14 functions describing the 
lowing traits can be formed: level of seismicity, seismic quiescence, variation of 
smicity, spatial concentration, clustering in time and space, contrast of seismicity 
adjacent areas and long term interaction. 

figure 2 shows a schematic representation of traits and the manner of their 
iluation along the time axis. Appropriate pattern recognition algorithms can be 
^eloped using some or all of these traits. Diagnosis of the time of increased 
)bability of occurrence of a strong earthquake (TIP) is made by applying the 
3ropriate recognition rules. After a TIP has been diagnosed, a strong earthquake 
nagnitude greater than or equal to a threshold magnitude Mq is expected to occur 
hin a few hundred kilometers over the next few years. This is called an 
ermediate-term prediction. A magnitude threshold Mq is determined for each 
ion on the basis of its past history of strong earthquakes, 
fo ensure uniform diagnostics of TIPs for different regions, the functions describing 
ious traits are normalized, as described below. The sequence of main shocks is 
t isolated, the number of aftershocks being considered as one of the characteristic 
ctions of the corresponding main shock. In order to smoothen the earthquake 
uence, their traits are first defined within a wide time window and epicentral area, 
xt, groups are formed corresponding to each of the functions, which comprise a 
of these functions evaluated for different arguments. For example, the group of 
function describing the intensity of earthquake flow consists of the values of the 
ction evaluated for different magnitude ranges. Each function is then quantized 
two or three discrete intervals, and various combinations of these are used to 
gnose a TIP. 

fwo related algorithms for such diagnosis of TIPs are available; one called CN, 
igned by Allen et al (1984) and the other, M8 designed by Keilis-Borok and 
sobokov (1984). Their application to different regions of the world is described by 
brielov et al (1986). Algorithm M8, used in the present study, is described below. 
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Figure 2. Schematic representation of the sequence of earthquake occurrence and evaluation 
of intrinsic traits for PR methods. Vertical lines on the time axis show the occurrence of 
earthquakes (the length being proportional to the magnitude). Various traits are evaluated 
along sliding time window at At steps. Broken vertical lines show this at t and t + At time 
instants. 


of sources and clustering of sources, measured by the maximal number of aftershocks 
of the ihain shocks. The territory under study is scanned by a network of overlapping 
squares. These squares are centred at appropriate points which show comparatively 
greater seismic activity with respect to their surroundings. The side length (L) of each 
square in degrees of the great circle is determined by the relation 

L(Mo) = exp(Mo — 5-6) + 1. 

The above functions are evaluated from the data of earthquake flow in these squares. 
Each of the first three functions is determined twice, for two different magnitude 
thresholds M i and M 2 . The set of values of a function evaluated for different arguments 
together constitute a ‘group’. The values and M 2 are determined by the criterion 
that the average annual number of such main shocks is greater than or equal to the 
predetermined numbers and ^2 respectively. These values of and n 2 are 
consistently adopted for all squares within a region and for all regions, in order to 
normalize the seismic activity. In defining the spatial concentration, the main shocks 




to be considered, lie within the magnitude range to (Mq — 0-2) while for clustering 
they lie in the range (Mq —2-0) to (Mq —0-1). The flow of main shocks is thus 
normalized with respect to Mq . The time window for which the functions are evaluated 
is 6 years for the first three functions and one year for clustering. The sliding step is 
chosen appropriately, say 6 months. 

A function is identified as being abnormally large when its value exceeds the 
prescribed value, say the 10% percentile of the range of values of the function. Pattern 
M8, is said to have occurred if at two successive instants of time K — 1 and K, at 
least a prescribed number H of individual functions attain abnormally large values, 
and there are at least a prescribed number G of groups which have at least 
one member that attains an abno-rmally large value. A TIP is then declared at the 
time instant K. The duration of TIP is a suitable multiple (say 10) of the sliding time 
step. The free parameters of algorithm M8 were determined by Keilis-Borok and 
Kosobokov (1984) by establishing TIP diagnosis in respect of the strongest earth¬ 
quakes of magnitude 8 or greater that occurred between 1970 and 1983 in 
different parts of the earth. The algorithm was later successfully applied to other 
regions as well, with minimal readaptation of parameters; only the magnitude 
threshold Mq had to be adjusted for every region (Gabrielov et al 1986). For example, 
Mq was 8 for Mexico and 6 for Western Canada. The success-to-failure score was 
no worse than in the original data-fitting (Keilis-Borok and Kosobokov 1984). Similar 
results were obtained in tests for algorithm CN (Gabrielov et al 1986) with a very 
low value of Mq equal to 4*5. 

However, as this test is retrospective, it cannot be considered to be conclusive. But 
it does support the concept of global self-similarity of earthquake flow (Sadovsky 
et al 1982; Gabrielov et al 1986) and warrants similar tests to be made for other 
regions. 


3. The region 

The region considered for this study is the Himalayan sector (figure 1) of the 
Alpide-Himalayan fold belt, including the Indus-Ganga-Brahmaputra foredeep to the 
south of the ranges and part of the Burmese arc. The evolution of the Himalaya is 
traced to the collision of the Indian plate with the Eurasian plate. Its present-day 
seismicity is attributed to continued convergence of the Indian plate towards the 
north by underthrusting beneath the Tibetan plateau. Owing to the gentle dip of this 
underthrusting plane which leads to a wide zone of frictional contact, the seismic 
activity along the belt is quite diffuse and mostly confined within a depth of about 
20 km. At its eastern extremity however, along the Arakan-Yoma fold belt, hypocentral 
depths range from shallow to intermediate (about 200 km). Although major earth¬ 
quakes periodically occur all along the Himalayan collision boundary to relieve the 
accumulating strains, some parts of the belt show more prolific seismicity, notable 
among which is the eastern segment. 

Four devastating earthquakes of magnitude greater than 8-4 have ruptured the 
various segments of the. 2400 km long Himalayan arc during the past 100 years, 
leaving seismic gaps of different lengths in between. Unfortunately, our knowledge 
about the seismicity of the Himalaya in preceding centuries is very sketchy although 
no earthquakes of comparable magnitude are known to have occurred since the 



sixieenin century, it is inereiore aiiiicuit to estimate me recurrence perioa or g 
earthquakes in the region, although earthquakes of magnitudes less than 5 are q 
frequent, interspersed occasionally by earthquakes of magnitude 5 to 6. Since 1 
however, the entire region seems to have entered a period of quiescence as fa 
earthquakes of magnitude 7 or greater are concerned. 


4. The data 

The earthquake catalogue selected for the present study was taken from the NC 
World Hypocenter Data File issue of 1987. More detailed catalogues based on f 
networks are still restricted to small areas and limited time-intervals. A more comp 
global catalogue is of course produced by the International Seismological Centr 
Newbury, England, but the long delay in dissemination of this catalogue mak( 
unsuitable for this study. 

The main shocks were isolated from the aftershocks by the algorithm describee 
Keilis-Borok et al (1980) using the space and time windows given in table 1. 

Epicentral locations of the main shocks are given in tables 2 and 3. The terrii 
lying between latitudes 21-0°N and 36-5°N and longitudes 65-6°E-99-0°E was care! 
scanned (figure 1). It is framed by the Hindu Kush mountains in the northwest farr 
for their strong seismicity at intermediate depths, and by the parameridional Sich’ 
seismic belt lying between 102°E and 105°E. 

The temporal chart of the main shocks given in table 4 depicts a magnitude-^ 
distribution of earthquakes in the various rectangular areas bounded by latiti 
17°N-45°N and longitudes 60°E-105°E. It may be recalled that algorithm M8 
only be applied to diagnose TIPs for the period which starts 7 years after the 
availability of a sufficiently complete catalogue, as some of the functions are evalus 
from data within a time window of six years prior to the instant of their evaluat 
Also, the data pertaining to the first 6 years are necessary to establish as to wl 
values of the functions are abnormally large. Table 4 shows that the catalogu 
relatively complete only for magnitudes greater than 4-5-5-0, and that too after 1 


Table 1- Windows for the identification of after¬ 
shocks (reproduced from Keilis-Borok et al 1980). 


M 

R{km) 

T (days) 

2-5 

19-5 

6-0 

3-0 

22-5 

11-5 

3-5 

260 

22-0 

4-0 

300 

42-0 

4-5 

35-0 

83-0 

5-0 

40-0 

155-0 

5-5 

47-0 

290-0 

6-0 

54-0 

510-0 

6-5 

610 

790-0 

7-0 

70-0 

915-0 

7-5 

81-0 

960-0 

8-0 

94-0 

985-0 


Table 2. Magnitude of the largest earthquake in each of the zones defined by the latitude and longitude. 
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Table 4. Number of earthquakes in the magnitude ranges and time windows in the region defined by latitude H-O^N to 45-0°N and longitude 60-0°E tc 
105 0 E. The numbers in the first column correspond to earthc|uakes with unknown magnitudes. 





5. Selection of Mq 

Four earthquake with M greater than 8 and a larger number with magnitudes ranging 
between 7 and 8 have occurred during the past 100 years in the Himalayan region, 
which is the area chosen for the present study (table 4). Considering the 2400 km 
length of the Himalayan collision zone, an equal number of such earthquakes could 
possibly occur in the future in other segments of the Himalayan belt, particularly in 
the region flanked by the rupture zones of these 4 earthquakes. However, since about 
1970, the region has been relatively quiet. During 1970-1987, only one earthquake 
of magnitude 7-0 or a little higher has taken place. One can, therefore, test algorithm 
M8 only for Mq greater than or equal to 7. For larger Mq, only false alarms, but 
not the failure to predict, could be considered. 


6. Diagnosis of TIPs 

The area under study was spanned by 18 squares centred at points shown in figure 3. 
For Mq = 7, the length (L) of the side of squares turns out to be 560 km, according 
to the formula stated earlier. The lower magnitude cut-off in the catalogue is, however, 
too high (5-0). Therefore, the average annual number of the main shocks is well below 
the required number 20. In some of the squares, this figure is as small as 3. Therefore, 
values of the threshold n had to be lowered. For the common values in all squares, 
we assumed — S and nj = 1, instead of 20 and 10, normally used in earlier studies. 
The other possible alternative of using ^2 = 2 was also tried. 

The threshold values of G and H chosen for recognizing the instant of TIPs, were 
4 and 6 respectively, the same as used in other parts of the world where M8 has been 
applied. Application of this algorithm to the Himalayan belt yeilded TIPs preceding 
strong earthquakes, in the square covering the India-Burma border region, as well 
as some false alarms in others. These results are shown in figure 3. They are reasonably 
acceptable and practically the same for both values of n 2 - The TIPs occupy 18% of 
the total time-space, and precede the only strong earthquake {M^ = 7T) that occurred 
on 29 May 1976, in the India-Burma border region (24-53°N 98-7rE). A TIP is also 
recognized in the Quetta region proceeding a strong earthquake of magnitude 7-0 
which occurred there (33-025°N, 66-315°E) on 3 October 1975. Some agencies reported 
this earthquake as having a magnitude of 6-7. 


7. Discussion of results 

In applying M8 to the present study, the criterion used for recognition of TIPs was 
the same as for other regions of the world. However, the parameters and n 2 were 
adjusted in the manner explained above for the purpose of defining magnitude 




Table 5. Number of earthquakes in the magnitude ranges and time windows in the Quetta region defined by latitude 29-0°N to 310°N and longitude 
6I-0°E to 63-0°E. The numbers in the first column correspond to earthquakes with unknown magnitudes. 




Figure 3. Big dots on the map show the centres of the overlapping squares studied i 
M8 algorithm. Darkened segments along the time axis corresponding to each of the 
indicate TIPs in that square. TIPs terminated by an arrow indicate occurrence of the pred 
event. 


thresholds. This enforced change in the values of and n 2 , makes it difficul 
directly support our results by those of earlier studies in other regions. On the o 
hand, this change hardly made the diagnosis easier. When n is small, the ranc 
variations of all functions are expected to increase, thereby increasing the err 
However, our experience in carrying out this exercise on Himalayan data did 
bring out such randomness. This, though contrary to expectation, is neverthe 
encouraging. For, this may be the result of a higher cut-off for lower magnit 
thresholds of detectability and a near completeness of the catalogue. The ana] 
was repeated by adopting 1967 as the cut-off data as the catalogue of lower magnil 



earthquake that occurred after that date was expected to be more complete. This 
exercise however, yielded only marginally better results. 

These results indicate that future monitoring of TIPs in the Himalayan region, 
using this algorithm may prove to be quite illuminating. It may however be necessary 
to develop region-specific parameter thresholds for various segments of the Himalayan 
belt. To do this, the available seismic data recorded by the national network of stations 
over the past half century need to be rigorously studied to produce a more 
representative and complete catalogue by including even such events which are not 
so well located. 
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rhe importance of strong motion seismology in India 
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Abstract. The high seismicity of portions of the Indian peninsula, together with the high 
density of population and industrial growth, results in a significant seismic risk in many parts 
of the subcontinent. Large construction projects throughout the peninsula require an 
adequate basis for earthquake-resistant design. Thus, as well as strong scientific arguments, 
there are major practical reasons why a substantial programme to record strong seismic 
ground motion should be carried out in India. This paper first reviews the history of strong 
motion instrumental recording, beginning with the important accelerograms obtained in the 
Koyna earthquake of 11 December 1969 through the recent increase in strong motion 
instrumentation, particularly in association with construction of large dams. It is argued that 
there is a pressing need for further extension of strong motion accelerograph coverage of 
India, especially along the seismically active regional thrust faults of the Himalayan region. 
Such programme expansion should follow deliberate strategies of site selection, designed to 
optimize the scientific and practical returns, given the requirements of minimum costs, 
reliable maintenance and accessible data. 

Keywords. Seismicity; strong motion seismology; strong motion instrumentation; ground 
motion recording. 


. Introduction 

rhere are three main reasons why a programme aimed at obtaining widespread 
ecordings of strong earthquakes in India is necessary. First, the large and growing 
>opulation and rapid construction of major engineered structures demand mitigation 
>f the high earthquake risk; secondly, such earthquake measurements contribute to 
he basic understanding of the complex tectonics of northern India; and thirdly, the 
inique continent-continent collision plate boundary of the Himalayas provides a 
aluable opportunity to further understanding of earthquake source physics and the 
leneration of strong ground motion by thrust faults. 

The first two of these propositions have been discussed by Krishna (1981) who 
tointed out that because Peninsular India contains some of the most seismically 
ctive areas of the world there are both broad scientific and important social reasons 
Dr adequate measurements of strong ground motions. As illustrated in figure 1, India 
as had many strong earthquakes since the last century, six of them having a 
lagnitude greater than 8, and 55 with a magnitude between 6-5 and 8-0. Most but not 
11 of these damaging earthquakes have been from strain release in the Himalaya 
egion, extending from Assam in the east to Kashmir in the north, with most of the 

aicmir* />nf>rorv takincr in nnrtViPQct TnHifl Ovf'rflll IlAS hpfin 3 Considerable 




Figure 1. Plot of some historical India earthquakes (adopted from Kaila and Sarkar 1! 

bridges through unequal soil settlement. Because of population and construct 
increase, the next major earthquake could have even greater damage and loss of 
than any previous Indian earthquake. Yet, due to a lack of instrumental observati 
of strong ground motion in India to date, no suitable attenuation laws for h 
amplitude seismic waves are available (Khattri et al 1983) and no reliable estimates 
local site effects. 

There are already available comprehensive earthquake catalogs (Bapat etal 
and maps of earthquake epicenters (Rao and Rao 1984) as well as seismic inten 
zoning maps and seismic risk maps of India (Kaila and Sarkar 1978; Kaila and 1 
1979). Some studies of major fault and seismicity patterns and relevant rec 
tectonics have also been published (e.g. Chandra 1977; Molnar and Pandey, 
Khattri et al all in this volume). This material provides a sound basis for planning 
extended strong motion instrument network with accompanying priorities 
probabilities of observational success. 


2. Seismicity 


Earthouakes occur in manv narts of India with verv variable recurrent rates (Rao 




ao 1984; Gupta et al 1986). Along the Himalayan thrust belt and surrounding areas 
ley are felt almost daily. The high seismicity is exemplified by the twenty earthquakes 
ith magnitudes greater than 5-8 and ten with magnitudes greater than 6-5 that 
:curred along the southern and eastern boundary of the Shillong Plateau, east of 
auki and north and northeast of Halflong during the last 70 years. 

The Dauki-Halflong fault zone (see figure 1), along the straight east-west portion of 
le Assam-Bangladesh border, has shown medium local seismicity in the last 70 years, 
rming a seismic gap between adjacent regions with larger recent earthquakes. In 
rms of the usual frequency-magnitude relation, this fault zone gives an a value of 3-8 
0-year period value) and a b value of 0-45. The low 6-value suggests the existence of 
gh shear stress in the region. 

In the meizoseismal area of the Great Assam earthquake of 1897, the probability in 
) years for ground accelerations greater than 0-2 (M > 6-5) is very high with an 
timated value of about 0-8. For this reason, a strong motion earthquake array was 
•oposed for this area at the International Workshop on Strong Motion Earthquake 
istrument Arrays held in Honolulu, Hawaii in 1978 (Iwan 1978). 

To the south, India is also subject to damaging earthquakes in some regions (see 
;ure 1). The earlier view that most of the Peninsula is aseismic has had to be revised, 
idia differs in this respect from most other large shield areas which have low seismic 
sk. Examples of recent damaging earthquakes are Koyna (10 December 1967), 
^adrachalam (13 April 1969), Broach (23 March 1970), Hyderabad (30 June 1983) 
id Dharamsala (26 April, 1986). One estimate of the h-value for Peninsula India is 
= 0-85 (with a — 4-4), which is typical of moderate seismic region (Rao and Rao 1984). 
istorical earthquakes in the peninsula are concentrated (1) to the northwest from 
)out Koyna through Surat to Kutch; (2) to the eastern coastal zone from about 
adras to Nagaran, and (3) in the southernmost region from about Travancore to 
mgalore. 


Present strong motion instrument network 

ne of the motivations for strong motion recorders in India was the construction of 
rge dams in India (Rau and Varma 1981). The first strong motion instruments of 
odern type were installed as early as 1963 at Koyna Dam in Maharashtra. Small 
rthquakes were located after the reservoir started filling in 1962. This part of India 
id been considered to be nearly non-seismic, and the occurrence of the damaging 
rthquake (magnitude 7) near Koyna Dam in 1967 was a surprise (Tandon and 
loudhury 1968; Chopra and Chakrabarti 1973). 

This experience stimulated various plans for an India National Network consisting 
about 100 strong motion accelerographs to be placed in appropriate sites to cover 
e entire country. The concept included a supplementation of locally built 
celerometers called RESA by 500 simplified instruments (termed SSR and similar in 
sign to the U.S.G.S. seismoscope). As of August 1983, 20 RES As and 200 SSRs had 
en installed in the field. These instruments were designed, fabricated and 
lintained by the Department of Earthquake Engineering at the University of 
)orkee. 

In addition, several river valley project sites in India have installed strong motion 
;trumentation both in the free field and at various elevations in the engineered 



Ramganga, Yamuna, Ukai, Kadana, Koyna and Iddiki (Rau and Varma 1981). In 
1987, about 28 such accelerographs were reported operational at various dams 
(Chandrasekaran 1987). 

In response to the 1978 recommendation from the Hawaii International Workshop, 
a joint Indian-U.S. project of 50 strong motion instrument stations, originally 
recommended as a large-scale strong motion array in the Assam region, was 
established in the Kangra region of the state of Himachal Pradesh (see figure 1). In 
1906, this area experienced an 8-b magnitude earthquake. Forty-six accelerometers 
were installed by August 1983. The Kangra array is roughly in the form of a 
rectangular grid, with principal north-west to south-east axis, spanning 1-5° in latitude 
and 2° in longitude. The array covers the main fault trends in this region. The 
elevation varies between 500 and 2500 meters. 

Somewhat later, a strong motion network at 45 sites called the Shillong strong 
motion array (see figure 1) began to be installed in Assam, also under the auspices of 
the University of Roorkee. This array covered the area between 25°N to 26°N latitude 
and 91°E to 93-5°E with station spacing about 30 km. This rather remote region 
presents considerable logistic difficulties in access and services. The installation and 
maintenance of the Shillong array reflects much credit and dedication on the 
sponsoring organizations. 

Other groups in India are also active with research based on strong motion 
recording. For example, at the National Geophysical Research Institute at 
Hyderabad, digital accelerometers have been deployed in seismically active areas and 
appropriate digital data‘processing algorithms are being developed. 

Parenthetically, we mention that the planned high-level bridge across the Jamuna 
(Brahmaputra) River in Bangladesh is being designed incorporating seismic resistance 
against substantial earthquakes. It is anticipated that part of the World Bank-funded 
project will be the installation of strong motion instruments in the East Bengal 
region. 


4. Recordings of ground motion 

Up to 1986, the only strong ground motions recorded by Indian accelerographs were 
those near the Koyna Dam 230 km southeast of Bombay on the Pre-Cambrian Shield, 
and no time histories of Himalayan earthquakes were available, although 
seismoscope records had been obtained for about 5 events in that region. The 
accelerograms recorded in the Koyna earthquake on December 11, 1967, with the 
peak horizontal acceleration of 0-63 g, represented an important class of recorded 
ground motions. These accelerations were the highest directly measured in the world 
up to that time. 

The persistent interest of Indian earthquake engineers and seismologists began to 
bear fruit in 1986 with clear analog recordings of strong ground acceleration in 
earthquakes in the northern seismic region by a number of stations of the Kangra 
array. Figure 2 reproduces an accelerogram and velocity and displacement ground 
motions recorded in the Himachal Pradesh earthquake, April 26, 1986. The 
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Figure 2. Recorded ground motions (N14E component) at the Dharamsala station in the 
Himachal Pradesh earthquake, April 26, 1986 (courtesy: A R Chandrasekaran). Epicentral 
distance about 5 km. 

publication of the set of these records with peak accelerations in excess of 0-2 g will be 
a valuable addition to the strong motion data base of India. 

On August 20, 1988, a major earthquake (Ms = 6-8) occurred in the Bihar-Nepal 
region, in approximately the same place as the devastating earthquake of January ,15, 
1934 (Ms = 8-4). The recent earthquake had high intensity over a substantial area 
(Richter 1958), with 800 loss of life from collapsed structures reported. It underscores 
the need for adequate strong motion instrumentation along the whole Himalaya 
region (Arya et al 1986). 

5. Future developments 

The Department of Science and Technology of the Government of India has reported 
an integrated project in Himalayan seismicity, particularly in the Himachal and 
Garhwal areas in the Eastern Himalayas covering Assam and Meghalaya. As well, 
monitoring of seismic activity is carried out by the India Meteorological Department 
and by other organizations such as the Department of Earthquake Engineering in 
Roorkee and the National Geophysical Research Institute. 

Construction of large dams in various River Valley Projects in India has led to 
additional seismic surveillance in various parts of the country (Srivastava and Dube 
1982). The minimum seismographic networks needed for engineering projects of 
various sizes have been discussed (Srivastava 1987), both of sensitive instruments for 



1 he impounding oi bhivaji bagar Lake m ly62 by Koyna Dam led to much interest 
in reservoir-induced seismicity in India (Gupta, 1985). The continued occurrence of 
earthquakes in its vicinity (Rastogi and Taiwan! 1980) as well as the new major dams 
such as Navagam Dam in Gujarat provides ample justification for adequate future 
strong motion as well as sensitive seismic network around such facilities. 

As mentioned in the introduction, the development of an enhanced all-India strong 
motion instrumentation network provides a singular opportunity to improve the 
basic knowledge of India earthquakes and to provide fundamental data for aseismic 
engineering design and seismic risk mitigation (Agbabian and Arya 1987). 
Continental-wide considerations and integration of both seismological and 
engineering needs (Bolt 1980) will help ensure that the observations of strong ground 
shaking are optimized and representative of various geological and soil conditions. 

Because of the unique continent-continent collision tectonics of the Himalayan 
thrust, strong motion recordings in this region will also provide important 
information for fundamental studies of earthquake mechanisms, not only in India, but 
in similar geological conditions around the world. For example, recent research 
suggests that thrust type earthquakes produce higher accelerations than the normal 
and strike-slip earthquakes. This prediction depends on the friction-rupture 
relationships for faulting; other reasonable models might predict lower accelerations 
for thrust faults. Data from the unique conditions of the Himalayan thrust should 
provide important constraints on the parameters for modelling major earthquake 
ruptures. 

Seismic gap theory should be given strong weight in consideration of the 
geographic strategy for deploying Indian strong motion arrays. It has now become 
clear from numerous successful predictions worldwide that seismic gaps usually 
represent regions of highest probability for very large earthquakes, although 
aftershock zones of earlier earthquakes often are the zones of highest short-term 
probability for smaller events. As the historic pattern of the occurrence of large events 
is now much better understood in relation to seismic gap theory (see other papers in 
this volume), it should play an important role in siting strategy. 

It is hardly necessary to stress the need for a national plan aimed at standardizing 
the instrumentation and vital maintenance, developing rapid and validated recording, 
processing facilities, and establishing efficient data dissemination procedures. 
Valuable recommendations, objectives and resolutions relating to these questions are 
listed in the Proceedings of the 1987 Indo-U.S. Workshop on earthquake disaster 
mitigation research (Agbabian and Arya 1987). 

Many decisions are required concerning the relative distribution of free-field 
instruments in relation to specific problems. Such instruments provide in time a sound 
physical basis for understanding the dynamics of faulting, the transmission of seismic 
waves through rocks and soils, possible nonlinear effects at soft-soil sites, topographic 
effects and the characteristic response of alluvial basins. Instrumentation of 
representative engineering structures such as buildings and dams is, of course, 
essential as well as installations to provide the seismic response of such “lifelines” as 
bridges and piping. 

The major projects involving the Kangra and Shillong arrays are now providing 



experience in India with large regional arrays of the type classified as “source 
mechanism and were propagation arrays” at the 1978 International Workshop report 
(Iwan 1978). These Indian arrays are well situated to measure ground motions from 
great earthquakes that may occur from fault slip along two segments of the major 
boundary thrust fault system. Consideration of the priorities for arrays of different 
classification (Bolt et al 1987) such as liquefaction and soil-structural arrays is now 
opportune. 

Objectives and guidelines for the deployment of strong ground motion arrays 
require appropriate indigenous studies but, in general, two criteria can be 
recommended: 

(i) Arrays should be sited with specific aims to measure well-defined aspects of strong 
seismic motion, and 

(ii) Siting should follow studies that optimize the scientific and engineering returns, 
given knowledge of the local three-dimensional geologic environment and the 
distribution of earthquake sources. 

In the last few years, considerable experience has been obtained in a number of 
countries, particularly the United States (California), Japan and Taiwan, with digitally 
recording instruments. The older analog instruments were thought (e.g. Iwan 1978) to 
provide more reliable and tested data acquisition in various harsh field conditions. 
The modern digital accelerometers, however, have proven to be comparable in 
reliability and because of their extended dynamic range and simplicity of data 
processing (now with regular PC tape or disc input), they must be seriously 
considered. A decision in each case depends on answering the questions: What is the 
expected maintenance of each system? And how will limitations on operational life 
affect the overall recording objectives in the region? 

Because large earthquakes may not occur for several years, there are special 
organizational and funding problems in maintaining long-range programmes for 
instrument maintenance, data analysis, archiving, distribution and data use. The only 
solution is stable, continuous commitment of government agencies and organizations 
with vigorous support from professional groups of seismologists and earthquake 
engineers. It is sometimes possible, as in the current California Strong Motion 
Instrumentation Programme (SMIP) to obtain significant financial support from 
mandatory small fees from large-scale building owners that are assessed at the time of 
building permit approval. (Currently this fee in California is 17 cents per 100 dollars 
construction cost, or 0-017%). 

It is crucial to pool strong motion data from every seismically active country of the 
world in a global effort to understand strong earthquake motions and mitigate seismic 
risk. An upgraded and integrated India strong motion network could contribute 
much to this international objective. The opportunity is particularly inviting 
following the recent creation of the Decade of Natural Disaster Reduction (1991- 
2000) by the United Nations. For a solution to the problems of earthquake hazard 
reduction, a key ingredient is more adequate sampling of strong seismic motions. 
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